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IMPROVED ESTIMATES OF THE AREA DENSITY 
OF CLOUD-TO-GROUND LIGHTNING OVER 
THE UNITED STATES 


Michael W. Maier 

Lightning Location and Protection, Inc. 
Tucson, Arizona 

and 

Jackie M. Piotrowicz 
NOAA/NHRL 

Coral Gables, Florida 


ABSTRACT 

A key element for any successful and cost-effective lightning protection 
design is an accurate estimate of the frequency of direct lightning strikes to 
ground. For example, an effective protection system for an electric power 
transmission line in Florida (an area of high lightning incidence) would be 
far too costly and not necessary for an area such as California (an area of 
low lightning incidence). In order to estimate the frequency of lightning 
strikes the local area density of lightning flashes to ground, i.e. the num¬ 
ber of strikes per unit area per month or per year, must be known. In the USA, 
ground flash density is often assumed to be a relatively simple function of the 
mean number of thunderstorm days per month or per year. Many design handbooks 
give formulas which relate ground flash density to thunderstorm days, but no 
formula is generally accepted as standard. All too often, the actual lightning 
strike incidence is considerably different from the lightning strike projections 
based on the thunderstorm day statistics. The resulting over or under design 
of lightning protection systems ultimately leads to greater costs. 

Here we show that the thunderstorm duration statistics are a much better 
indicator of true ground flash densities than the number of thunderstorm days. 

Maps of mean annual ground flash density are presented which have been inferred 
from long term thunderstorm duration observations. The thunderstorm duration 
data were obtained by the relatively dense network of existing aviation weather 
reporting stations. Approximately 450 stations are used, each with an uninter¬ 
rupted ' ,0 year period of record. The station density available here is essential¬ 
ly twice that of any previous thunderstorm frequency analysis of the United States. 

Our analyses show that maximum annual ground flash densities of 18/km 2 are 
found in the western interior of Florida from Naples to Cross City. High flash 
density values (> 12/km 2 ) are found over much of the Florida peninsula, the Flor¬ 
ida panhandle, and westward through the southern gulf states to eastern Texas. 

A large area of flash densities > 8/km 2 are found over most of Oklahoma, Kansas, 
Missouri, Arkansas, Louisiana, Mississippi, and Tennessee. High flash densities 
(> 8/km 2 ) in the western states are limited to south central Colorado, north¬ 
eastern New Mexico, northern Arizona and extreme southeastern Arizona. Pro¬ 
nounced flash density minima are found over the southern half of Texas and in 
the southeastern states from central Alabama northeastward through western 
Virginia. Most of the western and northeastern states have flash densities 
< 4/km 2 . Maps of the 95 0 i and 5" probability flash density values are also pre¬ 
sented. Verification of these maps awaits a nationwide direct measurement pro¬ 
gram or successful practical application. 


This paper was not available for incorporation into this book. 
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Lightning and Transient Protection of Radio Sys 


By 


Prof. S. Lundquist 
Institute of High Voltage Research 
Uppsala University, S-755 90 Uppsala, 
Sweden 


ABSTRACT 


Protection of important radio systems against electromagnetic transients 
lightning discharges and nuclear explosions (NEMP) nowadays is more or less a 
standard procedure. However, the problem of coordinating different protective 
measures and selection of the most suitable components still are problems dese 
further studies. Tests have therefore been performed to evaluate some importai 
parameters of individual components and of radio systems for frequencies up to 
75 MHz. At the testing both standard test pulses as well as fast rising pulses 
the nanosecond range has been used. Gas discharge tubes, cables etc. have bee" 
tested separately and as part of a system. In this way it is possible to study 
effect of different counterpoise and grounding systems, different cable lengths 
and different locations of the cables as well as the location and properties o r 
discharge tubes and semiconductor protection elements on the amplitude and ene: 
of the transient reaching the sensitive parts of the system. The results are . 
used to estimate the necessary intervals for routine maintenance and to predic 
MTBF for faults caused by lightning. A technique for field tests in actual in: 
lations will also be discussed. 








Coupling of Natural Atmospheric Interference Through Apertures 
By 

J. Hamm, W, Graf, and E. F. Vance 
SRI International 
Menlo Park, CA 94025 


ABSTRACT 

The coupling of electromagnetic energy through apertures has received a great 
deal of attention in recent years. This recent research has been both theoretical 
and experimental in nature and has been focussed on the problems associated with 
apertures which are small compared to a wavelength. 

The results of this research can be applied to the problem of determining the 
transient signals induced on conductors behind an aperture by lightning and static 
charging. While the problem of determining induced transients in complex systems 
(such as aircraft or ground-based facilities) cannot be solved exactly, the 
research of the past few years indicates that upper bounds can be determined for 
these transients. This paper reviews our recent work on aperture coupling and 
applies the results to models of airborne and ground-based facilities. 

The nature of this recent work is such that reasonable upper bounds can be 
determined without the use of extensive computer codes. The system designer can 
determine quickly whether or not a particular aperture presents a problem for elec¬ 
tronic systems located behind the aperture. The analysis performed in recent years 
has been of two different types: transmission line modeling of canonical circuits 
behind the aperture, and matrix formulation of Maxwell's equations. Both approches 
lead to relatively straightforward expressions for bounds on transient voltages 
and currents on conductors behind the aperture. Our laboratory measurements on 
test apertures are used to investigate the suitability of these bounds. 

The results of this research can also be applied to the problem of generating 
realistic specifications, standards, and test procedures for certifying systems 
which must operate in the natural electrical environment. 
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ATMOSPHERIC ELECTRICITY INTERACTS with aircraft 
to produce transient electric and magnetic 
fields at the outer surface of the aircraft. 
These surface fields and the signals Induced on 
wire penetrations of the surface act as 
equivalent sources for the transient electro¬ 
magnetic environment inside the aircraft. A 
companion paper 1 describes this problem and 
shows that wire penetrations of the surface are 
usually the most significant contributors to the 
internal transient environment. If all wire 
penetrations of an aircraft surface have been 
properly treated, the apertures on the exterior 
surface of the barrier can become significant 
points of entry for electromagnetic inter¬ 
ference. These apertures may be in the form of 
holes or openings on the surface or in the form 
of seams and joints. 

To solve the aperture problem, the system 
designer must: (1) determine quantitatively the 
effects of the apertures on the internal 
voltages and currents, and (2) develop measures 
for "closing" those apertures that affect system 
performance. 

THEORETICAL BACKGROUND 

The types of apertures encountered with 
typical aircraft are holes with maximum 
dimensions of a few meters, long narrow slots, 
structural seams and access hatches. For 
analysis purposes, these apertures can be 
considered "small;" that is, the characteristic 
dimension is smaller than a wavelength. Long 
narrow slots and seams are considered as 
distributed small apertures, and their coupling 
properties can be expressed per unit length. 
With the discussion being limited to "small" 
apertures as defined above, the following 
discussion of the theoretical background is 
generally valid for frequencies below 100 MHz. 
This is not a severe restriction, since the 
spectral content and coupling transfer functions 
appropriate for lightning are usually not 
important above 100 MHz. 

Two recent publications 2 ' 3 present excel¬ 
lent reviews of the literature describing 
aperture coupling. These problems ae difficult 
to analyze, especially the case of an aperture 
on a cavity containing equipment and wiring. 
The only configuration for which analytically 
rigorous results are available is that of a 
small circular aperture on a conducting plane of 
infinite extent. Numerically "accurate" results 
are available -for some cavity and aperture 
configurations. These were obtained by a 
numerical solution of integral equations in 
terms of the unknown aperture field compo¬ 
nents. The aperture fields were then considered 
as sources for the fields existing behind the 

The mechanisms of field leakage through an 
aperture are described in Figures 1 and 2. The 


normal electric field and the tangential 
magnetic field couple through the aperture to 
interact with conductors and wiring loops behind 
the plane. As the figures imply, this coupling 
problem is equivalent to the problem of electric 
and magnetic dipole sources located at the 
aperture with the aperture shorted. These 
sources, radiating in the presence of the 
conducting plane, produce the same fields in the 
shadow region as do the original external 
sources. Formally, the problem is reduced to 
finding the dipole moments of the equivalent 
sources and computing the fields produced by 
these sources. 

Figure 3 shows the general aperture 
coupling problem of interest for electromagnetic 
interference control. An external electro¬ 
magnetic field is assumed to be incident on a 
volume containing electronic equipment and 
wiring. The volume is assumed to be enclosed by 
a metallic surface containing holes and seams. 
The interaction of the incident field with the 
outer surface produces a surface surrent density 
J and a normal electric field component E n at 
tne apertures. These surface fields couple 
through the aperture to produce transients on 
conductors inside. The internal transients may 
be in the form of (1) surface current densities 
on equipment enclosures or (2) voltages and 
currents on exposed equipment wiring. The 
internal surface current densities Jj can 
interact through apertures on the enclosures 
with conductors at the next topological level, 
and the voltages and currents can drive these 
inner levels directly. 

The practical problems of interest usually 
concern volumes that are rather densely packed 
with equipment and wiring, where the wires are 
often routed in harnessed bundles. These 
coupling problems cannot be solved rigorously 
because of the complexity of the internal 
configuration of wires and conducting 
surfaces. Recent research into this problem has 
concentrated on determining upper bounds on the 
internal transients, and significant progress 
has been made in the past few years. Davis 5 has 
Investigated bounds on coupling through an 
aperture to a single-wire circuit loaded at each 
end, and other work has addressed the more 
important multiconductor problem. ,b Both 
problems are formulated In the same way. The 
aperture excitation is modeled by the usual 
electric and magnetic dipoles radiating at the 
aperture location with the aperture shorted, 
l.e. closed. The single-wire and multiconductor 
circuits are modeled as transmission lines above 
a ground plane, with impedance loads at each 
end. The transmission lines are driven by the 
fields of the equivalent dipole sources, and 
this excitation is modeled by Inserting 
equivalent voltage and current sources in the 
transmission line. Figure 4 describes this 
model, a key requirement of which is that the 


*Numbers in parentheses designate references 
listed at the end of the paper. 


$ 

£>; 


18-2 











circuit be separated from the aperture by at 
least a few aperture radii. 

The dipole moments can be bounded by 
bounding the aperture polarizabilities, and the 
induced currents and voltages can be bounded by 
applying the usual transmission line equations 
to the receiving circuit. For the single-wire 
receiving circuit, the mathematics is 
straightforward. Analysis of the multiconductor 
circuit requires a matrix formulation to 
generate bounds. The single-wire case is a 
special case of the more general matrix 
formulation. 

The aperture coupling experiments described 
in this paper have investigated coupling to 
circuits very near the aperture, as well as 
circuits distant from the aperture. The results 
of these experiments are used to verify recent 
theoretical models. 

MEASUREMENT CONFIGURATION 

A metal box was constructed with dimensions 
of 1.22 m by 1.83 m by 0.91 m. The box was 
excited by a parallel plate transmission line 
with a source and a matched load at the two ends 
of the line, as shown in Figure 5. The upper 
surface of the box was designed so that a 
variety of apertures could be located under the 
upper plate of the transmission line. This 
excitation resulted in a field distribution at 
the aperture as shown in Figure 6. The source 
driving the parallel plate line is shown in 
Figure 7. This pulser circuit supplied an 
incident electric field with a peak value of 18 
kV/m, with a risetime (zero to peak) of about 
10 ns, and a pulse duration of approximately 8 us 
(to 1/e of the peak). Figure 8 shows the 
electric field in the center of the parallel 
plate region. 

APERTURE VOLTAGE MEASUREMENTS 

The theoretical study by Casey 7 indicated 
that a small circular aperture can be considered 
as an inductance, and that the voltage induced 
across the aperture by magnetic field coupling 
is the product of the inductance and the time 
derivative of the external surface current 
intercepted by the aperture. Figure 9 shows the 
external current driven over the aperture and 
the resulting voltage measured across the 
aperture on the Inside of the test volume. The 
table in the figure Indicates the accuracy of 
Casey's quasi-static model for the aperture 
sizes used in the experiment. The waveshape of 
the aperture voltage Is very close to that of 
the time derivative of the external surface 
current. 

Another Important contribution of Casey's 
work is the development of a simple model for 
predicting the effectiveness of a mesh screen 
placed over the aperture. With the untreated 
aperture modeled as an inductance, a wire mesh 
cover can be modeled as a parallel Impedance 
element, as shown in Figure 10, where is 

the effective impedance of the mesh. This 


Impedance can be evaluated using the expressions 
developed by Casey. The table on Figure 10 
describes the effectiveness of various mtsh 
screens In reducing the aperture voltage and 
indicates the accuracy of Casey's model. 

These quasi-static models ar.- useful for a 
quantitative understanding of coupling through 
an aperture and are formally villd only for a 
circular aperture on an inf* .ire ground plane. 
The formal extension to noncircular apertures 
backed by partially filled cavities is a 
difficult problem that has not been solved. 
Certain bounding statements can be made, 
however, that relate to the usefulness of the 
results of the analysis of a circular aperture 
in an infinite plane. Both Casey and Davis 
state that the characteristics of noncircular 
apertures can be bounded from a knowledge of 
coupling through circular apertures, and that 
coupling through an aperture on an infinite 
plane is an upper bound on coupling to an 
aperture in a cavity. 

INTERNAL TRANSIENT MEASUREMENTS 

A single wire was inserted in the test 
volume in the configuration shown in Figure 
11. The dimensions W and D were varied, and the 
wire currents were measured under a variety of 
loading conditions. When the wire was short- 
circuited at each end, the observed current 
waveshape was similar to the waveshape of the 
incident electric and magnetic fields. When the 
wire was terminated in the characteristic 
impedance of the wire-to-surface transmission 
line (approximately 240 Q), the current observed 
on the wire had a waveshape similar to the 
derivative of the incident (external) fields. 
Figure 12 displays these measured currents for 
one set of wire locations. The tables in the 
figure contain the peak amplitudes of the 
currents for two values of the parameter W and 
for a range of values of the parameter D. 

The internal current measurement con¬ 
figurations were used with the analytical 
results of Davis in order to investigate the 
bounding process. 

Table 1 shows that the computed currents 
are indeed larger than the measured currents, 
generally by less than a factor of 10. The 
first entry in the table (D - 0.007, W - 0) 
describes a configuration in which the Internal 
wire is right in the aperture, a case for which 
the theory breaks down; i.e., the wire is well 
within an aperture radius of the aperture. 
However, the computed upper bounds are 
reasonable for the other cases listed for W * 
0. As the wire is moved away from the aperture 
(W = 0.23 m), the agreement Improves somewhat. 

SUMMARY 

This paper has reported on an experimental 
investigation of recent theoretical developments 
concerning aperture coupling. Continued 

development of bounding theories, especially for 
the multiconductor configuration, will allow 
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system designers to assess the Importance of 
apertures to the lightning immunity of air¬ 
craft. These theoretical developments may also 
be reflected In the structure of new electro¬ 
magnetic specifications and standards used for 
aircraft procurement. 
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FIGURE 6 FIELD VECTORS AND SURFACE CURRENT AT APERTURE 
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gths and the fields calculated inside the cab 
icated that the induced currents would 
bably be less than that which would result 
n strikes to the center dividing posts on the 
ishieId of the cab. 

The tests thus indicated no unusual strike- 
nts in this series of tests. 

HIGH CURRENT DAMAGE TESTS - The test 

erator for the high current tests consisted of 
large bank of capacitors charged in series 
allel and discharged in the same configuration 
the vehicle. The screen mat underneath the 
icle was carried up to the base of the current 
erator which was located on a truck bed in 
er to minimize inductance and with this 

angement currents exceeding 200,000 amperes 
Id be fired to Che vehicle. Discharges were 
ed to all the major strike points and no 
ets of the launch system occurred in any of 

test strikepoints with any of the test 
ikepoint locations. 

CLUSIONS 

In a series of tests on the LPC and the MLRS 
icle with the LPC installed no upset or damage 
urred whi:h would cause an inadvertent launch 
h the exception of the early tests before the 
tection bar was installed on Che LPC. Thus the 
tern withstood approximately 150 high voltage 
g arc lightning discharges and 50 discharges 
ween 150,000 and 200,000 amperes. This 
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ABSTRACT 

This paper presents a preliminary analy¬ 
sis of the data concerning triggered light¬ 
ning flashes obtained at the Langmuir Labora¬ 
tory in 1982. Discharge currents, visible 
light emitted by return strokes channel and 
nearby magnetic fields detected at ground are 
correlated. Measurements of distant VHF radi¬ 
ation by means of a spectrum analyser and a 
two dimensional interferometer location sys¬ 
tem at 300 MHz are discussed. 

DURING THE SUMMER OF 1982, an experiment on 
triggered lightning was held at the Langmuir 
Laboratory of the New-Mexico Institute of 
Mining and Technology (USA). This experiment 
was conducted by investigators from the US 
Air Force (FDL laboratory - WPAFB), New 
Mexico Tech, ONERA (Office National d’Etude 
et de Recherches A£rospatiales - France) and 
CEA (Commissariat A l'Energie Atomique - 
France). It had two main objectives : 

- studies of the behaviour of a hemispherical 
ended cylinder that simulated an aircraft 
fuselage when strucked by lightning ; 

- analysis of the electrical properties of 
triggered lightning. 

In addition to development and testing 
of the equipment associated with the instru¬ 
mented cylinder, the main experiments, to 
which this article is primarily related, are 
as follows : 

a) current measuring at flash foot ; 

b) detection of the visible light emitted by 
the lower part of return strokes channel ; 

c) measurement of magnetic field at ground 
produce by nearby flashes ; 

d) measurement of distant VHF radiation by 
means of a spectrum analyser and a two dimen¬ 
sional interferometer location system at 300 
MHz. Also, a high time resolved measurement 
of the electromagnetic radiation during a 
storm was performed. 


A preliminary analysis of some of the 
data procured during the summer is conducted 
in order to clarify the origin of VHF radi¬ 
ation during lightning. 

For that purpose we examine the proper¬ 
ties of short pulses observed during a storm 
and then dlscusse emission characteristics 
from triggered lightning as related to flash 
phase and as deduced from various measure¬ 
ments . 

EXPERIMENTS 

A diagram giving the location of all 
measuring equipment is shown in figure 1. 

TRIGGERING SITE - Current is Measured by 
means of a coaxial shunt connected to the 
triggering device. The resistance of this 
shunt is 0.001 ,2 and its rise time is 
shorter than 80 ns. The signal is recorded on 
a wide dynamic range. The bandwidth of the 
current recording equipment is DC to about 1 
MHz. 

Triggering Equipment - Lightning flashes 
are triggered beneath stormy clouds, using 
the well-known wire and rocket technique [1]* 
[2]. A small rocket, of the anti-hail type, 
unspools a steel wire, 0.2 mm diameter. The 
lower-part of the wire is connected to ground 
through the current shunt device. Triggered 
flashes ire obtained when the rocket reaches 
heights between 50 and 500 m depending of the 
electrostatic field strength created by the 
cloud. 

Light Sensor - A light probe is in¬ 
stalled about 300 m away from the triggering 
site ; it is directed toward the lower part 
of the channel that corresponds to the con¬ 
ducting wire which is vaporized at the begin¬ 
ning of the event. With an 8 degrees light 
acceptance angle, the measured optical signal 
♦Numbers in parentheses designate References 
at end of paper. 
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rs about 33 m of channel length. The 
i gives information in the visible range 
to 700 nm). 

Magnetic Sensor at Ground - A magnetic 
Dr is installed about 200 meters from the 
jering site. Its bandwidth is 100 Hz-20 
and the receiving signal is connected to 
ttery operated videotape recorder with a 
z bandwidth 13]. Its dynamic range is 
sted in order to record the effect of 
a current pulse passing through the flash 
nel. 

ELECTROMAGNETIC MEASURING SITE - Exper¬ 
ts are performed at the West-Knoll site 
ted about 1300 m south of the triggering 
(see figure 1). They include the follow- 
neasurements : 

Wide Band Measurement - The instrumental 
ce is made of a wide bande antenna con¬ 
ed to a fast transient digitizer (Tektro- 
7912 AD) which is controlled by a com- 
r (Tektronix 4052). The bandwidth of the 
nna is 20-200 MHz but we have verified 
it can be use to 500 MHz. Concerning the 
anse of the system to impulse signals, we 


enna with the recorder, only the first 10 
can be considered to represent the actual 
pe of the incoming pulse. 

The transient digitizer is used with an 
lifier that has a 600 MHz bandwidth and a 
ectable signal level of 1 mV. The computer 
programmed to record about one event per 

300 MHz Interferometer System [4] - The 


A/2 dipoles, installed at A/4 of a vertJ 
supporting mast. We have verified that 
coupling between antennae is negligible ; 

3 dB bandwidth of each channel is 370 k 
The sensitivity of each channel is ad jus 
following results obtained by other exp 
imental ists [5] : E.T. Pierce gathered 
results concerning different measureraer 
normalized at 10 km for a 1 kHz bandwidth, 
figure 2, we have indicated the threshold 
the interferometer and of the spectrum a 
lyser, in regards of others results wl 
magnitudes are gathered around the follow 
characteristic : 


WIDE-BAND MEASUREMENTS [6] - Th 

measurements are not synchronized with ot 
data and correspond to natural lightning 
well as triggered events. The results 
tained during the 82 campaign, bring out 
following points : 


whose shape and rise time are almost I 
ical regardless of the type of flash ; 

- the rise time is few ns ; 

- the pulses are emitted at a rate betw 
and 20 per ps. 

Figure 3 gives the shape of a ty 
pulse ; the table, shown on figure 4, 
cates the characteristics of 20 diff 
pulses obtained during a single storn 


which is a fraction 


RESULTS CONCERNING 


FLASHES - 
1 ightning 
Triggered 











B : 12h51mm46s, day 


46s, day 8.5.1982 ; event B : 12h51mm46s, day 
8.7.1982). 

Event A - The first radiation detected 
by the interferometer occurs 32 ms before the 
first return stroke and a few ms after the 
vaporization of the wire. Angular location of 
the VHF sources may not be obtained before 
the radiation associated with the first 
stroke, because the signal received is too 
weak to allow for a full reduction of data. 

Figure 5a shows the current versus time 
variation ; figure 5b shows the angular (el¬ 
evation-azimuth) locations of the sources 
corresponding to the flash. Suitable inter¬ 
ferometer information indicates several 
spatial continuous organizations of emitting 
sources. Associated with the first stroke is 
an upward propagation of sources, from the 
location of the triggering site. These move 
with a mean velocity of 4.10^ m.s -1 if 

it is assumed that all the events happen in a 
vertical plane containing the triggering 
site. Other organized locations of sources, 
indicate downward propagations that are prob¬ 
ably associated with dart leaders. These have 
a mean velocity between 10^ and 5.10^ 
m.s -1 . Figure 6 shows, for the first 
stroke and a dart leader, the time evolution 
of the elevation angle of the sources. These 
curves clearly indicate the upwards or down¬ 
wards propagation of the phenomenon. 

Event B - This negative triggered flash 
exhibits 16 large current pulses as recorded 
by the magnetic sensor close to the channel. 
However, during the same time, the optical 
sensor indicates 46 fast luminous pulses. 
After evaporation of the wire, the flash 
strikes the measuring equipment but the lower 
part of the channel is diverted so that we 
only have measurements before vaporization 
and for the four last strokes. The overall 
duration of the flash is 412 ms. 

The first stroke occurs about 12 ms 
after wire vaporization. The magnetic sensor 
at ground indicates a corresponding current 
of 16 kA with a rise time of 15 ps. The maxi¬ 
mum current amplitude during the flash, as 
deduced from this sensor, is that of the last 
stroke : 19 kA, rise time 1 ms. 

The radiation first detected at 300 MHz 
by the interferometer is simultaneous with 
the first stroke ; angular locations associ¬ 
ated with it are scattered above the trig¬ 
gering site. Most of the angular locations 
obtained during this flash are spatially 
scattered ; they correspond to bursts of a 
few hundred ms duration, that are generally 
associated with current or light pulses. On 
two cases, we observe a downwards organized 

propagation which may correspond to two dart 

leaders ; The velocities are lO^m.s -1 

and 8.10° m.s -1 . 

We have drawn on figure 7, current and 
light emitted by the flash versus time ; we 


angular location of three different downwards 
propagation sources defined on figure 7 
(points 0^, O2, 0^). 

VHF radiation is detected by the spec¬ 
trum analyser just before the beginning of 
current flow due to the upwards positive 
leader from the wire. Radiation is gathered 
in burst corresponding to light pulses. One 
of the most general results concerning this 
measurement is that the magnitude of the 
radiated field stays within a short dynamic 
range, whatever the actual component of the 
flash. Nevertheless, the temporal evolution 
of the radiation seems closely related to the 
flash to ground phase ; we have drawn, on 
figure 9, radiation associated with the 
stepped leader before the first stroke O4 
and, on figure 10, the radiation associated 
with the dart leader preceding the last 
stroke. 

The stepped leader radiation begins by a 
bursts of 77 ms duration, 1.9 ms before the 
stroke O3 (fig- 9). It is followed by short 
bursts of 5 to 10 ms width that are rather 
regularly spaced by blanks of 50 to 80 ms. 
The strongest radiation corresponds to the 
return stroke. 

The duration of radiation from this 
first return stroke is 14Q ps and corresponds 
to the duration of current deduced by the 
magnetic sensor. However, the main current 
pulse, as recorded by the current shunt, 
lasts only 40 ms. 

The dart leader before the last return 
stroke exhibits very different VHF character¬ 
istics. As shown in figure 10, radiation is 
more continuous and Increases slightly as the 
radiation sources (300 MHz) approach the 
ground in 600 ms time. All radiation ends at 
the start of the current pulse from the re¬ 
turn stroke. 

•DISCUSSION 

Data reduction of TRIP 82 measurements 
is still, in a preliminary stage and must be 
improved in the future. 

Nevertheless, it is now possible to 
point out some interesting facts concerning 
electromagnetic radiation produced by light¬ 
ning discharges. 

Regarding electromagnetic measurements, 
we are not able to see big differences be¬ 
tween triggered lightning and initially up¬ 
wards directed natural flashes. Anyway, the 
most significative result is that VHF radi¬ 
ation due to an upwards positive leader is 
weak (less than 50 MV/m.kHz at 10 km at 300 
MHz). It Is just detected by the spectrum 
analyser which sensitivity is greater than 
the interferometer one. 

We have observed, during natural or 
triggered lightning, short radiation pulses 
with a typically 2 to 5 ns rise time, and a 














to well known results on lightning (fig. 3) ; 
this strengthen the idea of the Impulsive 
nature of VHF radiation due to lightning. 

Results indicated by the interferometer 
have the same general characteristics as 
those concerning natural flashes observed 
during a previous campaign and presented in 
an other paper [7] : 

- radiation is, most of the time, constituted 
by pulses of less than 1 ms duration ; 

- those pulses are gathered within longer 
radiation bursts of few hundreds ms dur¬ 
ation ; 

- the corresponding locations of the radi¬ 

ation sources are either scattered in a wide 
angular range or gathered, to indicate a one 
way propagation which is associated with a 
mean velocity of few 10^ m.s - * for both 

dart leader and return stroke. These velocity 
values are not in disagreement with results 
obtained by other authors [8] [9]. 

Our 1982 measurements include new comp¬ 
lementary and simultaneous indications ob¬ 
tained with the VHF spectrum analyser. The 
main results are the following : 

-at any given frequency the strength of the 
radiation seems not to be strongly dependent 
on the flash component (leader, return 
stroke. ..) ; 

- the temporal variation of the radiating 
field is typical of the flash component as 
already mentioned by others [10] ; 

- VHF spectrum is analogue to the curve on 
figure 2, but the level obtained is a little 
less than the mean value. 

Comparison of two return strokes indi¬ 
cates that one is associated with propagating 
radiation sources but that these are com¬ 
pletely absent in the second stroke. In both 
cases, there is no VHF radiation after the 
strokes. This strongly suggests that the 
radiation that we observe is not directly 
associated with the current wave of the 
return stroke but with some other mechan- 

The continuous propagations, indicated 
by the interferometer, for certain bursts of 
pulses, show that the emission zone within 
the ms is probably limited to a few m to tens 
of m large. 

CONCLUSIONS 

Measurements of VHF radiation, simul¬ 
taneous with other measurements that identify 
each flash component, are basic for an under¬ 
standing of the radiation associated with 
1 ightnlng. 

A preliminary analysis of the data ob¬ 
tained during the summer of 1982 raises im¬ 
portant questions regarding the origin of the 
radiation recorded. It appears that strong 
VHF radiation is not always associated ex¬ 
clusively with the propagation and location 
of the return stroke. 

Triggered flashes, in particular the 
frequently obtained "anomalous" flashes that 


deviate from the evaporated wire, exhibit be¬ 
haviour that is close to that of naturally 
occuring upward initiated discharges. This 
makes triggered lightning a powerful tool for 
all lightning studies including more techni¬ 
cally oriented effects such as their interac¬ 
tions, with aircraft. 

All these considerations Imply the need 
to Improve the instrumentation in future 
program. A three dimensional, more sensitive, 
interferometer should be used in conjunction 
with a spectrum analyser and a wide-band 
antenna system. 
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LIGHTNING ATTACCENT PATTERNS AND FLIGHT CONDITIONS 
EXPERIENCED BY THE NASA F-106B AIRPLANE 
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During the NASA Langley Research Center 
Storm Hazards Program, thunderstorm penetrations 
have been made during three summers by an exten¬ 
sively instrumented F-106B research airplane in 
order to record direct lightning strike data and 
the associated flight conditions. The penetra¬ 
tions were made in Oklahoma and Virginia in con¬ 
junction with ground-based weather radar measure¬ 
ments by the NQAA-National Severe Storms Labora¬ 
tory and NASA Wallops Flight Center, respec¬ 
tively. In 1980 and 1981, the airplane received 
10 direct lightning strikes each year, with most 
penetrations occurring near the freezing level. 

In 1982, 239 penetrations were made during which 
156 direct lightning strikes occurred. Most of 
the penetrations in 1982 were made at higher 
altitudes and colder ambient temperatures than 
in 1980 and 1981. Following each flight, the 
airplane was thoroughly inspected for evidence 
of lightning attachment, and the individual 
lightning attachment spots were plotted on iso¬ 
metric projections of the airplane to identify 
swept flash patterns. 

This paper presents samples of the strike 
attachment patterns that were found, and dis¬ 
cusses the implications of the patterns with 
respect to airplane protection design. The 
altitudes, ambient temperatures, turbulence, and 
precipitation at which the strikes occurred are 
also summarized. The significance of these 
findings with respect to lightning avoidance by 
airplanes is also discussed. Finally, lightning 
strikes to several U.S. Air Force airplanes 
struck during routine operations are discussed 
in relation to the Storm Hazards flight 
experience. 

THE NASA LANGLEY RESEARCH CENTER STORM HAZARDS 
PROGRAM is being conducted to improve the state 
of the art of severe storm hazards protection, 
detection, and avoidance, and design of aircraft 
for those hazards which cannot reasonably be 
avoided. Although hazards such as turbulence 
and windshear are being studied, the primary 
emphasis of the Storm Hazards Program is being 
placed on lightning hazard research. Lightning 
is of special interest because the projected use 
of digital avionics systems and composite air¬ 
craft structures will require incorporating 
appropriate lightning protection design 


The flight program began in 1978 when a Twin 
Otter airplane equipped with an airborne light¬ 
ning locator system was flown on the periphery of 
thunderstorms in Oklahoma and Virginia (1).* 

The program continued in 1979 when operations 
began with an F-106B airplane, which was flown on 
the periphery of thunderstorms in Virginia (2). 
From 1980-1982 419 thunderstorm penetrations were 
made with the F-106B in Oklahoma and Virginia 
during which 176 direct lightning strikes were 
experienced. The lightning attachment patterns 
and flight conditions for 1980 are presented in 
(3). 

Since new aircraft are being designed with 
large areas of skin and structure to be made of 
composite materials, improved knowledge of the 
susceptibility of various parts of the aircraft 
surface to lightning strikes is necessary. For 
the purpose of establishing surfaces of different 
susceptibility to lightning strikes on aircraft, 
the Federal Aviation Administration has defined 
lightning strike zones (4). These definitions 
were later refined by Society of Automotive 
Engineers Committee AE4L (5) to accommodate the 
different lightning environments at forward and 
trailing-edge regions of the airplane. These 
definitions, however, do not establish the actual 
locations of the various zones on a particular 
airplane. At present this is accomplished by 
comparing new designs with actual experience 
(when available) of similar-shaped airplanes, or 
by tests in which scale models are subjected to 
simulated lightning strikes. Uncertainties still 
exist, and one purpose of the Storm Hazards Pro¬ 
gram is to clarify some of the more questionable 
aspects of establishing lightning strike zones. 
The purpose of this paper is to summarize the 
results of the Storm Hazards lightning- 
attachment-point analysis and to present the 
flight conditions associated with the direct 
lightning strikes experienced during 1980-1982. 

TEST EQUIPMENT AND PROCEDURES 

The fifteen airborne experiments on the 
F-106B in 1980-1982 are described in (6-10). 

Table 1 summarizes the experiments performed, 
and the locations of the data svstems are shown 
in Fig. 1. 
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TEST EQUIPMENT 

F-I06B Research Airplane - A thoroughly 
Instrumented F-106B "Delta Dart" airplane 
(Fig. 1) is used to make thunderstorm penetra¬ 
tions in the Storm Hazards Program. The two-seat 
cockpit and large J uternal weapons bay for carry¬ 
ing the research instrumentation systems make it 
well suited for thunderstorm research. Details 
on the F-106B airplane and the criteria used in 
choosing the airplane for this mission can be 
found in (2 and 3). 

An extensive "lightning hardening" program 
has been carried out on the airplane and the 
onboard data systems (2 and 11). Prior to each 
thunderstorm season, the lightning hardening 
procedures were verified by ground tests (11) in 
which simulated lightning currents and voltages 
of greater-than-average intensity were conducted 
through the airplane with the airplane manned 
and all systems operating. In addition, all 
systems were monitored during ground tests in 
which simulated lightning waveforms were injected 
into the airplane for calibration of the direct- 
strike lightning sensors (12). Finally, the 
paint has been stripped from the top and bottom 
wing surfaces of the integral wing fuel tanks to 
reduce lightning dwell times in those areas. 

Airborne Data Systems - The F-106B airplane 
is equipped with a number of data systems to 
measure the electromagnetic and environmental 
characteristics of thunderstorms during penetra¬ 
tions. The direct-strike lightning instrumenta¬ 
tion system (DLite) (6 and 13) documents the 
electromagnetic characteristics of direct 
lightning strikes and nearby lightning flashes 
at normal airplane flight altitudes. It con¬ 
sists of twelve electromagnetic sensors mounted 
on the surface of the airplane, a shielded 
recording system in the weapons bay, and a con¬ 
trol panel in the aft cockpit. The DLite record¬ 
ing system consists of two digital, expanded-’ 
memory, wide-band, transient-waveform recorders 
coupled to a 14-track analog tape recorder, and 
one wide-band, 2-channel, analog tape recorder. 

In addition, outputs of a current transformer 
attached to the nose boom and a shunt attached 
to the Kevlar/epoxy vertical fin cap were 
recorded on a Boeing Data Logger system mounted 
in the weapons bay (6 and 7). 

The lateral and vertical components of the 
ambient electric field were measured by a system 
of three field mills (6). The field mill data 
were recorded on the 14-track analog tape 
recorder in the shielded enclosure and were also 
telemetered in real-time to the mission control 
room at NASA Langley to help provide real-time 
guidance to areas of lightning activity. 

The lightning attachments to the airplane 
were filmed by two 16mm color movie cameras (6) 
running at a frame rate of 14 frames per second 
(the shutter is open 1/28 sec for each frame). 

One movie camera was installed forward of the 
windshield, focused on the nose boom. The sec¬ 
ond camera was mounted on the left side of the 
fuselage, looking aft. The field of view of 
this camera included the left wingtip and 


vertical tail as shown in Fig. 2. 

Environmental parameters and flight condi¬ 
tions were measured by the Aircraft Instrumenta¬ 
tion System (AIS) and the Inertial Navigation 
System (INS). Among the parameters measured by 
the AIS were the following: static pressure, 
dynamic pressure, angles of attack and sideslip, 
total air temperature, and normal acceleration. 
Interphone conversations and VHF radio trans¬ 
missions to the mission controllers were recorded 
on a separate tape track. The separate AIS 
14-track analog tape recorder also recorded the 
outputs of the INS, which included latitude, 
longitude, and heading. Time was provided to the 
DLite and AIS tape recorders by a battery- 
operated time-code generator in the AIS. The 
descriptions of the lightning flashes by the 
crew and contacts with the mission controllers 
were recorded on an onboard voice recorder which 
ran continually throughout the flight. 

Ground-Based Radars - For the research 
flights in Oklahoma (1980 and 1981), the NSSL 
Doppler radar at Norman (14) was used to measure 
the precipitation reflectivity data. Addition¬ 
ally, an incoherent 10-cm-wavelength surveil¬ 
lance radar (14) was used to provide air traffic 
control guidance to the airplane. 

For those flights which occurred in 
Virginia, four separate.radars were used. First, 
precipitation reflectivity was measured with the 
NASA Wallops Space Range Radar (SPANDAR). 

Second, an FPS-16 tracking radar at NASA 
Wallops was used to track the C-band transponder 
mounted on the airplane in order to provide the 
SPANDAR crew with real-time information on the 
location of the airplane. On those occasions 
when the INS was not used onboard the airplane, 
the FPS-16 radar data were used to produce plots 
of the airplane ground track. Third, the pre¬ 
cipitation data from the National Weather 
Service WSR-57 radar at Patuxent River, Md., 
were transmitted in real time to a color video 
display in the SPANDAR control room to assist the 
SPANDAR crew in providing real-time airplane- 
penetration guidance. Finally, a 70-cm-wave- 
length UHF radar at NASA Wallops was used to 
detect lightning channels in order to provide 
real-time lightning location information to the 
airplane (15). The seven ground-based experi¬ 
ments associated with the Storm Hazards Program 
in 1980-1982 are summarized in Table 2 and 
described in (15-18). 

TEST PROCEDURES 

Flight Procedures - The thunderstorm pene¬ 
tration procedures, given in detail in (11) are 
outlined briefly in this section. Two guide¬ 
lines adopted from previous thunderstorm pro¬ 
grams were that all flights be limited to day¬ 
light hours t. minimize the threat of flash 
blindness, and that penetrations not be made 
through storm areas having precipitation reflec¬ 
tivity contours over 50 dBZ to minimize the 
chances of encountering hail. In 1980 and 1981, 
the freezing level was chosen as the altitude 
for penetration whenever possible. In 1982, 
however, higher altitudes were chosen, with the 


26-2 








specific altitude determined from the indica¬ 
tions of lightning from the NASA Wallops UHF 
radar. Generally, the pilot set power and indi 
cated airspeed for 300 knots prior to storm 
entry and maintained a constant pitch attitude 
during the penetration, accepting the resulting 
altitude and speed excursions. This procedure 
was the best technique for flight control and 
provided a more accurate measurement of turbu¬ 
lence. The flight observer operated all the 
data systems from the rear cockpit, allowing th 
pilot to give his undivided attention to flying 
the airplane. At anv time during the mission. 


of Penetrations, Direct Strikes, 
am ncaiuy ciaahes - A thunderstorm penetration 
is defined as a data run through a thunderstorm, 
usually at a constant nominal altitude and head¬ 
ing, with the time duration computed from cloud 
entry to cloud exit. In 1982, the penetration 
data have been supplemented with two cloud 
passes, which are data runs through light, indis¬ 
tinct clouds with no discernible precipitation 
and turbulence and indefinite cloud entry and 
cloud exit conditions. The two cloud passes 
were included because three lightning events 
occurred during them. 

The lightning events experienced by the 
F-106B airplane are categorized as direct strikes 
or nearby flashes. A direct lightning strike in 
the Storm Hazards Program is defined as a light¬ 
ning flash which has been confirmed to have hit 
the airplane based on lightning attachment 
points, cockpit voice comments, lightning movies, 
or direct strike waveforms on the onboard light¬ 
ning data systems. A nearby flash is defined as 
a lightning flash which has triggered a sensor 
in the DLite system without attaching to the 
F-106B airplane. 

Following each flight, a determination was 
made of the number of direct lightning strikes 
and nearby flashes which had occurred. The pro¬ 
cedure was as follows: 

1. Conduct an external inspection of the 
airplane to check for lightning attachment 
points. 

2. Determine, from the voice transcripts, 
the times of direct strikes and possible nearby 
flashes called by the flight crew. 

3. Review movie camera films to determine 
the times of direct strikes seen by the two 
onboard cameras. 

4. Review Data Logger data for times of 

lid direct-strike waveforms from the nose boom 
rrent transformer and vertical fin cap shunt 
sed with Kevlar/epoxy fin cap only). 

5. Review DLite data for times of valid 
reforms from the sensors. 

6 . Compile a list of unique direct strike 
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of the airplane surf 
as follows: 

1. Inspect each airplane extremity for evi¬ 
dence of lightning attachment. The evidence was 
usually manifested by spots of molten and 
resolidified metal ranging in diameter from 0.01 
to 1.00 cm (0.004 to 0.4 in.),* usually sur¬ 
rounded by a region of discolored paint. 

2. Closely inspect all surfaces that lie aft 
of the nose boom for additional lightning attach¬ 
ment points that indicate the swept-flash path(s) 
taken following initial attachments to the nose 
beam. This often required the use of a 4 power 
magnifying glass because the diameters of some 

of the swept-flash attachment points were very 
small and the points were hard to distinguish 
from other blemishes on the airplane surface. 

3. Record the location of each attachment 
point and plot the points on isometric drawings 
of the airplane. 

4. Calculate the dwell time t that 
elapsed between attachment points, assuming the 
airplane was travelling at a constant velocity V 
of 182.9 m/sec (600 ft/sec), with the expression 
t = d/V, where d is the longitudinal distance 
between successive attachment points. 

5. Review and correlate findings with the 
pilot's and observer's observations. This step 
was particularly important in separating the 
attachment points produced by each strike after 
a flight in which more than one strike was 
received. 

6 . Postulate, based on the attachment 
points, films from the two onboard movie cameras, 
sensor triggers, and flight crew observations, 
the probable direction from which the strike 
initially approached the airplane, the initial 
and final attachment points, the swept-flash 
path(s), and the point(s) and direction(s) from 
which the flash exited the airplane. 

Reduction 
Parameters - St 


imber of 












The ambient temperature was determined from the 
total temperature measurement (19). The rela¬ 
tive intensities of turbulence and precipitation 
at the times of the lightning strikes and nearby 
flashes were extracted from the cockpit voice 
transcripts. 

DISCUSSION OF RESULTS 

FLIGHT CONDITIONS - The direct strikes and 
nearby flashes for the Storm Hazards '80-'82 
Programs are summarized by year in Table 3, in 
which the number of events having corresponding 
waveforms on the DLite and Data Logger systems, 
camera views, crew calls, and altitude and tem¬ 
perature data are indicated. As shown in 
Table 3, 10 direct strikes and 6 nearby flashes 
occurred in 1980; 10 direct strikes and 22 
nearby flashes occurred in 1981; and 156 direct 
strikes and 26 nearby flashes occurred in 1982. 
The numbers of thunderstorm penetrations are 
summarized in Table 3(c). 

Prior to beginning thunderstorm penetra¬ 
tions with the F-106B airplane in 1980, pub¬ 
lished aircraft strike statistics such as those 
shown in Fig. 3 (from (20)) were used to esti¬ 
mate the best penetration altitudes to use for 
maximizing direct-strike lightning data. The 
data in Fig. 3 show the altitudes at which air¬ 
craft were struck during routine operations, as 
compared with a typical cumulonimbus cloud. 

Based upon those data, most penetrations in 1980 
and 1981 were made at altitudes corresponding to 
ambient temperatures between ±10°C. 

The penetrations, strikes, and nearby 
flashes experienced during 1980-1982 are plotted 
as a function of altitude in 0.6-km (2000-ft) 
altitude bands in Fig. 4, and as a function of 
ambient temperature in 5°C temperature bands in 
Fig. 5. A total of 419 penetrations and two 
cloud passes were made at altitudes ranging from 
2.4 km to 11 km (8000 to 36,000 ft), with a mean 
penetration altitude of 6.5 km (21,400 ft). A 
mean penetration temperature was available for 
322 penetrations (see Table 3( c )), with mean 
penetration temperatures ranging from 10°C to 
-50°C, with an overall mean value of -18°C (see 
Fig. 5). The 419 penetrations and two cloud 
passes resulted in a total of 1208 minutes of 
thunderstorm duration (the 322 penetrations with 
mean ambient temperature values produced 980 
minutes of duration). The distributions of 
penetration time with altitude and temperature 
(Figs. 4 and 5) are very similar to the corre¬ 
sponding penetration distributions. 

Direct strikes were found at altitudes from 
4.3 km to 11.6 km (14,000 ft to 38,000 ft) with 
a mean value of 8.7 km (28,400 ft) (Fig. 4); the 
corresponding strike temperatures ranged from 
5°C to -50°C with a mean value of -32°C. 

Although the peak lightning strike activity was 
expected to occur at or near the freezing level 
(0°C) (20), most of the strikes occurred at 
altitudes above 7.6 km (25,000 ft) and at ambi¬ 
ent temperatures colder than -25°C. Although 
194 minutes of penetration time was spent at 


ambient temperatures between tl0°C, only 7 of 
the 176 strikes occurred at these ambient tem¬ 
peratures. In 1980 and 1981, when the majority 
of the penetration time was spent at ambient 
temperatures between ±10°C, only 5 of the 20 
strikes occurred at these temperatures. In 
fact, the peak strike rates per penetration and 
per minute of duration occurred not at the freez¬ 
ing level, but between the altitudes of 10.4 km 
and 11 km (34,000 ft and 36,000 ft) (Fig. 4), or 
between -40°C and -45°C (Fig. 5). Whereas the 
highest lightning strike rates were related to 
penetrations at higher altitudes and colder 
ambient temperatures, the nearby flashes 
occurred at altitudes from 4.3 km to 11 km 
(14,000 ft to 36,000 ft) (Fig. 4), corresponding 
to ambient temperatures from 0°C to -50°C 
(Fig. 5), with no apparent pattern of occurrence. 

In 1982, there were 156 lightning strikes 
compared to only 10 strikes per year during 1980 
and 1981. This order of magnitude increase in 
strikes was accomplished without an equal 
increase in the number of penetrations that were 
made or in the duration of time spent in the 
storms. Specifically, 239 penetrations and two 
cloud passes were made in 1982 in comparison to 
180 penetrations in 1980-81 (see Table 3(c)); 
and 749 minutes were spent in thunderstorms in 
1982, while 459 minutes were accumulated in 
1980-81. The success of the 1982 program was 
produced by concentrating on higher altitude 
penetrations in 1982 than before, and by using 
the NASA Wallops UHF radar to direct the F-106B 
airplane to areas of high lightning flash rates 
in 1982 (15). 

The lightning strike statistics gathered in 
the Storm HazatJs '80-'82 Programs are in gen¬ 
eral agreement with the results of several other 
research programs, including those of the U.S. 

Air Force Rough Rider Program (21). In this 
earlier flight program, the maximum lightning 
activity in thunderstorms was found to occur 
near 9.2 km (30,000 ft) altitude at about -40°C. 
In addition, ground-based instrumentation systems 
(22 and 23) have detected two layers of lightning 
activity in active thunderstorms. For the 
Florida thunderstorm studied in (22), the tem¬ 
perature of the lower boundary of the higher 
altitude layer of activity occurred at -39°C. 

This higher altitude layer of activity is 
described in (23) as a new class of small inter¬ 
cloud lightning activity which forms a canopy 
over the main updraft region of a storm, produc¬ 
ing a near continuum of lightning activity 
centered at 12 to 13 km (39,370 to 42,650 ft) 
height. Although there are some differences 
between the altitudes of peak lightning activity 
found in (23) and by the F-106B airplane, the 
differences may be due to the difference in the 
height of the cloud tops of the storms sampled, 
with tops in Oklahoma generally exceeding those 
of storms in Virginia. Finally, lightning 
echoes have been detected by the NASA Wallops 
UHF radar at the higher altitudes flown by the 
F-106B airplane in 1982 (15 and 24). Therefore, 
it stems likely that many of the Storm Hazards 
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'82 strikes were taken during flights through 
the upper lightning activity center of 
thunderstorms. 

The lightning strike statistics shown in 
this paper differ significantly from the pub¬ 
lished strike data for commercial aircraft, such 
as in (20) for several reasons. First, each set 
of data was accumulated from different meteoro¬ 
logical and operational conditions. The NASA 
data presented herein were taken from inten¬ 
tional penetrations through active thunder¬ 
storms, while the operational data include 
lightning events which occurred in a variety of 
storm conditions. In addition, the altitudes at 
which commercial aircraft penetrate thunder¬ 
storms were found in (21) to be related to three 
factors: avoidance of obvious thunderstorm 
areas; the duration of lower altitude holding 
patterns in bad terminal weather; and frequency 
of instrument flight rules (IFR) altitude assign¬ 
ments. In other words, commercial aircraft will 
normally deviate from course when enroute to 
avoid thunderstorms which have tops that reach 
to cruise altitudes, and only penetrate when 
required to do so in the terminal area, where 
typical assigned altitudes are near the freezing 
level. Therefore, the lightning event distribu¬ 
tion with temperature found in the Storm Hazards 
Program is the result of the higher percentage 
of time spent by the F-106B research airplane at 
altitudes with high rates of lightning activity 
above the freezing level, compared with the low 
percentage of time spent in storms at tnese alti¬ 
tudes by commercial aircraft in routine 
operations. 

The Storm Hazards data contain several 
implications for aircraft operations in and 
around thunderstorms. First, the use of light¬ 
ning strike statistics, such as those in Fig. 3, 
to avoid lightning strikes can lead to the oppo¬ 
site effect if the data are not fully understood. 
Published strike data have been used to prepare 
pilot briefings that indicate that aircraft 
should fly at altitudes where the temperature is 
colder than -10°C in order to avoid lightning 
strikes. Two examples of such operations are 
given in (2), in which strikes to three U.S. Air 
Force F-106A airplanes on routine missions are 
described. In one case, a pair of airplanes was 
at an altitude of 7.6 km (25,000 ft) and 7.9 km 
(26,000 ft) when the airplanes were struck in 
non-thunderstorm clouds. In the other case, a 
single F-106A airplane received multiple strikes 
while penetrating a line of thunderstorms at an 
altitude of 13.7 km (45,000 ft). In yet another 
(unpublished) case, a U.S. Air Force F-15A air¬ 
plane penetrated a thunderstorm in the same gen¬ 
eral area, altitude, and time as the NASA F-106B 
airplane. The NASA airplane took six inten¬ 
tional direct strikes (strikes 151-156 in 1982) 
without incident, while the F-15A sustained one 
direct strike which caused considerable damage 
to the airplane. None of the pilots of these 
USAF airplanes were anticipating a direct light¬ 
ning strike since they were using accepted 
lightning avoidance techniques. 


A second common misconception among pilots 
has been that lightning, turbulence, and pre¬ 
cipitation occur at the same place and time in 
thunderstorms. The Storm Hazards data contra¬ 
dict this belief. The number of direct strikes 
to the F-106B airplane are plotted as a function 
of the relative turbulence and precipitation 
intensities, as called by the flight crew, in 
Fig. 6. In the majority of cases, the direct 
strikes occurred in areas of the thunderstorms 
in which the crew characterized the turbulence 
and precipitation as negligible to light. How¬ 
ever, a few strikes occurred where the conditions 
were described as being heavy. The nearby flash 
data (not shown) are nearly identical to those 
shown in Fig. 6. In addition, three lightning 
events did not occur during penetrations. Two 
nearby flashes and a direct strike (nearby 
flashes 2 and 16 and strike 56 of 1982) occurred 
in what were called "cloud passes," in which the 
only notable feature was light clouds. In fact, 
three strikes occurred when the aircraft was at 
the edge of the clouds, or actually in the clear. 
Specifically, one direct lightning strike 
(strike 5 in 1980) occurred only 2 seconds after 
cloud entry, another (strike 77 in 1982) occurred 
only when the forward one third of the fuselage 
was in the cloud, and yet another (strike 3 in 
1981) occurred while the airplane was in clear 
air underneath an overhanging cloud. 

For safety reasons, the F-106B airplane was 
not flown into the reflectivity cores of storms 
where the precipitation reflectivity values 
exceeded 50 dBZ. Therefore, no comments can be 
made on the probability of direct lightning 
strikes occurring in such areas. However, the 
data show that, for the thunderstorm areas 
studied in the Storm Hazards Program to date, 
lightning strikes may be encountered at nearly 
all temperatures and altitudes by aircraft in 
the vicinity of thunderstorms, usually where the 
relative turbulence and precipitation intensi¬ 
ties are characterized as being negligible to 
light. Therefore, these data show that the 
presence and location of lightning do not neces¬ 
sarily indicate the presence and location of 
hazardous precipitation and turbulence. 

It is felt that most, if not all, of the 
direct strikes which have occurred to the F-106B 
airplane were intercloud strikes or upward 
branches of cloud to ground flashes, rather than 
the core channel of cloud to ground flashes, 
which conveys the entire flash charge to the 
ground. If this is the case, the peak ampli¬ 
tudes of current and charge should be relatively 
low in comparison to the published full-threat 
values (4 and 5) which are based on ground-based 
measurements of currents in the cloud to ground 
channel, and this has proven to be the case to 
date (6). However, some current rates of rise 
have proven to be very high (6), which implies a 
potential threat to aircraft electrical and 
avionics systems from lightning-induced 
voltages. 

Finally, there is strong evidence that a 
number of the direct strikes to the F-106B 
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airplane were triggered by the airplane (15), 
even when the airplane was in areas with high 
natural lightning flash rates. If this is the 
case, lightning detection devices may not be 
capable of warning an aircraft flight crew of a 
lightning threat. This fact, coupled with the 
poor correlation of lightning occurrence with 
hazardous precipitation and turbulence shown in 
the Storm Hazards Program, implies that light¬ 
ning warning/detection devices may not be suffi¬ 
cient as the sole cockpit device to gi*de air¬ 
craft flight crews around hazardous thunder¬ 
storm areas. 

LIGHTNING ATTACHMENT PATTERNS - In the pre¬ 
ceding section of this paper, the occurrence of 
aircraft lightning strikes with respect to the 
general features of the thunderstorm has been 
discussed. In this section, the discussion 
focuses on the locations of the lightning attach¬ 
ment points on the exterior of the F-106B air¬ 
plane. Four general patterns have been found in 
the lightning attachment patterns: swept 
flashes along the entire length of the fuselage; 
swept flashes which sweep back across the mid¬ 
span areas of the wing; strikes in which the 
initial entry and exit points occur at the same 
extremity; and strikes in which the entry and 
exit points are confined to the aft extremities. 
With any of these patterns, it is possible for 
the swept flash channel to rejoin behind the 
airplane after the airplane has flown through 
the channel as shown by the aft-facing lightning 
movie camera, and described below. 

The first lightning attachment point pat¬ 
tern is the typical nose-to-tail swept flash 
seen on many aircraft. The strike used as an 
example of this pattern (strike 16 of 1982) also 
is an example of a swept flash which rejoins 
behind the airplane after the airplane has flown 
through the channel. At the time of this 
strike, the crew called a strike down the right 
side of the fuselage. The lightning attachment 
points for strike 16 are shown plotted on an 
isometric view of the airplane in Fig. 7. The 
initial entry occurred at the nose boom with the 
initial exit occurring at the left wing tip. As 
the airplane flew through the channel, the entry 
point swept aft along the right side of the 
fuselage to the vertical tail and the exit chan¬ 
nel trailed aft from the wing tip. The general 
orientation of the channel shows the strike 
attaching to the airplane downward from the 
right, and exiting downward to the left off the 
wing tip. Four of the ten successive frames 
taken by the aft-facing movie camera are shown 
in Fig. 8. To aid in the interpretation of 
these data, it should be noted that the airplane 
travels approximately 15 m (50 ft) between the 
center time of each frame, and also travels 
7.5 m (25 ft) during the time the shutter is 
open for each frame. In Fig. 8(a), the entry 
point has already swept back along the entire 
length of the fuselage to the vertical tail, 
with the exit point trailing from the left wing 
tip. The view shown in Fig. 8(a) remains essen¬ 
tially unchanged for the first four frames 


(frame rate of 14 frames/sec). In the fifth 
frame. Fig. 8(b), the airplane apparently has 
flown out of the channel, as two pieces of the 
channel can be seen faintly behind the airplane. 
The sixth frame is blank, but in the seventh 
frame. Fig. 8(c), the lightning channel can be 
seen clearly behind the airplane. There is a 
displacement of the channel which was caused, 
apparently,by the airplane movement while the 
airplane was still a part of the channel. 

Frames 8 and 9 are blank. Finally, the tenth 
frame. Fig. 8(d), shows a burst in the vicinity 
of the afterburner which was counted as a 
restrike, although it may have been an entirely 
new strike. This film sequence shows that it is 
possible for an airplane to fly out of a live 
lightning channel. In this case, the airplane 
was out of the channel in about 0.3 sec 
(4 frames). 

The second pattern, in which swept flashes 
sweep back across the midspan areas of the wing, 
is illustrated in Fig. 9. In this figure light¬ 
ning attachment points for strike 2 of 1980 are 
shown plotted on an isometric drawing of the 
F-106B airplane. This strike is also an example 
of the third pattern, in which the initial entry 
and exit points occur at the same extremity 
(nose boom in this case) without involving any 
other extremity until the flash sweeps aft. The 
lightning strike scenario which follows was 
deduced from the geometry of the lightning 
attachment points and the pilot's comments at 
the time of the strike. (The lightning movie 
cameras were not installed during 1980.) The 
initial entry and exit points are presumed to 
have occurred on the nose boom. As the airplane 
moved forward, the entry point swept down the 
right side of the fuselage and beneath the right 
wing across the midspan area. The exit point, 
on the other hand, swept back to the left wing, 
where the flash branched. One branch swept down 
the leading edge to the wing tip and the other 
branch swept across the top of the wing in the 
midspan area. The general orientation of the 
channel gave an entry from below and to the 
right of the airplane, with an exit off the 
extremities upwards to the left. The exit 
points were presumed to be on the left wing 
because more severe damage was found on the left 
elevon and wing tip than was found on the trail¬ 
ing edge of the right elevon. A more detailed 
description of this strike can be found in (3). 
This is one of several strikes which have 
occurred in which the flash swept back across 
the midspan area of the wing with no upstream 
initial attachment point such as an engine 
nacelle. The flash would normally be expected 
to sweep back along the wing leading edge to the 
wing tip (see (25)). 

The fourth pattern of lightning attachment 
points that has been found consists of those 
strikes in which the entry and exit points are 
confined to the aft extremities of the airplane. 
A photograph of such a strike is shown in 
Fig. 10, taken from the aft-facing movie camera 
on the left side of the airplane. This strike 
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the physical damage found there. (The right 
wing tip is not in the field of view of the 
camera.) This strike is one of several filmed 
in 1982 in which the attachment points were con¬ 
fined to the vertical tail, wing tips, and after¬ 
burner. In fact, evidence of lightning attach¬ 
ment points have been found up to 36 cm (14 in.) 
inside the afterburner nozzle. 

The lightning attachment point data have 
provided further insights into the validity of 
the assumptions used in establishing lightning 
strike zones on aircraft by showing the manner 
in which an airplane interacts with a lightning 
strike channel, especially the manner in which 
flashes sweep aft from initial lightning attach¬ 
ment points. These data confirm that initial 
entry and exit points frequently occur at air¬ 
plane extremities, in this case the nose boom, 
the wing tips, the vertical-fin cap, and the 
afterburner. Swept-flash attachment points were 
observed along the full length of the fuselage, 
as is common in other airplanes of this general 
size, following initial strikes at the nose. 
Unexpectedly, several flashes swept aft across 
the midspan surfaces,of the delta wing. Swept- 
stroke attachments across unswept wings on air¬ 
planes without upstream attachment points such 
as engine nacelles or drop tanks are extremely 
rare (25), and only a few midspan strikes to 
other delta-wing airplanes have been reported. 

The swept-flash attachment paths and burn 
marks found in this program indicate that the 
midspan areas of swept-wing airplanes may be 
more susceptible to lightning attachment than is 
currently believed. If so, careful attention to 
lightning protection may be required for the 
external wing surfaces of integral wing fuel 
tanks in swept-wing airplanes now being designed, 
especially when composite materials and adhesive 
bonding will be used. 

In addition, the lightning attachment 
points found inside the afterburner indicate 
that instrumentation installed in the exhaust 
sections of engines cannot be presumed to be 
immune from lightning damage because of their 
location. Rather, lightning hardening proce¬ 
dures may be required for such installations. 

SUMMARY OF RESULTS 

During the NASA Langley Research Center 
Storm Hazards Program, 419 thunderstorm pene¬ 
trations were made from 1980-1982 with an F-106B 
airplane in order to record direct lightning 
strike data and the associated flight condi¬ 
tions. This study produced the following 
results: 

1. The mean strike altitude was 8.7 km 
(28,400 ft); the mean strike temperature was 


and altitudes in the vicinity of the storms, 
usually where the relative turbulence and pre¬ 
cipitation intensities are characterized as bei 
negligible to light. Therefore, the presence a 
location of lightning do not necessarily indica 
the presence and location of hazardous precipit 
tion and turbulence. 

3. The data confirm that initial entry and 
exit points of strikes frequently occur at air¬ 
plane extremities, in this case the nose boom, 
the wing tips, the vertical-fin cap, and the 
afterburner. Swept-flash attachment points wer 
observed along the full length of the fuselage, 
as is common in other airplanes of this general 
size, following initial strikes to the nose. 
Unexpectedly, several of the flashes swept aft 
across the midspan surface of the delta wing. 
Evidence of lightning attachment points were 
found up to 36 cm (14 in.) inside the 
afterburner. 
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(a) Location of electromagnetic sensors 


Lr uiiiduntru il sensors , 

(b) Location of additional research 
sensors and equipment 

Fig. 1 - NASA Langley Research Center Storm 
Hazards '82 research vehicle 


Plan view of F-106B airplane showing 
field of view coverage 



(b) Sketch of aft camera view 
Fig. 2 - Field of view of aft lightning movie 


Fig. 3 - Aircraft lightning strike incidents 
function of altitude. From (20) 
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Fig. 4 - Thunderstorm penetrations and lightning 
statistics as a function of pressure altitude 
for Storm Hazards '80-'82 
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Fig. 6 - Relationship of lightning strikes to 
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Fig. 7 - Lightning attachment points and strike 
scenario for strike 16 of 1982; flight 82-017; 
June 5, 1982; 18:39:46.8 GMT. (Points shown 
are a composite of the points - from strikes 




















Fig. 9 - Lightning attachment points and strike 
scenario for strike 2 of 1980; flight 80-019; 
June 17, 1980; 22:28:36 GMT. From (3) 


Fig. 11 - Lightning attachment points and strike 
scenario for strike 15 of 1982; flight 82-017; 
June 5, 1982; 18:35:27.8 GMT 



Fig. 10 - Photograph of direct strike to vertical 
tail of F-106B airplane from aft movie camera. 
Strike 15 of 1982; flight 82-017; June 5, 1982; 
18:35:27.8 GMT. First of three frames 
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ABSTRACT 

A program was undertaken to investi¬ 
gate the vortex geometry, particle distri¬ 
bution and electric field distribution in a 
trailing vortex shed from an airfoil tip 
placed in a dust laden flow. A differen¬ 
tial airfoil, positioned in a subsonic wind 
tunnel, generated the trailing vortex. Dust 
conditions, typical of a helicopter landing 
zone, were simulated by seeding the wind 
tunnel flow. 

A vaned probe was used to measure the 
vortex geometry. A particle impact probe 
and a miniature field meter were used to 
map the particle distribution and electric 
field distribution across the trailing vor¬ 
tex. The particles in the trailing vortex 
created a quasi-steady, electric field with 
contours of constant electric field similar 
in appearance to the particle distribution 
contours. 


INTRODUCTION 

Static charge buildup on aircraft has 
always been a nuisance with regard to in¬ 
flight transmission and reception of radio 
messages. In the case of a helicopter, 
static charge buildup is a severe problem 
because it also presents a personnel and 
material handling hazard. An individual 
loading sling loads from a hovering 
helicopter becomes an integral part of the 
discharge path for the static charge on the 


"This work was supported in part by the U.S. 
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surface of the helicopter. The resultant 
discharge causes violent shocks to the 
unloading personnel; and in situations 
where the stores are flammable or explosive 
arcing from the helicopter to these stores 
can result infires and/or explosions. 

Due to the hazard associated with 
electrostatic charging of helicopters, 
studies have been carried out to determine 
means of eliminating or dissipating of this 
build-up charge. Definitions of the pre¬ 
cise charge-generating mechanisms are still 
nebulous. There have been attempts to 
credit the large deposits of static charge 
on helicopters to ions in the engine ex¬ 
haust gases, to contact charging (tribo- 
electrification) occurring when suspended 
particles in the air (water droplets, mist, 
dust, snow flakes, ice crystals) impact the 
surface of the helicopter, and to movement 
of the helicopter through varying electro¬ 
static fields. 

A fundamental aspect of the fluid 
field associated with the flight of a heli¬ 
copter, which does not occur with fixed 
wing aircraft, is that the trailing vor¬ 
tices shed from each rotor tip flows 
downward in a helical path and bathes the 
surface of the helicopter. The trailing 
vortices shed from the tips of fixed wing 
aircraft, flow downstream and eventually 
become dissipated in the atmosphere. 

The purpose of this research was to 
study the charge-generating and separation 
capabilities of a trailing vortex system and 
to determine if the trailing vortex system 
can create a stable electric field and if 
so to determine the parameters which 
control the strength of this field. To ac¬ 
complish this, an experimental program was 
initiated to measure the strength, geome¬ 
try, and charging capabilities of a trail- 













ing vortex produced due to lift on an air- 
fact placed in the test section of a low- 
turbulence, subsonic wind tunnel. To mea¬ 
sure the rotational velocities present in 
the trailing vortex system, a non-rotating 
type vortex meter was designed and fabri¬ 
cated. A homogeneous mixture of particles 
was seeded in to the flow field in such a 
way as to simulate helicopter dusty landing 
zones. A miniature electric field mill was 
used to measure the electric field in 
seeded flow. An impact probe was developed 
which was used to map the particle density 
distribution throughout the flow field. 

Experimental measurement of the trail¬ 
ing vortex geometry and strength was made 
throughout a plane normal to the flow di¬ 
rection at several locations downstream of 
the trailing edge of the airfoil. Measure¬ 
ments made with the electric field mill and 
density probe were made at the same loca¬ 
tions as the vortex measurements. This 
permitted comparion of flow field geometry 
versus electrostatic charging. 

BACKGROUND ON ELECTROSTATIC CHARGING OF 
HELICOPTERS 

The hazards of static buildup on heli¬ 
copters were mentioned briefly in the in¬ 
troduction. To illustrate the seriousness 
of this problem, one has only to refer to a 
publication by M.E. Rodgers (1) where it is 
reported that rescue operations were being 
jeopardized by the frequent electrical 
shocks received by winchmen and survivors. 
In some cases severe shocks had been ex¬ 
perienced, and concern was expressed about 
the possible consequences to survivors in 
poor physical condition. Attempts to pro¬ 
vide conductive paths through the winch- 
men's clothing proved largely ineffective, 
and trailing wires for earthing purposes 
represented greater hazards than the static 
buildup under certain conditions. 

The various mechanisms leading to 
static electrification of helicopters in 
flight have been cited in References (1) 
and (2). These mechanisms are precipita¬ 
tion charging, spray electrification charg¬ 
ing by induction, and charging from engine 
exhaust gases. In Reference (3), the main 
rotor blades of helicopters are said to be 
the chief contributors to the electrostatic 
charge on helicopters. The remainder of 
this paragraph will be devoted to discuss¬ 
ing the details of these charging mecha¬ 
nisms. 

1. Participation charging — occurs 
when charge-carrying particles, such as 
dust, water droplets, snow flakes or ice 
crystals strike some part of the helicop¬ 
ters's surface. The charging currents 
depend upon the space charge density, in¬ 
terception surface area, and velocity of 
impact. 


2. Spray Electrification -- occurs 
due to the shattering of the surface liquid 
film of water droplets on impact with some 
part of the helicopter surface. The shat¬ 
tering of the droplet causes the electrical 
double layer to break and results in sepa¬ 
ration of parts having dissimilar charges. 

3. Charging by induction -- occurs 
due to the helicopter passing through 
regions of varying electric field. The 
absolute charge on the helicopter remains 
relatively unchanged while the potential 
field of the surrounding varies. 

4. Charging from the engine-exhaust 
gases takes place due to the thermally gen¬ 
erated ion pairs in the combustion chamber 
of jet engines. If ions of either polarity 
preferentially recombine, the remaining 
charges of like sign will yield a charging 
current. 

EXPERIMENTAL PROGRAM 

This section of this paper describes 
briefly the experimental equipment employed 
in the measurement of the vorticity con¬ 
tours, and the particle density contours. 
Most of the equipment discussed in this 
section was designed and fabricated speci¬ 
fically for the experiments which were per¬ 
formed during the course of this study. The 
equipment required to produce the seeded 
flow field and the trailing vortex system 
were a subsonic wind tunnel into which a 
differential airfoil was mounted, a parti¬ 
cle feed system, and a filter system. 

WIND TUNNEL AND DIFFERENTIAL AIRFOIL - 
A low turbulence subsonic wind tunnel was 
used to obtain a flow field for the experi¬ 
mental tests. The wind tunnel is 45 feet 
long (13.7 m) from inlet to the end of the 
turning box located at the tunnel exit and 
has the capability of developing test sec¬ 
tion velocities of 187 ft./sec. (57 m/sec). 
The test section is 10 feet long (3.05 m) 
with a 2 x 1.5 ft. (.61 x .46 m) cross 
section and is constructed of wood coated 
with non-pigmented varhish for dielectric 
properties. The trailing vortex system was 
generated in the wind tunnel through the 
use of a differential airfoil. A differ¬ 
ential airfoil is made up of two airfoil 
sections of equal span length and chord. 
It was mounted tip-to-tip across the wind 
tunnel test section. The airfoil sections 
were positioned so that the angle-of-attack 
of one section was equal and opposite to 
that of the other. The airfoils were sup¬ 
ported in the wind tunnel at their quarter- 
chord position by a steel support rod posi¬ 
tioned 2 feet (.61 m) aft of the contrac¬ 
tion section exit. The differential air¬ 
foil was used in this work because it pro¬ 
duces a stable vortex whose position in the 
wind tunnel test section is independent of 
the flow velocity, the airfoil angle of at¬ 
tack, and the distance aft of the airfoils. 
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PARTICLE FEED SYSTEM - The wind tunnel 
flow field was seeded with particles which 
are blown against the wind tunnel inlet 
screen by a nozzle. The particles were 
stored in a spring-loaded hopper capable of 
holding 50 pounds (22.7 kg) of particles. 
The hopper was equipped with a dial indica¬ 
tor to indicate the flow rate of particles 
from the hopper. The particles were drawn 
from the hopper by a venturi located be¬ 
neath the hopper. The spout of the hopper 
was connected to the throat of the venturi. 
Air passing through the venturi mixed with 
the particles and carried them to a mixing 
chamber in the back on the nozzle. 

The nozzle, which directed the parti¬ 
cles against the wind tunnel inlet screen, 
consisted of a cylindrical mixing chamber 
section, a converging-diverging section, 
and a hollow sting running down the center 
of the nozzle. At the end of the sting was 
mounted a deflecting cone. The mixing 
chamber received the mixture of air and 
particles from the venturi and also un¬ 
seeded air from the primary air supply. 
This mixture was exhausted through the con¬ 
verging-diverging section toward the wind 
tunnel inlet. 

A third air-supply line was attached 
to the hollow sting. Air passing down the 
sting was exhausted through jets located in 
the base of the cone. The main purpose of 
the cone was to reduce the core velocity of 
the nozzle exhaust. The air jets at the 
base of the cone helped to shape the nozzle 
exhaust plume. 

FILTER SYSTEM - The wind-tunnel flow 
exited into a large turning box which di¬ 
rected the flow of air and seed particles 
up toward the laboratory ceiling. A system 
to filter the particles from the flow with¬ 
out creating a severe pressure drop, was 
used to prevent the particles from blowing 
throughout the laboratory and to collect 
the particles for reuse. The system se¬ 
lected was a large I8-ft. high, 12-ft. dia¬ 
meter (5.5 by 3.7 m dia), filter bag which 
is made of cotton flannel. 

VORTEX PROBES - The vortex-sensing 
probe designed and fabricated for the 
measurement of the vorticity in the trail¬ 
ing vortex system was capable of determin¬ 
ing the geometry and magnitude of rotation 
of a trailing vortex by traversing a plane 
normal to the axial direction of the vortex. 
The non-rotating type vortex meter used 
senses the rotation of the fluid by means 
of strain gages. The output signal from 
the strain gages (either in millivolts or 
frequency) corresponds to a stress in the 
sensing element due to twist which is 
caused by the rotating fluid acting on the 
probe vanes. The angular motion in the 
fluid due to vorticity was sensed by the 
action of the flow over a set of vanes 
spaced as a cross. The vanes have dimen¬ 


sions of 1/2" x 1/4" (.64 cm x 1.27 cm). A 
strain gage bridge was used to pick up the 
torseon signals receivea. Figure 1 shows a 
typical probe-vane arrangement. 



Figure 1. Vortex Probe Schematic 


The calibration of the vortex probes 
was carried out by rotating the vortex 
probe in a uniform flow attained in the 
wind tunnel facility. The rotation of the 
probe caused torque to be produced on the 
vanes which caused a change in stress in 
the sensor. This change in stress changed 
the frequency of the radio signal being 
transmitted from the calibration system to 
a receiver located outside the wind tunnel. 
This frequency was visually displayed on an 
oscilloscope or monitored with a digital 
frequency counter. The probe was rotated 
oyer a range of angular velocities, and a 
plot of frequency versus angular velocity 
was obtained. Figure 2 presents results 
fro m the calibration. 



RPS 


Figure 2. Vortex Probe Calibration 




32-3 










TEST PARTICULATES - The particulates 
used for seeding the flow were poly-vinyl- 
chloride, (PVC) pellets, and sand. The 
poly-vinyl-chloride pellets were used to 
seed the wind tunnel flow for most of the 
experiments to avoid excessive wear of the 
equipment. 

Poly-vinyl-chloride (PVC-pellets) - 
The PVC pellets were spherical in shape, 
white in color and chemically inert. The 
physical properties were: 

Sp. Gr. (resin) .... 1.4 

Bulk density. 0.47 - 0.58 gm/cc 

Vol. resistivity . . . 10*^ ohm/cm 
Dielectric constant. . 3.7 

Specific Heat .... 0.25 cal/gm - °C 

The results of a sieve analysis for three 
samples of PVC as supplied and the results 
obtained during the course of th’s study 
are given for comparison in Table I. 



#120 PETRO SAND AND #3 GLASS - The 
#120 Petro Sand and 3# glass granules were 
generally spherical although the surface is 
very rough. The materials are 99.9/1 sili¬ 
con dioxide. The particle distribution for 
each material is given in Table:’. 



It was necessary to simulate tne con¬ 
ditions of a typical helicopter landing 
zone. Reference (5) was used to obtain 
typical helicopter dust concentration and 
particle size distribution of particulates 
recirculated through the rotor of a hover¬ 
ing helicopter at a typical helicopter 
landing zone. The PVC pellets and =3 glass 
had particle size distributions similar to 
that in the Phillips Drop Zone. The parti¬ 
cle size distribution for the ?120 sand 
fill between that of the terrain and the 
distribution obtained at a one-foot hover 
height. 

ELECTROSTATIC SENSING-DEVICE - This 
section describes the instrument used and 
results obtained from them in conducting 
measurement of the particles and the elec¬ 
tric field set-up by the swirling particles 
in the trailing vortex. A difficulty, as¬ 
sociated with using these instruments, was 
the lack of knowledge about the order-of- 
magnitude of electric field levels, voltage 
gradients and charging rates to be expected. 

ELECTRIC FIELD MILL - An electric 
field mill was constructed (6) which was 
used to measure the electric field distri¬ 
bution in a plane normal to the wind tunnel 
flow. The electric field meter sensor con¬ 
sisted of a set of fixed conductive vanes 
and a second set of rotating vanes, geome¬ 
trically similar to the first, which were 
grounded. The grounded set of vanes was 
positioned directly above the conducting 
vanes on an insulated shaft which prevented 
charge flow from one set of vanes to the 
other. The rotating vane periodically 
shielded and exposed the insulated probe to 
the electric field in which the field meter 
has been placed. 

The field meter used in this work was 
modified so it could operate in the seeded 
flow field in the wind tunnel. A grounded 
cone was attached to the field mill sensor 
in front of the rotating vane. This cone 
deflected the particles away from the con¬ 
ducting plate of the field mill sensor. 

The sensor plate output of the field 
mill' was connected to a model 610C Kiethley 
electrometer. The electrometer functioned 
as a pre-amplifier and impedance matching 
device. The output of the Kiethley elec¬ 
trometer then passed through a high-low 
band-pass filter to an amplifier. The out¬ 
put of this amplifier was then monitored 
with a digital counter for statistical 
analysis and a lightrbeam oscillograph for 
continuous recording of the field mill 
signal. 

The electric-field mill probe was in¬ 
serted into the x-y positioning device, and 
the output from the sensor recorded as it 
traversed vertically through the trailing 
vortex. Vertical traverses were run at 
several horizontal locations in the vortex. 
The field mill was also used to determine 
















the statistical distribution of voltage 
pulses recorded at discrete locations in 
the vortex. A drawing of the electric 
field mill is shown in Figure 3. 



ELECTROSTATIC DENSITY PROBES - To 
determine the density distribution of par¬ 
ticles throughout the wind tunnel flow 
field, electrostatic density probes were 
fabricated similar to those described in 
References (8, 9, and 10). 

The results of the previous designs of 
electrostatic density probes indicated that, 
with the paticulates used in this work, 
probes which require the passage of parti¬ 
cles through them, are susceptible to clog¬ 
ging. An impact probe used inthis work had 
a small brass cylinder, 0.125 inches in 
diameter and 0.1875 inch long (.05 cm dia 
and .074 cm). These dimensions were chosen 
to minimize disturbances to the flow and 
yet permit local measurements in regions of 
minimum particle density, such as the vor¬ 
tex core. The brass cylinder was soldered 
to 7 the center wire of a co-axial cable. A 
10 / ohm resistor was inserted into the 
probe lead between the probe tip and ground. 
Particles impacting the probe caused a cur¬ 
rent to flow through the resistor which re¬ 
sulted in a voltage pulse across the resis¬ 
tor. Each voltage pulse was amplified by 
an battery powered isolated amplifier with 
a 20 mega-ohm input impedance and an over¬ 
all gain of 1,000. The output of the am¬ 
plifier was connected to a digital counter 
which recorded the individual pulses and 
displayed the running total. The back¬ 
ground noise level of the amplifier was 
less than a microvolt. The amplifier was 
contained within a grounded metal case to 
eliminate pickup of stray signals. 

CALIBRATION OF THE IMPACT PROBE - The 
impact probe was initially tested in the 
exhaust of a small nozzle. Particles were 
introduced upstream of the nozzle and 
through impacts with the duct walls became 
charged. To increase the probability of 


particle impacts prior to exiting from the 
nozzel, pieces of screen were placed in the 
duct to increase the turbulence and in¬ 
crease the probability of particles impact¬ 
ing the duct and/or the screens. The ap¬ 
proximate number of particles per sample 
was determined from the sample distribu¬ 
tions. The weight of each sample used in 
these calibration tests was 0.2 gm. As a 
check on the number of particles per sample 
a ring probe soldered to the impact probe 
was located at the exit of the nozzle. The 
diameter of the circular loop was slightly 
larger than the exist diameter of the noz¬ 
zle. The particles were introduced into 
the calibrations system over a 4.0 minute 
time interval. Eighteen samples of each of 
the three particulates were used for these 
tests. The total number of particles 
sensed by the ring probe was monitored by a 
digital counter. The results of these 
tests indicated the number of particles per 
sample sensed by the ring probe agreed with 
the number calculated using the measured 
distribution for the PVC pellets and was 
more than 50% low for the #3 glass and #120 
sand samples. The error between the number 
of particles sensed by the ring probe and 
that calculated from the distributions of 
the sand samples was probably due to low 
charge levels residing on the sand parti¬ 
cles. 

The circular wire was removed from the 
cylindrical probe tip and the cylindrical 
impact probe was then positioned 0.5 inches 
(1.27 cm) from the nozzle exit in line with 
the nozzle exhaust. The pulses induced in 
the probe due to particles impacting the 
probe and particles passing close to the 
probe surface were recorded. Histograms of 
the results of these tests are also shown 
in Figure 4. The mean value and standard 
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[_ Figure 4. Particle Count Histogram #3 j 

Glass 

deviation of the measurements are indicated 
in each figure. The results of these tests 
were used in determining the effective dia¬ 
meter of the probe. 

X-Y POSITIONER - The x-y positioner is 
an electromechanical device capable of lo¬ 
cating probes accurately within the wind 
tunnel test section. The x-y positioner 
located the probes in a plane aft of the 
airfoil trailing edge. The positioning de¬ 
vice was moveable and could be positioned 
at any point in the wind tunnel test sec¬ 
tion. The probe support hub of the device 
was capable of supporting any and all of 
the instruments used in the course of this 
work. Separate power trains drove the 
probe support vertically or horizonally to 
a positioning accuracy of better than 0.01 
inch (.004 cm). 

EXPERIMENTAL RESULTS 

EXPERIMENTAL PROCEDURE - The experi¬ 
mental measurements obtained cover three 
separate phases; vortex data, electric 
field distribution, and particle density 
distribution. All the results discussed in 
this section regarding electric field mea¬ 
surements and density distributions are for 
PVC seed particles only. 


The vortex probe was used to measure 
the vorticity distribution at several loca¬ 
tions aft of the airfoil trailing edge. The 
probe was moved to the wind tunnel center 
line which coincides approximately with the 
vortex center. Frequency readings were 
taken over ten second intervals. At least 
two and often three frequency readings were 
taken at each point. The probe was then 
repositioned at various locations in a 
plane normal to the wind tunnel flow. 
After recording the frequency readings at 
the desired locations in the vortex, the 
probe was then moved far from the vortex, 
but not into the wind tunnel boundary 
layer, and a background reading obtained. 
To make a vorticity map at another location 
aft of the airfoils, the tunnel was shut 
down. Then the x-y positioner was reposi¬ 
tioned at the desired location aft of the 
airfoils and the vortex probe repositioned 
on the wind tunnel center line. The tunnel 
flow velocity was brought back up to the 
desired velocity and, as previously out¬ 
lined, frequency readings recorded at 
various locations in a plane normal to the 
flow. Vortex data was obtained for two 
flow velocities and two angles of attack at 
various locations downstream from the air¬ 
foils. 

The electric field meter tests were 
performed after the vorticity measurements 
were completed. The electric field measure¬ 
ments were made at five transverse loca¬ 
tions in the plane of the trailing vortex. 
The field meter was positioned at a trans¬ 
verse location, then continuously moved 
vertically through the vortex; the output 
of the field meter was recorded continu¬ 
ously on a light-beam oscillograph. The 
field meter was then positioned at five 
vertical locations over which the continu¬ 
ous trace was taken. At each vertical lo¬ 
cation, the magnitude of the electric field 
was monitored. The output of the field 
meter was connected to a digital counter 
with a variable voltage threshold. By vary¬ 
ing the counter threshold level, the volt¬ 
age pulses below a specified level did not 
register on the counter. The results of 
the statistical measurements were plotted 
on probability graph paper and from this 
plot the mean field strength was obtained. 

The particle density distribution was 
the last phase of the experimental tests. 
The baseline particle distribution was 
measured without airfoils in the wind tun¬ 
nel. The mass flow indicated by the parti¬ 
cle density p.robe was compared against the 
know mass flow from the particle feed sys¬ 
tem. The particle density probe was then 
used to map the density distribution 
throughout the trailing vortex at several 
locations aft of the trailing edge of the 
airfoil. 








VORTEX PROBE RESULTS - The vortex 
probe results are shown in Figures 5 and 6. 
The center of the wind tunnel is at O-hori- 
zontal and O-vertical. The results are 
plotted as lines of constant angular velo¬ 
city, which is equal to half the vorticity. 
Although the wind tunnel flow was stopped 
completely from test to test, the vortex 
location and geometry was not noticeably 
altered. 



Figure 5. Angular Rotation Contours - rps 



Figure 6. Angular Rotation Contours - rps 


The vorticity measurements downstream 
of the 0012 differential airfoil at ?7 
degrees angle of attack were run with the 
wind tunnel flow at 115 fps (35 m/s). Care 
was taken to insure that the airfoils were 
set at equal and opposite angles of attack. 

These contour plots shown in Figures 5 
and S indicate that the geometry of the 
vortex is well defined as close as 12 
inches (.305 m) aft of the airfoils. Study 
of all data taken revealed that the differ¬ 
ential airfoils produce a stable vortex 
whose position in the wind tunnel test 
section is independent of the angle of at¬ 


tack, the flow velocity, and distance down¬ 
stream of the airfoils within the range of 
the test. 

ELECTRIC FIELD METER RESULTS - Elec¬ 
tric field data was taken at Z r 12 inches 
(.305 m) and Z = 24 inches (.61 m). In 
each case, the angle of attack of the air¬ 
foils were set at ±7 degrees and i 12 de¬ 
grees. When the meter was located in the 
proper axial position, it was then posi¬ 
tioned at the desired horizontal and trans¬ 
verse location with the x-y positioner. The 
wind tunnel flow was then brought up to 
speed and the particle feed system turned 
on. The electric field meter then tra¬ 
versed a vertical path through a portion of 
the vortex, beginning 3.5 inches (8.9 cm) 
above the vortex center line and ending 3.5 
inches (8.9 cm) below. At six locations 
spaced vertically across the vortex (4.0 
and 1.5 inches (10.2 and 3.8 cm) above, at 
the center, 1.5 and 2.5 inches (3.8 and 6.4 
cm) below, the electric field meter output 
was monitored statistically. Due to the 
limited capacity of the particle feed 
system, the flow of particles was stopped 
as the meter was repositioned between sta¬ 
tistical readings and during transverse 
movement of the meter. It was necessary to 
shut down the wind tunnel flow completely 
several times during the mapping of the 
electric field in order to collect the par¬ 
ticles and refill the particle-feed hopper. 

The flow field in the absence of the 
airfoils was traversed vertically with the 
electric field meter at five horizontal lo¬ 
cations across the wind tunnel test section. 
Visual observation of the oscillograph 
traces for these traverses indicated a 
slight decrease in the electric field level 
near the south wall of the wind tunnel test 
section. 

The output signal of the electric 
field meter was fed into suitable ampli¬ 
fiers and to the oscillograph or the digi¬ 
tal counter. The oscillograph was used for 
obtaining continuous traces. The digital 
counter was used during the statistical 
monitoring of the electric field. The mean 
value of electric field was obtained from 
probability plots. 

The electric field data and the volt¬ 
age distribution curves foro< = i7°. Z = '2 
and 23 inches is shown in Figures 7 and 3. 
The complete electric-field contours were 
obtained by linear interpolation between 
the locations where the electric field mea¬ 
surements were made. The dotted portions 
of the contours were added on the basis of 
trends indicated by the other contours. 

PARTICLE DENSITY DISTRIBUTION - The 
particle density probe was traversed ver¬ 
tically and horizontally in a plane normal 
to the wind tunnel flow. At each location, 
the number of particles impacting the probe 
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was recorded over two ten-second intervals. 
The average reading was then taken as the 
particles per second passing that point in 
the wind tunnel. 

The result of these measurements are 
shown in Figure 9 which is a plot of the 
particle density distribution without the 
airfoils in the test section. The particle 
density distribution, with no airfoils in 
the test section, indicates a density 



Figure 7. Electric Field Distribution 

(Volts/m) = ±7° Z = 12 inches 

gradient from the top of test section 
toward the bottom. This is to be expected 
due to the gravita¬ 
tional force on each particle. It must be 
emphasized that the gravitational effect on 
the particles occurs during the time the 
particles are leaving the feed nozzle plume 
and being accelerated in the wind tunnel 
inlet from a velocity of 13 ft/sec (4.0 
m/sec) to a velocity of 115 ft/sec (35 
m/sec). After the particles enter the test 
section inlet, the effect of gravity is 
negligible compared to the drag force due 
to the flow velocity. The average particle 
density across the test section is 250 par¬ 
ticles per second. With the differential 
airfoil installed in place, the particle 
density probe was traversed vertically and 
horizontally over a plane aft of the air¬ 
foil. At each location, particle counts 
were taken twice over ten second intervals. 



Figure 8. Electric Field Distribute 
(Volts/m) = ±7° Z = 23 1 


Figure 9. Particle Density Distribut 








The readings were averaged. The density 
distribution at Z = 12 inches aft (.305 m) 
of the airfoils for ±7 deqrees anqle of 
attack are shown in Figures 10 and 11. The 
particle density contours were obtained by 
linear interpretation between measured 
values. 

Due to the limited capacity of the 
particle feed system, the wind tunnel had 
to be shut down several times during the 
measurement of the density distribution 
over the entire plane of the flow field in 
order to collect the particles and refill 
the particle feed hopper. Each time this 
was done, the last monitored point in the 
flow recounted to insure there was no vari¬ 
ation in the particle feed rate or vortex 
geometry. 

EXPERIMENTAL TESTS WITH SAND IN FLOW 
FIELD - Measurements of the electric field, 
density distribution and particle charge 
were made in the wind tunnel flow seeded 
with sand. The procedure for making the 
measurements is the same as discussed pre¬ 
viously for tests run with the wind tunnel 
flow seeded with PVC pellets. Due to the 
abrasive characteristic of the sand, only 
one vertical trace through the core of the 
vortex with the field meter and density 
probe was made for each type of sand. The 
magnitude of the electric field was ana- 



Figure 10. Particle Density Distribution 


lyzed statistically at five locations over 
which the continuous trace had been run. 
Particle charge was measured at several lo¬ 
cations over the same vertical path through 
the vortex core. All measurements made 
with sand in the flow field were taken 12 
inches (.305 m) aft of the airfoil. 



Figure 11. Particle Density Distribution 


A plot of the particle density and 
electric field strength, 12 inches (.305 m) 
downstream of the airfoil with the wind 
tunnel flow seeded with #120 petro sand, is 
shown in Figure 12. The magnitude of the 
electric field was measured statistically 
at five locations on a vertical path 
through the midpoint of the trailing vortex 
The flow velocity is 115 fps (35 m/sec) 
with the differential airfoil set at ±7 
degrees. 

DISCUSSION OF EXPERIMENTAL RESULTS 

VORTEX PROBE RESULTS - The vorticity 
contours are shown in Figures £ and 6 at a 
plane 12 inches aft of the airfoils for a 
velocity of 115 ft/sec (35 m/sec) with the 
airfoils set at ±7 degrees angle-of-attack. 
It should be noted that the magnitude of 
the vorticity contours are specified in 
terms of angular rotation and that this 
value should be doubled to obtain the mag¬ 
nitude of the vorticity. 

ELECTRIC FIELD METER RESULTS - The 
electric field contours for the conditions 
°< = ±7°, 1 = 12 inches and Z = 23 inches 
(.305 m and .58 m) are shown in Figure 7 
and 9 respectively. The wind tunnel flow 
was seeded with PVC pellets for which m = 
1.67 lb/min (.76 kg/min). The shape of the 
electric field contours at Z = 12 appears 
similar to the vorticity contours for a 
vortex which is not completely rolled up. 
The electric field contours at Z = 23 indi¬ 
cate a weakening of the field strength re¬ 
lative to that at Z = 12 in the vicinity of 
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Figure 12. Electric Field and Particle 
Density - Sand 

the center of the vortex. This is to be 
expected since the charge-carrying parti¬ 
cles, which are creating the electric field, 
are moving with time (or axial distance) 
away from the center of the vortex. This 
movement, of the particles out of the vor¬ 
tex center, is due to the centrifugal force 
on the particles caused by the fluids rota¬ 
tional motion. It is interesting to note 
that the electric field intensity at z = 23 
increases in the upper left-hand region. 
This indicates that the particles, in this 
region at Z = 12 are still under the influ¬ 
ence of vorticity bearing fluid which has 
not completely rolled up into the trailing 
vortex. However, this is not in agreement 
with the vorticity results which indicate 
that, for the condition = ±7°, the shed 
vorticity has completely rolled up at a 
distance Z = 12. A possible explanation 
for this occurrence is that the particle 
being heavier than air tends to lag 
slightly behind the fluid due to their 
inertia and, therefore, roll up into the 
vortex at a slightly slower rate than the 
fluid particles. 

The magnitude of the electric field at 
Z = 12 inches was plotted against radial 
distance along rays A, B, B, D, and E, 


Figure 7. The results indicated that where 
the flow has rolled up, the electric field 
distribution is nearly axisymmetric. 

The electric field contours for<X= ±7° 
at Z = 23 is shown in Figures. These con¬ 
tours are similar to what one would expect 
from flow of charged particles in a free 
vortex. When the charged particles reach 
their equilibrium radius, where they are no 
longer under noticeable influence of the 
rotating fluid, the particle distribution 
that should result is one with most of the 
particles in the outer region of the vortex 
and progressively fewer particles nearer 
the vortex center. The individual parti¬ 
cles carry electrostatic charge; hence, the 
electric field would be expected to in¬ 
crease as one moves away from the vortex 
center. 

PARTICLE DENSITY PROBE RESULTS - The 
particle density probe results are shown in 
Figures 3 through 11 as contours of 
particle density (in particles per second). 
The background density recorded with no 
airfoils in the test section indicates a 
homogeneours mixture of PVC pellets in 
horizontal planes within the wind tunnel 
test section but some stratification from 
top to bottom of the flow. This vertical 
density variationis probably due to the low 
speed at which the particles enter the wind 
tunnel prior to the acceleration into the 
wind tunnel test section. A second possi¬ 
bility is that the plume from the feed noz¬ 
zle may have been directed too low on the 
wind tunnel inlet screen. The number of 
particles traveling through the entire 
plane of the 3 ft. sq. (0.8 sq.m.) test 
section was calculated to be 0.8 x 10' 
particles per second. This number compares 
well against the number of particles being 
fed into the test section from the feed 
hopper which is 0.57 x 10' an average 
particle diameter (200//). 

The particle density contours for the 
conditions** = ±7 degrees at Z = 12 and 
Z - 23 are shown in FigureslO and 11 . The 
contours at Z = 12 have an appearance simi¬ 
lar to vorticity contours taken somewhat 
closer to the airfoils. The density con¬ 
tours at Z = 23 are similar in shape to the 
vorticity contours obtained when the vortex 
is completely rolled up, Figures, however, 
enclosing a much larger area. 

The variation of the density along 
radial lines was plotted versus distance 
for several angular positions. The varia¬ 
tion of density was nearly linear and 
indicated the distribution of particulates 
in these test could be assumed axisymme¬ 
tric. 

RESULTS OF TESTS WITH #120 SAND - The 
electric field measurements and particle 
density variations for the #120 petro sand 
is given in Figure 12. The electric field 
profile for the #120 petro sand shows that 
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the electric field varies from 1,800 
volts/m at the vortex center to 2,250 
volts/m outside the vortex core. The par¬ 
ticle density distribution shows that the 
particles are being thrown out of the 
vortex core due to centrifugal force. The 
particle density distribution indicates a 
slight stratification of particle density 
similar to that with PVC pellets in the 
flow. The electric field is somewhat lower 
than that which occurred with the PVC pel¬ 
lets. This is probably due to the differ¬ 
ence in the electrical characteristics of 
the two materials. The major variation of 
electric field intensity occurs near the 
vortex core which is also where the major 
variation in particle density occurs. This 
result agrees with the that obtained from 
tests with PVC pellets. The field gradi¬ 
ents for the PVC tests were greatest near 
the edge of the vortex core. Figure 7, 
where the density gradient is greatest. 
The electric field gradient in the #120 
petro sand is greatest near the vortex core 
also where the largest variation in density 
gradient occurs. This pattern tends to 
verify that the electric field meter is 
actually measuring the electric field set 
up by the flowing particles. 

CONCLUSIONS 

The major conclusions which can be 
drawn from the results of this study are: 

The differential airfoil generated a 
stable vortex whose position in the wind 
tunnel seemed independent of the wind tun¬ 
nel flow velocity, angle of attack and dis¬ 
tance downstream. The core of the vortex 
appeared to be turbulent. The trailing 
vortex system in a flow seeded with par¬ 
ticles established a quasi-static electric 
field distribution due to the charged par¬ 
ticles. The electric field and particle 
density distributions downstream of the 
airfoils are approximately linear with 
radial distance from vortex center and 
roughly axisymmetric at each downstream 
station in the region of measurements. 

The shape of the curves of vorticity 
and particles distribution were quite simi¬ 
lar which indicates a direct dependence of 
the particle distribution on the vorticity 
field. The shape of the electric field was 
quite similar to that of the particle 
distribution. The electric field in turn 
depends on the charged particle distribu¬ 
tion and hence in turn on the vorticity 
distribution. The results of the work do 
confirm the hypothesis that charge separa¬ 
tion and an electric field do indeed result 
as a consequence of the trailing vorticies 
created by lifting surfaces in the presence 
of dust or particulate laden air. 


It should be noted that in corrolary 
work selective charge deposition does 
indeed occur with bodies placed downstream 
in a particulate loaded vortex (11). 
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ABSTRACT 

An investigation of the effects of high 
amplitude simulated lightning currents on the 
tensile strength of T300/5208 graphite/epoxy 
has been conducted. Tests were conducted on both 
unidirectional [0], and quasi-isotropic, [+45, 
-45, 0, 90] s , laminates. Current, ranging from 
3.4 kA to 48.0 kA, was conducted through tensile 
test specimens and the specimens were then 
tested to measure residual strength. Some test 
specimens were conditioned to a moisture level 
of approximately one percent prior to being 
subjected to high amplitude current. In addi¬ 
tion, some specimens were strained to 4000 ycra/ 
cm prior to being subjected to the test current. 
Residual strengths are compared with baseline 
strengths and the results are reported herein. 

LIGHTNING STRIKES TO AIRCRAFT are a major con¬ 
cern because of the potential damage to avion¬ 
ics, fuel systems, and the airframe. Thus, the 
aircraft must be capable of sustaining lightning 
strikes without incurring significant damage. 

The effects of lightning on aircraft can be 
catastrophic. Localized damage due to intense 
heat can occur; however, widespread damage may 
also occur due to current dissipation away from 
the strike attachment point. Airplane damage 
from lightning strikes is described in 
references (1-4)*. Helicopters are also suscep¬ 
tible to lightning strike damage, and some of 
the associated problems are discussed in refer¬ 
ence (5). The certification peak current of 
200 kA and the associated waveform are described 
in reference (6). 

The introduction of advanced composites as 
an aircraft structural material has stimulated 
interest in aircraft lightning protection. The 
high strength, high stiffness, and low density 
of composite materials make them attractive for 
aircraft structural applications. However, the 
electrical properties of composites and metals 
are significantly different. Graphite/epoxy 
(Gr/E) has a relatively high electrical resisti¬ 
vity compared to aluminum. The resistivity of 
Gr/E is 6.0 x 10“ ^ ohm-cm, whereas the resisti¬ 
vity of aluminum is 2.8 x 10”^ ohm-cm. The 
energy absorbed from a lightning current is a 
function of the resistivity of the material 
through which it flows. Therefore, Gr/E will 


(7-11). The visible effects of natural lightn 
on two composite materials are reported in ref 
ence (12). However, information on the rosidu 
strength of composites after simulated lightni 
strikes is limited (ref. 11). 

The purpose of this paper is to report on 
the effects of simulated lightning currents on 
the tensile strength of graphite/epoxv composi 
materials. The test currents were conducted 
through the specimens with no arc gap in the c 
cuit. These tests, therefore, were different 
from tests conducted by arcing the current acr 
an air gap to a test panel and then tensile te 
ing specimens machined from the panel. Tensil 
test specimens were exposed to current levels 
ranging from 3.4 to 48.0 kA to determine the o 
set of strength degradation. Tests were also 
ducted to determine the effects of moisture an 
tensile load on the residual tensile strength 
after a simulated lightning strike. 

Identification of commercial products in 
this report is to adequately describe the mate 
als and does not constitute official endorseme 
expressed or implied, of such products or manu 
facturers by the National Aeronautics and Spac 
Administration. All units of measure used in 
report are given in the Internetional System c 
Units. 

TEST SPECIMENS 

All of the specimens were machined from 8 
ply laminates fabricated from T300/5208 Gr/E p 
impregnated tape. Specimens with unidirection 
[0], and quasi-isotropic, [+45, -45, 0, 90] s , 
orientations were evaluated. The specimen gee 
try is shown in (fig. 1). Each specimen was 3 
cm long to accommodate electrical clamps and 
mechanical grips. A summary of the conditione 
test specimens is given in (table I). 

EXPERIMENTAL PROCEDURES 

Each specimen was mounted in the series R 
circuit represented in (fig. 2(a)). The test | 
cedure consisted of conducting a high amplitud 
simulated lightning current through each speci 
men. A representative test waveform used in t 
experiments is shown in (fig. 2(b)). A speciir 
mounted in the test circuit is shown in (fig. 

In order to determine the effects of mois 





to a strain level of 4000 pcm/ 
:*sentative of the present de- 
jr composite aircraft struc- 
ul 14). The load frame was 
the circuit and the specimen 

of each current waveform was 
oscilloscope trace. The data 
photographs were used to cal- 
jrrent and the associated action 
:ion integral is defined as 
fi 2 dt (1) 

instantaneous current and t is 
id with that current. The action 
>4re of the energy associated 
ll discharge. 

isile tests were completed, sam- 
specimens were used to deter- 
rature rise from resistance 

1 oxidation of the resin in the 
imples were evaluated by diffuse 
is described in reference (15). 
is test was to determine if any 
ion could be explained by 
composite matrix. 

US I ON 

il C.r/Ii specimens were tested 

2 strength data. Although this 
c aircraft structural laminate, 
ton was evaluated as a control 

failure of unidirectional 
r-dominate instead of matrix- 
2 case for cross-ply laminates 
it a obtained for these tests 
(fig. 5). There was no strength 
a current density of 155 kA/cm? 


The residual strengths of the quasi-isotropic 
laminate tests are presented in (fig. 8). The 
effects of moisture are readily apparent. The dry 
specimens were the least affected by current in 
that the onset of strength degradation was at an 
action integral of 5.0 kA^-s. The as-received 
specimens exhibited an onset of strength degrada¬ 
tion at about 2.5 kA^-s and a 31 percent degrada¬ 
tion at 4.25 kA^-s. For the specimens with 1 
percent moisture, the onset of strength degrada¬ 
tion was at about 2.0 kA2-s and a 31 percent de¬ 
gradation at 3.0 kA^-s. 

The data in (fig. 8) indicate that the re¬ 
sidual strength was significantly affected for 
the specimens that were loaded to a strain level 
of 4000 i;cm/cm and simultaneously subjected to 
high amplitude current. The loaded specimens 
exhibited a 6 percent strength degradation at an 
action integral of 0.25 kA^-s and a 31 percent 
degradation at about 1.7 IcA^-s. These results 
indicate that strain increases the susceptibility 
of Or/E composites to damage by simulated light¬ 
ning current. 

The strength data from (fig. 8) are plotted 
as a function of current density in (fig. 9). The 
onset of strength degradation occurs at a current 
density of 55 kA/cm^ which is in agreement with 
data in Reference (7). At this current density, 
55 kA/cm , moisture did not cause a strength re¬ 
duction; however, the specimens strained to 4000 
ucm/cm had a strength reduction of 31 percent. 

The matrix is exposed to higher shear loads 
in the quasi-isotropic laminate than it is in the 
unidirectional laminate. A comparison of the 
residual tensile strength data for the two lami¬ 
nates indicates that a quasi-isotropic laminate, 
(fig. 9), is more susceptible to current damage 
than, is a unidirectional laminate, (fig. 5). 












quasi-isotropic specimens, excluding the loaded 
specimens, was degraded at action integral values 
ranging from 2.0 to 5.5 kA2-s. 

The results of the diffuse reflectance-FTIR 
tests were inconclusive. Some of the specimens 
had an increase in oxygen while others did not. 
The measurements were made over the specimen 
thickness which may have masked possible matrix 
oxidation at the matrix-fiber interface. Such 
oxidation would destroy the bond but may not 
have been detected by the FTIR. More detailed 
analysis and photomicrography are required to 
determine if the failure mechanism is mechani¬ 
cal cracking, oxidation, or a combination of the 

Using the assumption of a 200 kA lightning 
strike to an 8-ply Gr/E laminate, the curve 
faired through the as-received specimen data in 
(fig. 9) was used to develop the plots in 
(fig. 11). The curves for the 16- and 48-ply 
laminates were generated from the 8-ply curve. 

The change in the number of plies changed the 
laminate thickness while all other parameters 
remained constant. This simply changed the 
current density for a given radial distance 
from the assumed lightning strike attachment 
point. In (fig. 11), the tensile strength is 
shown as a function of radial distance from the 
point of the assumed lightning attachment point. 
Although the current flow is dependent on ply 
orientation and lightning attachment points of 
the aircraft, the current was assumed to be con¬ 
ducted uniformly away from the attachment point. 

The results indicate that, for the 8-ply 
(0.12 cm thick) laminate, the residual strength 
is zero out to a radial distance of 2.75 cm. 
Approximately a 10 percent strength degradation 
is indicated at a radial distance of 3.5 cm. 

This radius of damage is less than the 7-cm ra¬ 
dius for a 100 kA strike used in the calcula¬ 
tions presented in reference (10). The action 
integral density reported in reference (10), 
however, is 20 times greater than that indica¬ 
ted in this report to cause strength degradation. 

Data presented in reference (7) indicate 
that for a 100 kA simulated lightning strike, 
matrix cracking may occur 3 cm from the attach¬ 
ment point. This current density level is 
similar to that of the 16-ply curve in (fig. 11). 
The 16-ply curve indicates that strength damage 
would be initiated approximately 4 cm from the 
attachment point. This may indicate that the 
residual tensile strength was degraded at a 
radial distance greater than the visible damage 
reported in reference (7). 

SUMMARY OF RESULTS 

An investigation of the effects of high 
amplitude simulated lightning currents on the 
tensile strength of unidirectional and quasi¬ 
isotropic graphite/epoxy laminates has been con¬ 
ducted. The results of this investigation are 
as follows: 

1. Onset of strength degradation for the 
8-ply unidirectional laminates occurred at an 
action integral value greater than 18 kA^-s. 


2. Onset of strength degradation for the 
8-ply quasi-isotropic laminates occurred at 
action integral values of 5.0, 2.5, and 2.0 
kA^-s, respectively, for the dry, as-received, 
and moisture conditioned specimens. 

3. Onset of strength degradation for the 
quasi-isotropic laminates occurred at signifi¬ 
cantly lower values of action integral than in 
arcing tests reported by previous investigators. 

4. The strength degradation for the quasi¬ 
isotropic laminates loaded to 4000 prm/cm was 
higher than the strength degradation experienced 
by the unloaded laminates. 
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Fig. 9 - Effects of current on the residual 
tensile strength of dry, as-received, 
moisture conditioned, and loaded quasi¬ 
isotropic T300/5208 graphite/epoxy laminates 
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Fig. 10 - Estimated visible damage to un¬ 
protected graphite composites. Figure 
7.39, p. 238, of reference 1 



Fig. 11 - Tensile strength degradatii 
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ABSTRACT 

The critical energy required for the positive streamer-spark transition has been measured 
using an electrode geometry which allows the streamer to propagate in a uniform field. 

The energies are lmJ/cm for breakdown paths in the range 2-11 cm. Much longer delay times 
have been recorded than those observed in short point-plane gaps and these correlate well 
with the dark current. The occurrence of these long delay times allows the possibility that 
the molecular vibrational energy may contribute to the breakdown 






FROM AN ELECTROSTATIC POINT OF VIEW, aircraft 
are commonly regarded as small capacitors 
subject to a charging current, produced by 
p-static or other means, as a result of which 
the aircraft potential is raised to tens or 
hundreds of kilovolts whilst total stored 
energy remains small. The corona which ensues 
limits the ultimate potential achieved and 
this may occur by design at well-located 
passive dischargers cr otherwise. Tanner and 
Nanevicz (l)*and Boulay (2) have demonstrated 
how the electromagnetic coupling between the 
corona occurring at strategically placed 
discharges and avionics can be minimised. 

The aircraft is usually negatively charged so 
these discharges exhibit negative corona 
displaying characteristic high frequency 
Trichel pulses. By contrast, Nanevicz (3) 
in a comprehensive review has identified the 
source of a number of spurious coronas many 
of which may develop into an arc-like 
discharge. As Taillet (4) has indicated, 
these transient discharges exhibit much 
larger currents i and rates of change di/dt 
than simple corona; he suggests 10*A and 
10®As- 1 respectively in times ~100 ns. Such 
events can give rise to very strong coupling. 

The so-called streamering that can be seen 
at a dielectric surface is an example. In 
this case, negative charge accumulates on the 
surface creating large tangential surface 
fields. Positive streamers can then develop 
at adjacent conducting sections of the 
airframe and propagate into the negative 
charge centre. It is the transition of the 
low current, low conductivity streamer 
channel into an arc or surface discharge that 
gives rise to the large and rapidly varying 
currents referred to by Taillet. Some of the 
critical parameters required for this 
transition are considered here. 

BACKGROUND 

Positive streamers often represent the 
first stage of the electrical breakdown of an 
insulating gaseous medium. Unlike their 
negative counterpart, they propagate in 
relatively modest ambient fields (4kV. cm- 1 
at STP) forming tree-like structures as they 
do so. Fig.l. shows a Lichtenberg figure of 
such a structure; this appears to be a faint 
glow to the naked eye which cannot resolve the 
individual branches. Once formed, these 
branches or channels have a radius *'-'30pm and 
a free electron density 10 14 cm- 3 
corresponding to a resistance lMT'cm- 1 . Because 
of this low channel conductivity the linking 
of two electrodes or charge centres by a 
streamer system is not the sole requirement 
for breakdown. In order for this to occur a 
further increase in electron density to 
10 17 cm- 3 is required otherwise, particularly 
in electronegative gases such as air/oxygen, a 
reduction in density and, therefore, conductivity 
takes place because of electron attachment 


(a-n < o). 

A number of criteria have been proffered 
as being crucial in the transformation of a 
streamer channel to a spark i.e. the glow-to- 
arc transition. These include a minimum 
streamer charge in long gaps: Renardiers 
Group (5) and short gaps, Gibert .and Dupuy(-6>, 
a maximum current (i in Fig. 3.). Marode (7) 
and critical streamed energy, Sigmond and 
Goldman (8). This last criterion is related 
to and supported by the hydrodynamic model of 
Marode, Bastien and Bakker (9) who suggest 
that a rise in temperature of the channel, 
produced by the energy released, leads to a 
decrease in neutral gas density at constant 
E-field so that E/N (N = number density) is 
increased to the point where a-n > o and the 
subsequent increase in impact ionisation leads 
to the required electron density for breakdown - 
The critical streamer energies reported by 
Sigmond and Goldman were obtained using small 
point-plane gaps (< 15mm) and single channel 
repetitive breakdown streamers. The 
preliminary measurements reported here involve 
a somewhat different electrode configuration 
which extends the streamer range to 11 cm. and 
attempts to simulate streamer growth in a 
uniform field. 


EXPERIMENTAL 

The electrode arrangement is shown in 
Fig. 2 . Positive streamers develop at the 
needle point located at the centre of a small 
hole in the ground electrode and propagate 
across the gap to the upper cathode. The 
electrodes of 42 cm. diameter provide a uniform 
field in the gap which is variable up to 13 cm- 
the surrounding chamber gives provision for 
gas and pressure variations. A 50fi system 
connected between point x and ground monitors 
the conduction currents flowing in the streame: 
system. For breakdown studies these currents 
may be examined directly with an oscilloscope 
or indirectly using a 20 MHz Transient 
Recorder. The latter is particularly useful 
in the study of air breakdown when pre¬ 
breakdown streamers can occur prior to the 
breakdown streamer itself so that the initial 
streamer or associated luminosity cannot be 
used as a trigger. In such cases, the pre¬ 
trigger mode allows the use of the breakdown 
current as a trigger so that the required 
pre-breakdown currents are recorded. 
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RESULTS AND DISCUSSION 

The time dependence of the breakdown 
streamer currents in gaps up to 11 cm. all 
exhibit the same basic profile shown 
schematically in Fig. 3. The positive primary 
streamer travelling at velocities > 5 x 10 7 cm/s 
crosses the gap in a time which is usually 
small compared to the breakdown time T . The 
low conductivity channel, or more accurately 
because of streamer branching, channels which 
now link the two electrodes give rise to a 
current maximum, i , which occurs after a few 
hundred nanosecondi followed by a decay 
governed by the electron attachment rate 
coefficient v of the gas; a large v will provide 
a large negative ion density in the channel. 

The subsequent non-luminous period preceding 
the breakdown is characterised by a dark current 
i which, in the larger channels, may flow for 
a time t < T fi . The measured values of i imply 
a channel resistance varying from 0.2Mfl/cm in 
short gaps to 3Mft/cm in the larger gaps. The 
breakdown behaviour in the present electrode 
geometry can be compared and contrasted with 
that observed in short point-plane gaps referred 
to earlier. 

CRITICAL ENERGY FOR STREAMER-SPARK TRANSITION 
Because of the uniform field geometry employed, 
the energy per unit length of gap W is 



The values of W obtained in gaps up to 5 cm 
using room air at a pressure of 300 torr are 
shown in Fig. 4. as a function of gap size; 
the minimum breakdown field for all gaps is 
4kV/cm. 

The applied field was increased manually 
so that, in the absence of a starter electron 
near the needle point, the recorded breakdown 
field could exceed the minimum. Apart from the 
smallest gap of 2-3 cm. the unit critical 
energy is approximately constant at ~ImJ/cm. 
This value can be compared with those of 
Sigmond and Goldman (8) using point-plane 
electrodes and humid air where they found 
critical energies of 70 yJ/cm at 1 cm rising 
to 110 yj/cm at 1.4 cm suggesting a parabolic 
relationship between total energy and gap size. 
The present figure of 280 uJ/cm at 2 cm is 
consistent with this behaviour but clearly as 
the gap increases a linear rather than a 
parabolic relationship develops. An explanation 
of this initially rapid rise in unit critical 
energy may be provided by the nature of the 
streamer systems in the two cases. Whereas 
Sigmond and Goldman were confident they were 
using single channel, repetitive streamers the 
present measurements were made, as Fig. 1. makes 
clear, with a heavily branched system. Clearly, 
once branching occurs, the average energy input 
will rise because some will be dissipated in 
channels that do not form part of the path of 
the final arc. Tang (10) has reported an almost 


linear increase in streamer numbers with 
propagation distance in a uniform field which 
would support the concept of a constant unit 
critical energy. Since branching is 
unavoidable in most practical circumstances 
any realistic estimate of the breakdown energy 
required must include the associated losses. 
Although not included in Fig. 4. some recent 
measurements at 11 cm. confirm the working 
value of ImJ/cm. 

PEAK CURRENT i AND BREAKDOWN TIME T . 
Marode (11), again forking with small 
(1-1.5 cm) point-plane gaps, established a 
definitive relationship between i and T g 
which takes the form P 



where i is a minimum current for breakdown 
and A,a°constants. (A = 8 x 10- 3 , a = 1.62 
provided T is measured in ns). This has not 
been confirmed by the present measurements, 
indeed the breakdown time is independent of i , 
which is almost constant and independent of ^ 
gap size with a value of 60 - 70 mA. 

If one ignores the s,mall charge injected 
prior to the attainment of peak current then 
a good approximation to the unit critical 
energy W may be obtained as 


- _5_e 


i 


- 2 . 1 ' 


(-vT b ) 


1 assuming a constant value for W 
then D i v should be independent of T . For 
air, v ^vl0 6 s- 1 so that, provided tne breakdown 
time is greater than~l ps, a constant i 
provides some support for the energy criteria. 
Using the values given above W 0.25mJ/cm., 
in reasonable agreement with the measured value 
bearing in mind the dependence of W on the value 
of v which itself is a sensitive function of E/N 
at the magnitude relevant to this experiment 
(i.e. E/N~50 Td). Of course, if T < 1/v then 
with the same assumption of constant W,the 
reverse argument applies and i will be 
dependent on Tg. Marode's expefimental results 
(equation 2) confirm this behavious, so that is 
is reasonable to conclude that the maximum 
current criteria is a manifestation of that 
based on energy. 

DARK CURRENT i AND BREAKDOWN T . Although 
the unit critial energy is independent of the 
'ion in the break- 


gap size, considerable 
down time is apparent; 
increases with gap size 
times of 1-2 ys, as obsi 
Goldman and others, but 
or more are common. If the gap 
then shorter breakdown times result. 

These times show a marked dependence 
the magnitude of the dark current i (Fig 
Although the dark currents tend to 6c grouped 
according to gap size with large current; 


The 2 cm gaps exhibit 
rved by Sigmond and 
at 5 cm times of 50 ys 
overvolted 


5.) 
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usually occurring in the smaller gaps, the 
occasional short time for a given gap is 
associated with an above average dark current 
(e.g. points A in Fig. 5.) 

The origin of the dark current is not 
entirely clear although it must be electronic 
since ion mobilities are too low to sustain 
currents of a few milliamps. The prevailing 
experimental conditions suggest electron drift 
velocities of —>5 x 10 6 cm/sec so any electrons 
created in the channel during the primary 
streamer transit will either have been swept 
out by the ambient field or attached to an 
electronegative molecule, such as O2 to form 
a negative ion. Thus the electrons which 
constitute the dark current must either be 
generated within the channel itself, possibly 
by a detachment process, or injected from the 
cathode. Marode (7), drawing an analogy 
between a streamer channel and a glow discharge, 
suggests a high field cathode fall may be 
responsible for electron injection from the 
cathode. A small number of secondary electrons 
emitted at the cathode are multiplied by 
collisional ionisation before entering the 
channel proper; he finds a cathode fall of. 
•~400V over 3p. However, in the absence of 
a cathode (e.g. surface streamers on a wind¬ 
shield or the corona-leader structure in long 
gaps) some internal mechanism, such as 
detachment, must be’viable to ensure the 
breakdown. Gallimberti (12) has shown that 
the detachment lifetime of electrons from 
0 r " decreases from 50ys at 700K to lps at 1000K 
wnen the degree of detachment is about 5% of 
of the negative ion population and rises 
rapidly with temperature. Marode et al (9) 
and Bastien and Marode (13)have calculated 
the anticipated channel temperatute in specific 
cases to be several hundred degrees which is 
sufficient to allow some contribution by 
detachment to the dark current. Of course, 
maintaining a dark current of a few mA does 
not explain the onset of the breakdown. The 
channel conditions during the dark period 
suggest an electron density of between 10 13 
-10 35 cm.- 3 whilst 10 17 cm.- 3 is considered 
the minimum for breakdown. The most 
promising explanation for this sudden increase 
is that of Rogoff (11) for hydrogen and Marode, 
Bastien and Bakker (9 ) for air which depends 
upon a decrease of the channel gas density N 
following an expansion of the neutral species 
due to the energy released in the channel. At 
constant field this causes an increase in E/n 
so that a-Ti > o and collisional ionisation leads 
to the final current growth. For this to occur 
a certain minimum energy must be supplied to 
the channel which provides at least qualitative 
justification for the critical energy concept. 
Ignoring losses due to radiation and conduction 
a simple estimate of the rise in channel 
temperature AT is 


where a is channel radius, ; gas density and 
s specific heat, fis the fraction of the 
available energy which contributes to the 
translational and rotational energy of the gas 
molecules and hence the gas temperature. In the 
short gaps for which calculations have been 
made, f< 0.05; the bulk of the remainder being 
stored in the vibrational states of the molecule. 
With f = 1 and a = 30y equation (4) yields 
AT ~1500K using the minimum value for W of 
50pJ/cm. so that when the vibrational energy is 
made available a substantial temperature rise 
should occur. The time constant for vibrational 
relaxation T >> T 0 for short gaps which 
effectively prevents the vibrational energy 
contributing to the breakdown process. At 600K 
this relaxation time is ~ 50ys and falls with 
increasing temperature which suggests that large 
amounts of energy may be dumped in the channel 
provided dark currents continue to flow for such 
a time and thus maintain the channel temperature. 
As seen from Fig. 5. this is possible for 
streamer channels >3cm. 

CONCLUSIONS 

The energy required to give a positive 
streamer-spark transition has been measured and, 
for breakdown paths greater than 3 cm. is 
proportional to the length of the path; at 
300 torr the energy required is ■— lmj/cm. 

Breakdown time T increases with length and 
correlateswell with the magnitude of the dark 
current. For T fi < v- 1 the breakdown time 
depends on the amplitude of the initial streamer- 
current but is independent of this current for 
T > v- 1 . This suggests that the shorter 
breakdown times occur when there is initially 
a sufficiently large energy input to achieve 
the required rise in channel temperature other¬ 
wise the breakdown is delayed until the 
molecular vibrational energy is available. 
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lame no Ming properties as well. 

jor question that has been raised is 
>sible over-pressures from a lightning- 
■tonation wave which might exceed the 
of a honeycomb type arrester. The 
I survey and analysis of shock- 
interaction, as presented in the 
ind illustrated in (Fig. 4), suggest 
over-driven detonation wave will 
ittennate down to a detonation wave. 

:ent study of the effect of shock 
waves , in flame arresters have 
that some arresters do indeed prevent 
; of detonation waves (6). The flame 
lo appear to be effective in providing 
protection if they are properly 
with no gaps around the outside of the 
lich could permit flame passage around 
le. Presumably the detonation wave 
by a lightning shock wave would decay, 
le initial pressures exceeded that of 
it ion wave, into a normal, detonation 
>r some distance of propagation, 
rhis point needs to he experimentally 
in view of the extreme pressures and 
>orted in aircraft fuel vents from this 
l of lightning shock waves and 
(l). 


- tne upper corresponding to detonation waves 
the lower corresponding to deflagration wa 
(flames). The upper range is limited by Rank 
lines labelled 1 and 3 and the lower range 
limited by Rankine lines labelled 4 and 6. 
part EF of Hugoniot curve implies unphysical in 
fluxes and is hence without physical interest. 

Consider the intermediate Rankine 1 
labelled 2. It cuts the Hugoniot curve at 

points, A and B, suggesting two possi 

detonation states. Point B denotes a w 
detonation whereas A denotes a strong detonati 
As the slope of line 2 is reduced, A and 
approach one another to finally coincide at I 
the mass flux approaches that corresponding 
line 3. Point I is known as Chapman-Jouget (C 
detonation state. Similarly, when the mass f 
is low, for any Rankine line falling between 
and 6 there occur two possible (C and 
deflagration states. The lower state is known 
weak deflagration and the upper one, strong, 
the mass flux is gradually increased, th 
finally approach the C-J flame state at point 
From entropy considerations, it could he sh 
that states of weak detonation and deflagrat 
corresponding to sections of the Hugoniot cur 
between I and E and below II respectively, 
unstable and are physically unattainable. 

In the hydrodynamic or Zeldovich-von Neum 





e ways of creating 
distinguished: the 

:h» rase of a flame 
and flu* creation of a 
action of a powerful 
rase the concentration 
ieous detonation may 

rond ca«e, the limits 
>n range. These are, 
sely referred to as 
i and forced ignition. 

has recently been 
>d Ryes on (6). They 
ral data on the 
flowing mixtures of 
> diameter piping at 

• of It meters. Also, 
lit fraction of the 
lit. ion. We note that 
: with the recognition 


'Fig. 4) for a closer 
Lly allowed detonation 
? C-J condition (point 


Hydrodynamic phenomena which have 
aeter of an explosion also arise in aii 
lection with spark discharges. These 
lied hy Mandel'shtam and his coworkers I 


is formed in the air between the electrodes 
immediately after the discharge. As a result of 
Joule heal ing the air in this column is heated tci 
temperatures of the order of many tens of 
thousands of degrees and is strongly ionized (at 
least singly). Propelled by thermal expansion, 
the column expands and acts on the surrounging 
air like a piston, sending into it a cylindrical 
shock wave. (After the shock wave travels a 
distance far exceeding the column length, it 
gradually becomes spherical.) 

In the early stages the air density 
distributions with respect to radius have a 
character appropriate to a cylindrical blast 
wave. Dolgov and Mandel'shtam (II) measured 
these distributions at successive times by 
interferometric methods, according to the 
following electrical circuit parameters in the 
experiment: C = 0.25 pf, L = 2 jih, V = ID kv and 
a discharge gap of 5 mm. The conclusion was 
that, at the beginning when the shock wave is 
still strong (the velocity > about 2 km/sec) the 
density distribution corresponds to that in a 
strong line explosion (14). At a later stage, 
the wave weakens and the counter pressure of the 
air ahead of the shock wave begins to have an 
effect. In this case the density at the 
periphery of the blast wave does not differ 
greatly from standard density, while in the 
central regions the density is very low. Since 
the pressure behind the blast front is, roughly 
speaking, equalized spatially, the temperature in 
the central regions is very high. This central, 
strongly rarefied and high-temperature core is 
the electrically conducting column. The average 
density of air in the column is approximately 10- 
3 of that in the undisturbed air, and the average 
temperature is approximately 40,000 degrees K. 
Spectroscopic measurements give consistent 
results. (Thermodynamic equilibrium in the 








ml.* heat generation, whirl, i, tie term I net! by 
inlf-period of the discharge, is comparable 
(.•ho time over which the shock wave is 
■erf, With this taken into account, the rate 
energy generation was determined 
mentally From shock trajectory data, 


'estimated in the laboratory, hut which take 
ice on a much larger scale, are encountered in 
storm. hight.ning is simply an electrical 

ifhargo am! the thunder is produced leading to 
shock wave which is formed and which 

;eneratee into an nrousi'ic *tnv of loro,, 
fauces. Thivlyuk and Mandel'shtam (I) 
sun,d hy spectroscopic means the average 
peratnve in a lightning column, and found it 
be approximately equal to 20,000 degrees K. 
s value is in agreement with calculations 
od on formulas given hy Rraginskii. (1.6) 
uming for typical, values of the current and 
ie the values 30 kamp and 100-1000 usee (the 
ins of the lightning column i.s 10 cm), 
(mates of the pressure behind the shock front 
’« found to he such that at distances of the 
er of several meter? the thunder can be quite 


It is possible to give a rough quantitative 
timate to the extent of front propagation after 
ich the counter pressure effect begins to 
aken the leading shock wave, In line with the 
cond-order cylindrical bln9t wave decay'theory 
‘+)> ve have in the physical problem a 
■irac ter is t i c length, 

1/2 


where p is the effective lightning blast energy 
pur unit length and P0 i.s the atmospheric air 
pressure. Typically E is on the order of 30 to 
300 Joules/cm for a lightning strike. The 
corresponding length scale R i.s on the order of 
70-200 cms. Until the propagation expansion is 
beyond R], the shockwave remains strong. 

Strong radial expansion of a lightning 
column ‘In a detonahle mixture car, he ‘ 
approximately likened to a pistonlike driving 
pulse. The principal physical phenomenon is 
hydrodynnmic in nature. In the early phases when 
thn detonation wave i.s heavily overdriven, the 
heat released at the detonation wave front is 
unimportant when compared with that deposited in 


behind, If the predetonnt(on mixture has a mass 
density and a Istant chemical energy available 
for release W per unit mass, then, a second 
physical length scale 

1/2 

R - ( F, / w ) , 

2 

say, exists in the problem, At this extent of 
cylindrical expansion, the total energy releaaed 
chemically through detonation becomes comparable 
in importance to the enrgy deposited by the 
lightning pulse. Thus, it may be regarded as the 
minimum distance it would take for a storing 
Mast, wave to gradually settle down to a mode of 
stable detonation propagation. Typical values 
appropriate for a dctonablc hydrocarbon mixture 
at standard atmospheric pressure and temperature 
lead to range of 10-23 cms for R . 

2 

We now reemphasize the relative limits of 
explosion, detonation and flame in mixture 
concentration in a practical situation. From the 
qualitative exhibit of (Fig. 5), we recognize 
tiiat only a fraction of inflammable mixture 
concentrations are inherently capable of shock- 
triggered detonation. Within the latter class, 
only a portion are inherently capable of self¬ 
explosion (flame to detonation transition). 
Conversely, a mixture free of self-explosion 
possibilities may still be vulnerable to the 
detonative initiation of a driving shock wave. 
For aircraft fuel vent lines, there is the 
question of possible over-pressures from a 
lightning-induced detonation wave which might 
exceed those of the simple C-J detonation wave 
for which the conventional flame arresters have 
been tested. In light of the foregoing 
discussions in the present section, several 
aspects of that question can be addressed: 

1. When. atmospheric conditions are such that 
there is spatial gradation of fuel-mixture 
concentration along the tube, leaner toward the 
vent opening, there is the possibility of flame 
ignition by shock wave somewhere inside the tube. 
If die tube length is large enough, greater than 
a few. meters, say, there further exists the 
possibility of having detonation waves generated 
by self-ignition. 

2. If concentration around the vent opening is 
already detonahle, a detonation wave can be 
expected.as a result of a lightning-induced shock 
wave emitted from a strike point within meters 
away from the vent, In this case, a near-stable 
detonation wave may propagate into the vent tube 
through the opening. 

3. If the strike point is very near (within a 
quarter of a meter or so, say,) and if the 
atmospheric conditions of detonahiHty happen to 
prevail out to a few meters from the vent 
opening, then the detonation wave resulting from 
the lightning-induced column expansion can be 
highly overdriven. It is conceivable that the 
process of slowing _down to a stable C-J 
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before wave entry into the vent tube. In this 
case, the detonation wave would continue to 

some additional distance. It is important to 
note that data on the effectiveness of flame 
arresters as supressors of overdriven detonation 
waves are lacking in the literature. 

In conclusion, we observe that the 
probability of development and propagation of 
lightning-initiated detonation waves in the vent 
tube, even overdriven ones, cannot be discounted 
outright. Questions concerning possible schemes 
of prevention or suppression in an existing vent 
tube-flame arrester system warrant more detailed 
investigations. 
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AN OVERVIEW OF THE F-14A AND F/A-18A LIGHTNING TESTS 


R. A. Persia and C. C. Easterbrook 
Electro Magnetic Applications, Incorporated, Denver, Colorado 


ABSTRACT 


During the spring and summer months of 1982, various low 
and high level current injection tests were accomplished on the 
F-14A and F/A-18A aircraft. The intent of this paper is to present 
an overview of the tests and a summary of some of the more signifi¬ 
cant results. The low level tests included swept CW, a fast (—300 
risetime test, and a nuclear EMP related test referred to as the 
Surface Current Injection Technique (SCIT). High level pulses up 
to 173 kA with a 10-90% risetime of 1.7 Msec were injected into a 
powered up active aircraft. Shock excitation tests were also 
performed with voltages on the aircraft on the order of 1 MV. 
Comparisons between lew and high level, fast and slow risetime, 
and CW and pulse responses are summarized. In addition, because 
the two aircraft afford an opportunity to compare metal and carbon 
fiber comoosite (CFCl construction, comparisons between the two 









UHF INTERFEROMETRIC IMAGING OF LIGHTNING 


by P. Richard, A. Delannoy, G. Labaune, P. Laroche 

Office National d'Etudeset de Recherches A6rospatiales (ONERA) 
BP 72 92322 CHATILLON CEDEX 


ABSTRACT 

ONERA has operated in Ivory C oast in 1981 
and in the U nited States in 1982 a lightning im aging 
system. This system is an UHF interferometer which 
gives the location of the emitting sources all along 
the various phases of the lightning discharge. The 
corresponding set-up is described with emphasis on 
the means used for increasing the accuracy. 
Examples of results obtained both with natural and 
triggered lightning strokes, and correlation of these 
results with E field m easure m ents will be presented. 

AN ELECTROMAGNETIC LIGHTNING IMAGING 
SYSTEM using interferometric measurements was 
developed at ONERA (Office National d'Etudes et 
de Recherches Aerospatiales). This equipment was 
first experimented in 1981 during the COPT 
(Tropical Deep Convection) campaign in Ivory 
Coast, in complement to a field meters network for 
monitoring the global electrical activity of storms. 
In 1982, the system was used at the Langmuir 
Laboratory, Socorro, New Mexico for the 
observation of triggered lightning flashes. 

We will present here typical results obtained 
during these two campaigns, with particular 
emphasis on the phenomenology; the analysis has 
been obtained through correlations with E field 
m easure m ents. 

DESCRIPTION OF EXPERIMENTS 

INTERFEROMETRIC SYSTEM - In its princi¬ 
ple, the interferometric location is based on the 
m easure m ents of the phase difference of an incident 
plane wave upon a couple of antennae. The angular 
location of a source is obtained with two 
interferometric couples constituted by a three 
antennae system. The interferometer is built with 
two complementary systems: a large one with a 
10X baseline for an accurate but ambiguous location 
of sources (since phase measurement is obtained 
with a2> ambiguity), and a small system with a X 
baseline for ambiguity removal. At our receiving 
frequency (300 MHz), the baselines are thus 1 and 10 
m long (Fig. 1). 


The system optimization has been performed 
taking into account the instrumental and source 
parameters fll :: . A basic assumption is the "point 
source hypothesis" i.e that within our time 
resolution, the source can be considered as a unique 
point source within to the angular resolution. It is 
anticipated that this condition is fullfllled with a 1 
ps time resolution. 

Among instrumental parameters influencing 
phase measurements accuracy, particular emphasis 
was devoted to the coupling between antennae for 
the small system, to instrumental phase errors 
associated with cables and electronics and to 
baseline dimensional uncertainties. All these errors 
were reduced at the data processing phase either 
through corrections (for coupling), or calibration for 
other instrumental errors. The nominal resulting 
angular accuracy is 5.10 _ 3 in direction cosine. 

The principal source parameters are the signal 
nature and the source size when not a point source. 
Provided the condition (bA>oXD/X)<< 1, it can be 
shown that phase errors due to the loss of 
correlation between antennae, depending on signal 
nature (noise waveform, pulses), is negligible 
compared to other contributions (b is the receiving 
bandwidth, v o the central frequency and DA the 
baseline size in wavelength unit). Extended sources, 
within a /us, would be plotted as a point source at 
the source cente- of m ass. 

Independence with respect to the signal nature 
is one of the most important feature of the 
interfero m etrlc m easure m ent : w hile other location 
methods such as hyperbolic location or DTOA 
(Difference in Time Of Arrival) depends on 
amplitude and shape of the signal for correlations 
leading to delay measurements [2, 3, 5]» the 

iiterfero metric location uses only the phase 
information of signals, thus permitting a high 
dynamic acquisition (80 dB in our case through the 
use of limiter amplifiers) leading in a continuous 
way, unlike other methods, to the source location 
with a 1 us time resolution. 

An example of data processing is presented in 
figure 2 ; we can see the different steps, from the 
raw phase data (in 3ine and cosine), to the image 
reconstruction ; we see in particular the 
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improvement when reconstruction uses the large 
system data. 

Performances limitations are : 

-sensitivity limited to 1 (iV/m.kHz at 300 

MHz 

-problems of ambiguity : still the small 
system has three remaining ambiguities ; in a few 
cases, ambiguity removal was not possible when we 
observed scattered successive sources all over the 
celestial sphere, 

- lack of distance information ; nevertheless in 
some cases, this information is brought back by 
charge location fro m the E field sensor network, 
which is used for the calculation of propagation 
velocities of emission phenomena, (the minimum 
velocity inferred fro m this distance m easure m ent is 
given in each case). 

In 1982, the system has been improved by 
removing the ambiguity ; - a (V2) baseline was used, 
and in 1983, the sensitivity is being increased up to 
10~2 pV/m.kHz ; a real time digital data acquisition 
and pre-processing sub-system is being added ; 
3-dimensional location will be possible using two 
10 km apart interferometric systems, both giving 
angular images of lightning. 

E FIELD SENSORS NETWORK - A ten stations 
network for E field measurements at ground level 
was implemented, during the COPT experiment, 
within a 32 x 32 km area (Fig. 3). Each station 
comprised a field mill sensor with a ±30 kV/m 
dynamic range, and a capacitive antenna with data 
acquisition at 3 different sensitivities : ± 30, ± 300 
and ± 3000 V/m (the damping time constant had 
been set to 10 s). All signals were sampled at a 
800 Hz rate. These measurements have been mainly 
used during two phases of the lightning flash : first, 
for the analysis of the slow charge evolution within 
thunderstorm clouds ; secondly, for the modelization 
of lightning discharges either with unipolar model 
for cloud-to-ground stroke, or dipolar model for 
intracloud activity. [ 6 ] 

The results have been correlated to the 
interferometric locations of electromagnetic 
sources at 300 MHz; we will essentially present 
here the observed phenomena related to the 
different phases of the lightning discharge. 

2. RESULTS ANALYSIS 

For the COPT experiment data processing, our 
attention was focused on two events representative 
of the different types of thunderstorms in this 
region : a storm resulting from diurnal convection, 
and a squall line storm. A total number of 23 
discharges has been studied, most of them of the 
intracloud type. Though we cannot pretend to an 
exhaustive pheno m enologlc description of lightning 
with such a limited number of events, some 
characteristic electromagnetic phenomenon have 
been observed in the different discharge phases 
leading to interesting co m parlsons with E field 
variations and neutralized charges locations. 

GENERAL R E M A R KS - The signal appears to 
be made of a series of pulses, with pulse width 
usually not exceeding the time resolution (1 ps). 


According to the pulse rate, the emission phases can 
be classified on two types : either low rate pulse 
radiation, or high rate pulses bursts. 

Low Rate Pulses Radiation (Fig. 4) is 
characterized by less than 15 pulses/ms ; this 
radiation can last up to several hundred ms ; it is 
mostly associated with spread out sources 
presenting low drift angular velocity, and are 
observed usually during breakdown phases of both 
cloud-to-ground and intracloud discharges. They are 
simultaneous to slow E field variations of large 
amplitude (up to several tens of kV/m) 
corresponding to intracloud neutralizations. 

Pulses Bursts (Fig. 4) - At 300 MHz, these 
phenomena appear as the most intense radiation ; 
they usually last several hundred (is with a mean 
time between pulses of 2 to 4 (is. These bursts are 
associated with fast propagation of sources (10^ to 
7.10? m/s) over large distances (up to 10 km) either 
within the cloud or between cloud and ground. All 
propagations are observed in only one direction at a 
time. They occur in all phases of a lightning 
(breakdown, leader, dart leader, intracloud), and are 
always associated with small E field steps of the 
order of 100 V/m corresponding to small amounts of 
charge transfer or neutralization. When associated 
to large E field slow changes, they frequently occur 
after the largest variations. 

TYPICAL OBSERVATIONS ASSOCIATED TO 
THE DIFFERENT PHASES 0F THE DISCHARGE- 

Intradoud Activity - Both types of radiation 
are observed during this type of discharges. During 
the first few hundreds ms, we observe a low rate 
radiation simultaneous to large E field variations. 
The associated U H F sources locations are spread out 
over distances which can be estimated from a km to 
several km and present apparently slow drift 
velocities (■v 10 1 * - 105 m/s). Usually, bursts are 
observed in the last hundreds ms of the discharge, 
following large E field variations ; the associated 
sources propagate at high speed (>10? m/s) within 
the volume concerned by the previous breakdown 
processes (Fig. 5). 

From their similitudes in electromagnetic and 
E field signatures and chronology, these phenomena 
can be compared to the K process 17 , 10] associated 
to recoil streamers. It is interesting nevertheless to 
point out the difference, of an order of magnitude, 
between es*imated velocities for electromagnetic 
sources ( -10? m/s) and velocities currently 
presented for recoil streamers (1—4.10^ m/s) T 7]. 
Similar phenomena qualified as "Q noise" [8] or 
"solitary pulses" |4| have been observed and also 
present velocities up to 3-107 m /s. 

Breakdown Processes and Leader in Cloud-to- 
Ground Flashes - In precursor to a cloud-to-ground 
flash, processes have been observed to last from a 
few tens ms to several hundreds ms. The shorter 
events can be associated to a leader activity with a 
preponderance of short bursts, this activity being 
preceded for long precursors by processes usually 
within the cloud and typical in some cases of 
intracloud discharges. 

In one of this case of cloud-to-ground flash, 
we observe a 625 m 3 intracloud activity preceding a 
25 ms leader. The first 350 ms presented a low rate 
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pulse radiation (10-15 pulses/ms) associated with a 
spread-out activity expending throughout the cloud ; 
corresponding slow E field variations (total of 
30 kV/m) are associated to large vertical dipolar 
neutralizations between altitudes of 4 km (negative 
charges) and 7 km (positive charges). The last 
275 ms are a superposition of the intracloud 
discharge and breakdown processes prior to the 
cloud-to-ground flash. During the weakening of the 
low rate activity, we observe 4 bursts (Fig. 6) 
corresponding to rapid propagations of sources over 
large distances (probably between 5 and 10 km). All 
propagations are observed to converge towards the 
charge that will be neutralized by the first return 
stroke. Final paths of bursts 1, 3 and 4 seem to 
follow the same downwards trajectories ; respective 
velocities are 1C)6, 5.10? and 1.5.10? m/s. These 
downwards trajectories seem to be limited within 
the cloud. A last burst (4) corresponds to the 
beginning of the leader phase ; the pulse emission 
rate increases during this phase and is associated 
with a large nu m ber of sources scattering fro m the 
breakdown volume. Due to the large extent of the 
activity over the remaining ambiguities of the small 
interferometric system, it was not possible in this 
particular case to differentiate radiation within the 
cloud from possible downwards radiation associated 
with the leader. 

Other cloud-to-ground discharges are not 
superimposed on such an intracloud activity, and it 
was possible to locate electromagnetic sources 
during leader phases. Prior to that, we observed 
either short or long breakdown ; an example of the 
last type is a 300 ms activity previously to a 50 ms 
leader (Fig. 7). The radiation is of the short burst 
type, several tens of ms apart, and is localized along 
a vertical path. In the first hundred of ms, sources 
are shown to be close to the ground, below the cloud 
base, and present in one case, 230 ms before the 
return stroke, a propagation to ground with a 
velocity of 7.107 m/s. The overall topology of 
breakdown associated electromagnetic sources will 
then be t..e place of leader and subsequent dart 
leader sources propagations. Figure 7 presents 
E field variations on four stations and the associated 
electromagnetic sources locations for the 
breakdown and leader phases. Radiation during the 
leader is obr srved only during E steps and 
corresponds fo' the first burst (b) to a downwards 
propagation to the ground (2.107 m/s) and for the 
following bursts either to localized in-cloud 
radiation or upwards propagations (2.107 m/s). E 
field steps can be modelized as long intracloud 
horizontal dipolar discharges driving negative 
charges away from the location of first stroke 
(similar observation is presented in 19l» intersteps 
E field evolutions do not give any detectable 
emission, they correspond to similar but slower 
charge transfers. Charges models indicate that 
during all this leader, horizontal discharges are 
predominant compared to vertical transfer usually 
associated to leader advance. The low time 
resolution of E field measurements cannot give us 
Informations on the chronology, within the ms, of 
charges transfer in the cloud and propagation of 
UHF sources below the cloud. Nevertheless, in this 
case, observations suggest a different image of a 


precursor phenomenon where the breakdown activity 
is already associated to UHF sources defining parts 
of a cloud-to-ground channel. This channel is then 
the place, during leader phase, of propagating 
sources at high speed ( > 107 m/s) well above usual 
leader mechanisms velocities, and associated to the 
transfer of charges within the cloud. 

In a third case of cloud-to-ground flash, no 
noticeable breakdown activity can be observed. 
E field measurements indicate the start of a leader 
10 ms prior to the return stroke ; though the 
distance of the activity (12 km) and the resulting 
imprecision on interferometric location at low 
altitude, electromagnetic sources present a 
propagation towards ground in the first 3 ms of the 
leader activity, subsequent sources for the 
remaining 7 ms are mostly located in altitude 
(Fig. 8). E field measurements show an usual leader 
phenomenon lasting 10 ms [11]. Timing difference 
again suggests that UHF radiation mechanisms are 
not strictly related to a leader progression. In these 
3 examples of cloud-to-ground lightning, we have 
noticed from weather radar data that every case of 
low elevation activity, below the cloud base, is 
associated to heavy precipitations. 

Dart Leader - Dart leader radiation is of the 
burst type ; all data either fro m triggered or natural 
cloud-to-ground lightning flashes show a 100 to 
300 jjs long phase of sources, propagation to ground, 
with velocities between 107 and 5.10’? m/s ; dart 
leader radiation does not come from the cloud but 
from sources descending towards ground (Fig. 9). 
The return stroke is sometimes observed with a 
delay (up to several hundred jjs) ; the sources are 
localized within the cloud, usually away from the 
leader starting point; this suggests that the 
observed delay (or part of it) is associated with the 
upwards and intracloud propagation time of return 
stroke 


Return Stroke - At 300 MHz, the return 
strokes radiation present a large dynamic, going 
from no noticeable radiation to levels comparable to 
burst type radiation. The study of its location shows 
either spread out sources, or more compact groups 
of sources away from locations of precursor 
activity. Some of these phenomena could be 
explained, from an instrumental point of view, if we 
were in presence of strong slm ultaneous sources all 
along the return stroke channel; their location 
would in this case appear to be a point-source, their 
center of mass, moving apparently in an incoherent 
way. Nevertheless, return stroke radiation is not 
continuous and usually not strong, and thus 
instrumental effect cannot explain all the obser¬ 
vations. Such lack of spatial correlation between 
return stroke associated VHF-UHF sources and 
precursors has already been mentioned in the case 
of other location methods [4]. 

While in two cases of first stroke sources 
appear spread out around leader trajectories, in 
most other cases or first or subsequent strokes, they 
are located within the cloud. This suggest that UHF 
emission during return stroke does not come 
essentially from the current wave, but is induced by 
the effect of this wave either around the channel or 
within the cloud (Fig. 9). 
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CONCLUSION 
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MAGNETIC FORCE AND SHOCK DAMACE RESULTING FROM LIGHTNING STRIKES 


ABSTRACT 

The problem of aircraft radome protection 
is of particular importance as this is an area 
where considerable damage has been experienced. 
The probability of a strike to a radome is 


ance or that work are described here in which 
the movement of arcs were studied as well as 
the forces they produced. This has also lead 
to a theoretical appraisal of the impulsive 
nature of magnetic forces generated by lightning 
currents and the effective mechanical forces 
they produce. This theory has been used to 
compare predicted forces with actual forces 
calculated from the deflection of a simple 
mechanical system over a limited range of 
operating conditions. 

The paper starts with a review of the 
various mechanical effects which can result 
from lightning strikes. 

MAGNETIC FORCE AND SHOCK EFFECTS may be 
encountered in lightning interactions with 
aircraft from a number of sources and through 
several mechanisms. The major contributions 
to the problems are: 

1) Magnetic forces (including pinch 
effects and magnetic interaction). 

2) Mechanical shock waves associated 
with arc attachments. 

3) High current arcs in confined and 
semi-confined spaces. 

4) Fusing conductors. 

5) Hydraulic shock in fuel tanks.. 

MAGNETIC FORCES - These forces arise from 

the interaction of magnetic fields generated 
by lightning current pulses flowing in conduc- 


braid 
surface by a tag 
n the di/dt phase 
couraged by indue 


possible creating maximum magneti 
similar situation nearly always a 
a narrow conductor (e.g. a divert 
attached to a larger structure (e 
fuselage). These magnetic forces 


a icxJKA, £ x xu a s puise, involving compressive 
forces of ~320MJ/m * 1 2 3 4 5 (21 tons/in 2 ) . Magnetic 
interaction forces are the subject of the main 
part of this paper. 

MECHANICAL SHOCK WAVES ASSOCIATED WITH THE 
ARC ATTACHMENT - This effect is usually observed 
at the lightning arc attachment ppints on the 
airfract. It is caused by the explosive expansio 
of gas in the arc channel in the longitudinal 
direction producing a shock wave during the early 
part of the first return stroke. Resulting damag 
can vary from denting 2mm thick aluminium sheet 
to splitting 0.5mm thick aluminium or breaking a 
hole through carbon fibre panels less than about 
O.7mm thick. The mechanism is peak current and 
action integral dependent, hence only component A 
of the lightning waveform is significant. 

HIGH CURRENT ARCS IN CONFINED AND SEMI- 
CONFINED SPACES - High impulse current arcs 
produce radial shock waves and very large gas 
overpressure through rapid heating and expansion. 
In confined spaces this can cause considerable 
damage as in cases of punctured radomes etc. 

It can also cause debonding or blow-outs when 
such arcing occurs in voids in insulating layers 
or in glue lines. The severity of the problem 
is dependent on peak current and action integral 
so component A of the lightning waveform is again 
the main culprit. 

FUSING CONDUCTORS - A conductor embedded 
in a non-conducting material will fuse if suffi- 

























cient carbon is present, cause the structure 
to split open. 

A similar effect occurs when a conductor 
is confined only on one side as in the case of 
a radome diverter strip. If a fuseable strip 
is used as a diverter strip, on fusing while 
carrying a lightning current impulse, the 
surface to which the strip is affixed may be 
split. This effect appears to depend on action 
integral, peak current, the quantity of material 
vaporised and on the mechanical properties of 
the supporting material. 

HYDRAULIC SHOCK IN FUEL TANKS - As observed 
above, an explosively fusing conductor produces 
a shock wave in the surrounding gas. Similarly, 
a fusing conductor immersed in a fluid will pro¬ 
duce a shock wave in the fluid and this may in 
certain circumstances have great potential for 
causing damage. Unlike fuel vapour explosion 
or detonation, this effect is more likely to 
occur when fuel tanks are full. 

The remaining part of this paper is 
concerned with the first of these damage mecha¬ 
nisms as described above, and is chiefly directed 
to the problem of magnetically induced forces 
associated with radome protection. The results 
should, however, reflect on any impulse magnetic 
force problems. 

ARC INTERACTION WITH RADOME DIVERTER STRIPS 

The problem of magnetic interaction 
between a lightning arc channel and lightning 
conductor strips has been discussed previously 
by A W Hanson in "Lightning Protection Techniques 
for Radomes having Forward Mounted Pitots",(1). 
There he considers a step leader attaching to 
the pitot boom when there will be some delay 
between the initial leader attachment and the 
time of arrival of the first return stroke. 

During this delay the aircraft will have moved 
forward and hence sweeping of the leader channels 
will occur. The interval between the initial 
leader attachment and arrival of the return 
stroke could vary from 0.5ms to 60ms according 
to the altitude of the aircraft: the altitude 
determines how much further the step leader must 
propagate before the channel is complete for 
the first return stroke to occur. With aircraft 
speeds in the range 66-200m/s the aircraft may 
have travelled between 0.033m and 12m. Hence 
it is possible that a leader be swept over the 
surface of a radome prior to the first return 
stroke in which time re-attachment may be impos¬ 
sible because of the dielectric surface. Parallel 
or near parallel current paths are created by 
the arc and the radome protection strip carrying 

current back to the fuselage. The results of 

this effect (illustrated in Fig. 1) are believed 
to have been seen by the author in at least one 
damaged radome. 

In these circumstances magnetic forces 
will exist between the arc channel and the 
current carrying conductor pushing them apart, 
(Fig. 2). Effective forces as high as 22kN/m 
(15001b/ft) with a conductor spacing of 33mm 
have been produced in laboratory experiments 
using solid conductors and lightning simulation 


current waveforms. As one of the conductors in 
practice may be an arc, a key question is then 
whether the arc moves appreciably in the time- 
scale of interest as this might greatly reduce 
the effective force between an arc and a fixed 
conductor. If an arc does display some "inertia", 
what will be the time resolved displacement of 
the arc and hence what reaction does it provide 
for the forces acting on a protection strip 
attached to a radome wall?. Experiments were 
devised to investigate these questions. 

ARC MOVEMENT 

To investigate this phenomenon it was 
necessary to generate long arcs with high (lOOkA) 
currents. Initiating such arcs is a problem 
because high current capacitor banks do not in 
general produce very high voltages. The simp¬ 
lest method is to use a fine fuse wire of metal 
or carbon across the gap where an arc is required, 
A few kilovolts are then enough to initiate the 
breakdown by fusing the fuse wire in the early, 
low current phase of the discharge. Early tests 
described in (2) on the lOOkV capacitor bank at 
CLSU (shown diagrammatically in Fig. 3) and 
using an Im|con image converting camera framing 
at 2.5 x 10 frames/sec (40ps between images) 
showed problems in recording fuse wire initiated 
arcs. This was due to the intense illumination 
produced by the fusing of the wire and subse¬ 
quently in the vapour from it. Other methods 
of triggering long arcs had then to be sought. 
Sophisticated techniques using high powered 
lasers do exist (3) but the cost and complexity 
are prohibitive. Further tests compared trigger 
electrode initiated arcs with fuse wire using 
the circuit shown in Fig. 4. The Imacon photo¬ 
graphs of this experiment are shown in Fig. 5. 

As can be seen from the photographs, the lower 
arcs show much more detail than the upper fuse 
wire initiated arcs. The important feature to 
observe from the pictures is that the arc does 
not move or break-up significantly until frame 
5 or 6, i.e. 200-240ys after the striking of 
the arc. 

FORCES PRODUCED 

Direct investigations were then made of 
the forces produced by the arc with the arrange¬ 
ment shown in Fig. 6. It was very similar to 
that used in Fig. 4, but a single fuse wire gap 

was used for convenience. The important para¬ 

meter to monitor here was the movement of the 
parallel return conductor for which an optical 
technique was used. 

This displacement monitor was required to 
measure the movement of the return conductor 
fairly accurately and with microsecond time 
resolution. This had to be achieved in the 
presence of large magnetic and electric fields. 

A system using a focused light beam was deve¬ 
loped (4). This comprised a microscope lamp, 
two condenser lenses and a large area photo¬ 
sensitive diode with additional circuitry to 
give a voltage output proportional to the 
amount of light falling on the surface of the 
diode (Fig. 7). 
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Basically the idea was to put an opaque 
flag of minimal inertia on the return conductor 
and arrange the light beam to focus to a point 
just beyond the tip of the flag. The second 
condenser then focused the beam onto the photo¬ 
diode surface. Experiment established that the 
output of the photodiode circuit was linearly 
proportional to the observation of the light 
beam by the flag over the distances of interest. 
The sensitivity of the device could be altered 
by moving rhe focus of the beam towards or away 
from the flag, and the dynamic range by adjusting 
the intensity of the light source. 

Tests were done using the above apparatus 
applying currents with peaks ranging from 60kA 
to 85kA, with various load impedances to produce 
different fundamental frequences and both damped, 
sinusoidal and clamped unidirectional waveforms. 

A typical damped sinusoidal current waveform is 
shown in Fig. 8a with the displacement in the 
lower trace, while a clamped waveform and dis¬ 
placement curve is in Fig. 8b. 

RESULTS OF EXPERIMENTS 

The results are presented in Table 1. They 
have been assembled in ascending values of 
capacitor voltage and added circuit inductance. 
Results from using a solid conductor in place 
of the arc have been tabulated adjacent to the 
arc test results for direct comparison. Direct 
comparison is not straightforward because of 
the difference in total circuit resistance 
between arc and solid conductor tests and 
between clamped amd ringing waveforms. This 
is more fully discussed later. 

The "effective" force on the return cond¬ 
uctor may be expected to depend on the square 
of the current (of the force between two 
current varying conductors) and on the time 
for which the current pulse lasts. The value 
of "effective" force is defined as the force 
which, when applied slowly and continuously, 
would produce the same strain in the mechanical 
system as the transient force. This is calcu¬ 
lated using standard mechanical deflection 
equations. Accordingly, the maximum deflection 
has been plotted againstjthe action integral 
of the current pulse (Ji dt) for each test 
discharge and this is shown in Fig. 9. 

It should be noted that for the ringing 
waveforms, the test circuit included the fault 
limiting resistance (R in Fig. 3) in the capa¬ 
citor bank which dominited the change from the 
static resistance of a solid conductor to the 
time varying resistance of the arc. Thus the 
two lines representing ringing waveform tests 
are very close. In the case of the clamped 
discharges, the fault limiting resistance of 
the bank is removed at peak current so the 
decay is dominated by the load resistance. 

Thus to obtain the same action integral in the 
arc as in the solid conductor, a higher current 
was used for the arc tests. 

From the graph it may be concluded that: 

1) For similar current waveshapes, the 
effects of an arc and solid conductor on the 
movement of a return conductor are of a similar 


order. Bonce the arc has an effect much the 
same as a mechanically stiff conductor over 
these timescales. 

2) Changes in the current waveshape produce 
widely differing effective forces for similar 
action integral. Hence action integral does 
not combine the important parameters in a suit¬ 
able way to produce an overall proportionality. 

3) For a given current waveshape, the 
action integral does display approximate pro¬ 
portionality to return conductor deflection. 

4) To obtain overall proportionality the 
effect of waveshape must be included with the 
action integral parameter. 

THEORETICAL CALCULATION OF CONTINUOUS 
EQUIVALENT FORCE 

The waveform produced by lightning can be 
approximated to a double exponential with a 
fast rise time of the order of a few micro¬ 
seconds and a much larger decay time T from 
lOOys to several milliseconds. 

The type of forces resulting from light¬ 
ning type waveforms have been discussed by 
T E James and J Phillpott in (5). Conductors 
will be subject to high initial impulsive 
forces, followed by a slower and much smaller 
decaying component. The maximum mechanical 
stress will occur at the maximum deflection 
when all the kinetic energy is converted into 
strain energy. 

If we consider a simple mechanical system 
which has an angular frequency ’ w' subject to 
a rapidly applied force P (being the maximum 
magnetic force at peak cu?rent) which decays 
exponentially in a time r/2. This is related 
to what is termed an "effective transient 
force" P which would produce the sane stress 
in the mechanical system when applied slowly 
but continuously. 

A relation can be derived between the 
ratio of the effective transient force (P ) 
and the peak applied force (P ), i.e. P /F , 
and the product of the natural angular frequency 
of the mechanical system 'u'and the exponential 
decay time of the current pulse T, i.e. ut. 

This relationship for a fast rise time 
and exponentially decaying pulse was originally 
derived in 1959 by T E James and subsequently 
by R H Evans in (6) from simple mechanical 
considerations and is reproduced here as Fig.10. 

Using this curve, the effective mechanical 
forces for fast, intermediate, and long types 
of lightning pulses with the characteristics 
shown in Table 2 have been composed. This is 
shown in Table 3. The values of force P Q are 
relative and have been normalised to a value 
of 1.0 for a peak current of 50kA. Value of 
the angular frequency 'u' of the mechanical 
structure have been assumed to lie between 
10^ rads/s for a heavy flexible system and 
lO 4 rads/s for a lighter Btiff one. It can 
be seen from Table 3 that the maximum effective 
force is produced by the fast high current 
component due to the much greater peak force 
P which is produced. 

More general conclusions may be drawn 













from the curve in Fig. 10: it is clear that in 
the curved central area of the graph that wx 
critically affects the effective force and in 
a non-linear function, hence effective force 
will not be a function of action integral 
alone. At the bottom end of the graph, effec¬ 
tive force is approximately a linear function 
of ux so it may be expected that mechanical 
systems with very low u will respond according 
to action integral levels. At the top end of 
the graph the effective force ratio tends to 
a constant factor of 2, so above a certain 
level of lit one may expect effective force 
to scale with peak current only. 

COMPARISON BETWEEN THEORETICAL FORCE 
AND MEASURED FORCE 

In the case of the clamped circuit test 
discharges it is possible to compare the 
effective force calculated from the relation¬ 
ship plotted in Fig. 10 with the effective 
force calculated from the measured displace¬ 
ment of the conductor. This latter calculation 
gives the uniform loading along a simple sup¬ 
ported beam as used for the return conductor, 
which would produce the same maximum deflection 
as measured while passing the high current 
impulse through the circuit. The natural 
angular frequency (w) of the conductor may be 
found directly from the frequency of the dis¬ 
placement oscillations. 

The results of these calculations are 
shown in Table 4. It can be seen that there 
is reasonable agreement between the practical 
results and the theoretical calculations. 

Figure 11 shows displacement calculated 
from the theoretical equivalent force (P ) 
plotted against peak current squared for each 
of the values of T in Table 4 plus some addi¬ 
tional values (ui is kept constant). Also 
plotted are the measured displacements numbered 
to refer to Table 4. Agreement is quite good 
as some inaccuracies would be expected for the 
following reasons: 

1) The waveforms used only approximated 
to exponential decay especially in the case of 
arc conductors. 

2) The mechanical model of the arc/rod 
and return conductor was very simplified assu¬ 
ming simple support and uniform loading 
throughout. 

3) No allowance was made of arc movement 
at long timescales. 

CONCLUSIONS 

During the experiments very large forces 
were recorded resulting from the magnetic 
interaction between arcs and current carrying 
conductors in close proximity. These were up 
to 70% of the forces between two solid conduc¬ 
tors at similar spacing and carrying currents 
of similar magnitudes and waveshape. From 
this it may be concluded that similar design 
considerations must be given to forces arising 
from possible magnetic interaction with arcs 
as is given to magnetic interaction between 
solid conductors. 


Calculations of forces using the theory 
derived for exponentially decaying impulse 
currents have given reasonable agreement with 
those calculated from conductor displacement 
in experiments. It is reasonable to conclude 
therefore that the exponential decay theory 
may be used to indicate the level of forces 
that might be expected in practice at least 
over a limited range of mechanical systems. 

The experimental results confirm that 
in the range of the ux product used (0.2 to 
2 rads) action integral alone is insufficient 
as a measure of mechanical deflection; there¬ 
fore waveshape also must be taken into consi¬ 
deration. However the results indicate that 
for ux < 0.2 the deflection would be a linear 
function of action integral. This is supported 
by the theory as the graph of effective force 
against ux illustrates. In general, to predict 
peak magnetic forces it is necessary to know 
peak current, current waveshape and the anular 
frequency of the mechanical system. 
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Table 2 - Maximum Effective Parameters for Composite Current Waveform 



Table 3 - Effective Transient Magnetic Force P 


41=10 41=10 41=10 41=10 


0.75 1.73 0.75 1.73 


Reproduced from "Simulation of Lightning Strikes to Aircraft” (reference 5) 
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The Lightning and Electrostatic Control Effort 
For the Apache/AH-64 Helicopter 

By 

G. A. Booth and H. M. Hoffart 
Hughes Helicopters Inc. 

Culver City, CA 90230 


ABSTRACT 

The primary flight pattern of the AH-64 Helicopter will subject the vehicle 
primarily to cloud-to-ground lightning discharges as well as cloud-to-cloud. Other 
considerations are the tribo-electric effects. To effect control, a concerted 
review of the vehicles lightning and charging aspects was performed. Deleterious 
effects were examined and corrective measures described in this paper will assure 
greater personnel and vehicle safety while complying with the vehicles mission 
requirements. 

The AH-64 Apache Helicopter uses a four blade main rotor and a four bladed 
tail rotor driven by two gas turbine engines to meet the required flight and per¬ 
formance characteristics. The main and tail rotors, in addition to sections of 
the vehicles body, include composite materials in their design and construction. 
Composite materials are also used for structural support as well as for aerodynamic 
covers for some of the avionics equipment as well as for the gas turbine engines. 
The extensive use of composite materials in the AH-64 vehicle design impacted on 
the control efforts for lightning attraction and discharging plus the problems 
associated with possible lightning penetration into on-board avionics equipment. 
Other areas associated with electrostatic discharging included electromagnetic 
interference/electromagnetic compatibility requirements and effective eletrical 
bonding control for all electrical-electronic units used in the vehicle. To effec¬ 
tively implement the above areas, whatever control measures were exercised, 
the control effort impact was considered for each of the areas outlined above. 

Wick static dischargers are used extensively. Including blade tip static discharge 
wicks, to dissipate vehicle electric potential charges to limit vehicle potential 
levels to only a few thousand volts. Blade wicks also reduce pulse discharge 
levels. The lower peak energy levels per pulse can effect a 20dB to as much as 
lOOdB of quieting for the vehicles communications equipment. Thus, an improvement 
in electromagnetic compatibility is realized, though the initial effort pertained 
to the reduction of vehicle electrification potential levels. Streamer currents 
are also reduced because of the use of multiple wick static dischargers. Elec¬ 
trical bonding designed into the structure and surfaces and effectively implemented 
will reduce arcing that would otherwise be present due to high potential levels and 
high level streamer currents. 


The paper will therefore delineate the controls and measures employed to 
reduce the effects of lightning, static electrical charging and the interaction of 
electromagnetic interference control and electrical bonding. 
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The Sandia Lightning Simulator 
By 

James C. Bushnell and J. G. Kostas 
Sandia National Laboratories 
Albuquerque, New Mexico 87185 


ABSTRACT 

Sandia has developed a lightning simulator to produce currents equivalent to 
those of severe natural lightning, (99 percentile), on external loads on the order 
of 8 microhenries. The simulator can produce from 1 to 4 pulses, with the first 
pulse having peak current of 200 klloampere, rising to this value in 2 microseconds 
with an e-folding decay time on the order of 75 microseconds. The interpulse 
spacing may be varied over a wide range. A continuing current on the order of 300 
amperes for 1 second may also be included. The total action Integral is as high 
3 x 10^ ampere 2 - second, and total charge transfer 300 coulombs. The design and 
implementation of this facility and associated Instrumentation will be discussed. 
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Rectangular Pulse Criteria 

o Each generic interface circuit must 
withstand four unipolar rectangular 
pulses in each polarity. 

o Each test pulse shall have an open 
circuit voltage of at least 500 volts 
and a source impedance of less than 
100 ohms. 

o Pulse width shall be at least 5 micro¬ 
seconds in duration with a rise and 
fall time not greater than 100 nano¬ 
seconds . 

o Each test of rectangular pulses shall 
be applied between each line and 
equipment ground (common-mode) and also 
between the wires of any parallel pairs 
(transverse-mode). 

o Four samples of each generic test 
circuit, where feasible, shall pass 
this criteria or hardening techniques 
shall be evaluated. 

Damped Sinusoid Criteria 

o Each generic interface circuit must 
withstand 15 damped sinusoid voltage 
pulses in a time period of one second 
or less. This step was used primarily 
as a comparison to the 500-volt rec¬ 
tangular pulse criteria; the energy con¬ 
tent of one rectangular pulse is roughly 
equivalent to that of 15 damped sinusoid 
pulses. 

o Each test pulse shall be a damped sinu¬ 
soid waveform with a frequency of 1 
Megahertz (+50%), an initial peak voltage 
of 500 volts (+10%), and a second peak 
voltage of greater than or equal to -300 
volts with a generator source impedance 
of 100 ohms at 1 Megahertz. 

o Initial rise-time shall be less than or 
equal to 100 nanoseconds. 

o Each series of 15 voltage transient 

spikes shall be applied between each line 
and equipment ground (common-mode) and 
also between the wires of any parallel 
pairs (transverse-mode). 

o Four samples of each generic test 

circuit, where feasible, shall pass this 
criteria with power applied or hardening 
techniques shall be evaluated. 


Figure 1 illustrates the difference in the 
applied pulse criteria waveforms. The similarity 
in high frequency spectrum content and amplitude 
of each pulse is shown in figure 2. The major 
difference between waveforms is the dc component 
of the rectangular pulse which usually produced 
the thermal breakdown in the devices tested. 




Fig. 2 - Relative spectral content of applied 
pulse criteria 
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TEST SETUP AND PROCEDURE 

Each generic circuit was mounted in a 
printed circuit board and plugged into an inter¬ 
face box which provided access to circuit power, 
the pulse generator, and monitoring points. This 
interface box was then connected to the appro¬ 
priate pulse generator test load as shown in 
figure 3. After each test, this arrangement was 
used to evaluate the pass/fail status of the cir¬ 
cuit. 

The schematic of figure 3 identifies the 
voltages monitored with a 100-MHz, differential 
oscilloscope. Vq was used to establish the 
applied pulse amplitude while Vl identified the 
voltage developed across the generic circuit 
interface. An 80-ohm current limiting resistor, 
which provided the source impedance, was used to 
stabilize the generator output and maintain the 
proper waveform throughout the test load dynamic 
range. This test arrangement made the pulse 
generator relatively insensitive to test circuit 
impedances, allowing the generator to see 31 
ohms for a short circuit test load and 50 ohms 
for an open circuit test load. 

Since V G was a known applied quantity (i.e., 
the applied pulse), all pertinent circuit quanti¬ 
ties could be determined from the transient re¬ 


known transient voltage of V G and the measured 
transient voltage of or V R were numerically 
combined and integrated by trapazoidal methods 
to evaluate the energy dissipated over the pulse 
duration. These voltage and energy manipula¬ 
tions were easily applied for the rectangular 
pulse; however, phase changes and waveform dis¬ 
tortions created by the test circuit made the 
damped sinusoid pulse calculations too tedious 
except for general approximations. 

In general, the basic test procedure in¬ 
volved applying the four-pulses/500-volt rec¬ 
tangular criteria in each polarity and, when 
device quantities permitted, applying the 15- 
pulse sinusoid criteria. This step was accom¬ 
plished for each interface junction at 100-volt 
increments up to 500 volts. Pretest and post¬ 
test functional checks were performed after each 
pulse to determine the pass/fail status of the 
device. If a failure occurred, appropriate 
protection was devised based on allowable design 
constraints and testing was repeated. 

After successfully demonstrating the required 
damage survivability, the circuit was tested for 
functional upset using the sinusoid criteria. 
Upset responses could not be eliminated on most 


sponse of or V R . The decision as to which of 
the two voltages would be measured was based on 
which voltage was anticipated to have the 
smallest peak amplitude since this would provide 
the best resolution of the test voltages. The 


circuits due to design constraints; consequently, 
these upsets were reported for evaluation to the 
flight control designers and provisions were made 
for system level tests in a lightning pulse 
chamber. 



1 


Fig. 3 - Basic test configurate 















































TEST SUMMARY AND RECOMMENDATIONS 

Test results on a few of the interface 
circuits are summarized in Table 1. In each 
case, device damage protection was achieved with 
little impact on circuit design. Most of the 
circuits used in the DFCS were protected with 
series resistance and/or shunt capacitance 
usually already included in the design. 

Transient suppression devices were only re¬ 
quired on circuits where current or frequency 
requirements prohibited resistors and/or 
capacitors. In a few cases, devices were 
able to withstand the criteria pulses in one 
or both polarities. In these situations, 
internal damping or suppression accounted 
for the survivability. 

All devices produced some response to the 
damped sinusoid criteria. Selection of different 
protection components modified the duration and 
severity of response. In general, the output 
response of a device was a variation on the 
typical output signal and lasted the duration of 
the lightning transient. The exception to this 
statement were those devices which could be 
driven into saturation, such as operational 
amplifiers. In these cases the saturation re¬ 
covery time and circuit loading determined the 
duration of device response. 

During the course of the program many ob¬ 
servations were made concerning generic circuit 
testing. These observations and some associated 
recommendations are listed below: 

1) Sudden polarity reversals from the 
sinusoid produced no observable degrada¬ 
tion as originally expected even though 
devices failed both positive and 
negative rectangular pulses. The pre¬ 
dominate rectangular failures were 
usually identified by device breakdown 
late in the pulse, indicating thermal 
failure. Sinusoid failures were rare 
but usually occurred on the leading edge 
of the transient and was duplicated in 
the rectangular pulse. 

2) Comparison of the power on/off states 
revealed that solid state devices 
failed at lower generator voltages when 
power was applied. Typical circuits 
showed a reduction in load impedances 
with power applied which resulted in 
higher current levels into the device. 

3) A decoupling capacitor was required in 
most cases to isolate the device from 
the low pulse-generator impedances. 

With device i.iputs decoupled, the cir¬ 
cuit under test was pulsed in each state 
of operation. 


4) "Equivalent" replacement parts do not 
have equivalent responses to lightning 
transients. Different manufacturers oft 
use variations in device technologies, 
chip layout, and manufacturing technique 

Variations in functional response to the 
lightning transients were the most commo 
difference between replacement parts. 
Telltale clues for differences in damage 
characteristics were individual gate and 
total package power handling capability, 
maximum voltage ratings, and average 1if 
ratings. 

5) Devices with built-in voltage suppressio 
in one or preferably both polarities 
should be selected. Internal voltage 
suppression reduced or eliminated dis¬ 
crete component requirements for lightn¬ 
ing protection and could reduce the 
device saturation time experienced with 
operational amplifiers. 

6) Devices which have been tested for 
static discharge requirements (2) may 
provide a good source of lightning 
hardened devices. Static discharge 
tests usually involve much higher 
voltages for shorter durations than the 
lightning test but many of the failure 
mechanisms would be common to both. 

7) Devices which have higher trigger 
thresholds for state changes and/or use 
higher operating voltages should be 
selected. These devices were more 
capable of dissipating the energy from 
coupled lightning transients. 

8) Avoid the use of multi-port devices in 
which the exposed interfaces share a 
common node within the device as with 
the LSI 2793310-1. In this case a 
failure of one device input resulted in 
the loss of four to eight input dis¬ 
cretes . 

9) High speed devices with low failure 
voltages and/or wide bandwidth require¬ 
ments (TTL, ECL, etc...) are difficult 
to protect from lightning transients. 
Suppression devices for these circuits 
must operate very close to maximum volt¬ 
age tolerances; any significant overshoo 
from suppression delays can quickly de¬ 
grade or damage the device. Wide band¬ 
width requirements minimize the use of 
filters and series resistances at the re 
quired values because they would prevent 
proper device operation. 
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The best alternative is to avoid this 
type of device in the I/O interface; how¬ 
ever, if nothing else will satisfy system 
requirements then the I/O lines should 
make use of good twisted-shielded pairs 
or coaxial cables. Additional testing 
is required in this area of device tech¬ 
nology to investigate available options 
and their degree of protection. 

CONCLUSIONS 

The generic circuit testing of the DFCS was 
successfully completed with all interface cir¬ 
cuits protected. No outstanding protection tech¬ 
niques were required and many devices were in¬ 
herently "hard" by the nature of the circuit de¬ 
sign. In fact, most circuits are expected to 
show damage survivability to two or three times 
the levels applied simply due to the interface 
components selected as part of normal design 
considerations. 

Interface circuit responses to a lightning 
transient are expected to cause some problems; 
however, these problems only extend to the fact 
that a response occurs. Normal system design 
usually accounts for unexpected transients at 
interfaces by using various hardware and software 
filtering techniques, redundancy, and persis- 
tance. The final concern then is to assure that 
secondary coupling beyond the interface circuit 
does not occur above normal logic thresholds. 
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LIGHINING TESTS OF AN ELECTRONIC ENGINE CONTROL 
WITH ASSOCIATED WIRING 

J. Itobb, D. Loncrini, H. Qgasian and P. Geren 


ABSTRACT 

Impulse and swept CW Lightning tests have been carried out 
as part of a developmental program on a new Electronic Engine 
Control (EEC) for the next generation of transport aircraft. The 
tests were made using an experimental engine as a test bed, and 
were carried out with the controller and with a variety of wiring 
harnesses. 

Digital electonic engine controls are generally used on 
modem high bypass engines which have relatively poor shielding. 
The controller is mounted on the fan case and because of the 
requirements for easy maintnance and quick access, only a few 
attachment points, such as hinges and latches, can be used for 
the fan cowl. This lack of bonding points greatly reduces its 
shielding effectiveness. 

The tests on the Electonic Engine Control (EEC) included 
(a) both impulse and swept CW tests with open and short circuit 

measurements of the induced voltages and currents at the 

controller to determine the basic coupling parameters, (b) low 
level tests with dtrany loads at the controller to provide 
estimates of the voltages vrtiich would actually be seen at the 
pins of the EEC as a function of lightning current anplitude and 
rate of raise and (c) upset and damage tests of an operating EEC. 
Bulk cable currents and magnetic fields about the engine were 
also measured for correlation with the cable conductor 
measurements. 

Single and double shielded cables and cables with both 

pigtail and 360 degree backshell to shield connections were 

studied. The impulse tests were performed with an impulse 
generator connected both as a simple capacitor bank and as a high 
voltage Marx impulse generator. The drive current was a 230-kHz 
damped sinusoid with currents ranging from 200 to 23,000 amperes 
crest current. 

The basic problem and the solution may be summarized 

as follows. With the fan cowls which have only a few 

attachment points, the EEC and the wiring are in effect on the 
electromagnetic exterior of the aircraft aid therefore the extra 
25 to 45 db of shielding, which would be provided by the aircraft 
skin if the wiring were inside the aircraft or were under a 

peripherally bonded metal engine cowl, must be obtained for the 

exposed wiring on the engine. This can be done by either using 
double shielding or by adding peripheral bonding and improved 
shielding to a composite fan cowl. 

















selves to the most frequent case encountered 
on aircraft, l.e. that of the negative "cor¬ 
ona", schematize the overall process in the 
following way : 

- first phase (fig. la) - A conductor 
(C) possessing a point A is subjected to an 
electric field E for example measured at a 
reference point L. When the field reaches at 
A a threshold E Q , ionization by electron 
collision occurs in the high local field 
near, the point. The positive and negative 
charges thus created, separate under the 
influence of the field and form an electric 
dipole of moment ql. The positive charge 
follows the field towards the point A, the 
negative one being carried along towards the 
exterior of the point ; 

- intermediate phase (fig. lb) - The 
positive ions are captured by the conductor 
(C), and secondary emission leads to new 
electronic avalanches. Negative ions are 
created by attachment of electrons to elec¬ 
tronegative molecules in the low field region 
away from the point. These move towards the 
outside and lead to the increase in a nega¬ 
tive space charge. The process continues as 
long as the creation of electrons-ions pairs 
is maintained ; 

- final phase (fig. lc) - The negative 
space charge becomes strong enough to reduce 
the electric field at’ the point A to values 
lower than E 0 , 

This puts an end to the avalanche pro¬ 
cess. Positive charge is then totally cap¬ 
tured by the conductor (C) and negative space 
charge is still moving towards the outside at 
a speed u_E where (u- is the mobility of 
negative ions. 

Another current pulse may reappear when 
the negative charge drifts far enough from 
the point A for the field reduction effect to 

The whole process is illustrated by the 
help of two oscillograms of figure 2 rep¬ 
resenting the pulses associated with a corona 
discharge whose mean current i m is of 70 
pA ; the first and intermediate phases corre¬ 
spond to the zone of quick variation of the 
pulse (zone I, fig. 2a) ;the final phase 
being visualized by the slow variation of the 
latter. The oscillogram (b) provides the re¬ 
currence of these pulses which is here close 
to 500 kHz. 

product of the frequency f by the charge q 
liberated in each pulse. It is roughly rep¬ 
resented by a variation law of this form : 

im = K(F. - F. a )2 

where K is a constant at a given pressure, 
E s being the value of K at the threshold ot 
the corona discharge. 

From a simplified radiation emission 
standpoint the "corona" discharge can he as¬ 
similated to an electric dipole of moment ql, 
where q and 1 are variable parameters in time 
and whose fundamental frequency of excitation 


ELECTROMAGNETIC COUPLING BETWEEN AN 
ANTENNA AND AN ELECTRIC DISCHARGE - We have 
represented in figure 3 a schematic configur¬ 
ation of electromagnetic coupling between an 
antenna and an electric discharge : the el¬ 
ement (1) is a conductor possessing a prefer¬ 
ential zone of "corona" discharge occurrence 
(point A), the element (2) is the considered 
antenna. This configuration can be taken into 
account either for the evaluation of the 
electromagnetic characteristics of a dis¬ 
charger in the laboratory, or for the deter¬ 
mination of interference of a reception sys¬ 
tem on aircraft with or without potential 
dischargers ; in this last case the element 

(1) will represent the aircraft structure 
itself. 

The parasite short-circuit current of 
the antenna i cc ( . ) is given by [1] 

[2] : 

i C c () = L () - H LM ( ’> 


where :l(.) Is a coupling function rep¬ 
resenting the transmission between a refer¬ 
ence point L of the domain D of the corona 
discharge and the antenna. H L ^(■) is an 
independent function characterizing the 
properties of the noise source Inside the 
domain D and evaluated at the same point L. 

Tanner [1] has expressed these two func¬ 
tions in terms of the following two re¬ 
lations : 


( 2 ) 


e L ( ) Is the component of the radiation 
field normal at (1) at the point L where a 
voltage V( .) is applied between the conduc¬ 
tors (1) and (2) 


.(.) ■ 






j M ( )dD 


(3) 


Jm( ) is the current density vector of fre¬ 
quency for any point M of the discharge and 
F-h( • ) the corresponding field at M when the 
tension is applied between (1) and (2). 

An almost static condition applies in¬ 
side the domain 


..,( . ) 


= k, „ 


(M 


u« is unit vector tangent to Em( ) at 
point M, and k^x a scalar. 

Under these conditions (3) becomes : 


the 


(5) 


The function ) will be determined 
either for a configuration of tests in the 
laboratory or, on the aircraft, between a 
point L of the structure and the antenna 
being considered. Hlm( ) can be con¬ 
sidered as an Intrinsic function of the dis¬ 
charger to be evaluated in the laboratory or 
on the aircraft. 
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OPTIMIZATION OF THE DISCHARGER FUNCTION 
FROM THE ELECTROMAGNETIC RADIATION STAND 
POINT [2] - Optimization Implies minimizing 
isequently, the study 
Lows several possi- 


thogonal configuration 

c) to reduce to negligible values the 
frequency dependent components of the current 
j( ■)» which induces us to consider a continu¬ 
ous discharge regime ; 

d) to increase in the domain D the num¬ 
ber N of elementary discharges so that, for 
the efficiency of a givin discharge, the re¬ 
spective amplitudes of the N components 
jjj(-) can be reduced due to the fact that 

he corresponding pulses are not simul- 
aneous. 

Reduction of the scalar ” The 

principle of the method consists in moving 
the zone D from the discharge of the element 
(1) in order to reduce the inductive capaci¬ 
tive effect of the space charge in motion. On 
reconsidering the schematization of the coup¬ 
ling in figure 3b, we can examine the case in 
which the extremity of the conductor is iso¬ 
lated in such a way that a third conductor is 
created at the level of which the discharge 
appears. In order to maintain on (3) the 
presence of a "corona" discharge, it is of 
paramount importance to keep up between (1) 
and (3) the electrostatic potential ; this 
condition is achieved by inserting, between 
these two elements, an appropriate resistance 
R. 

If we compare again the "corona" dis¬ 
charge with an electric dipole f, we notice 
that the quantity of charges displaced by ca¬ 
pacitive induction on (1) is reduced, or, in 
other terms, for the same tension V( u ) ap¬ 
plied between (1) and (2), the electric field 
Ej|( ) is reduced. 

In order to guarantee the radiofrequency 
decoupling mentioned, the resistance Inserted 
between (1) and (3) can be optimized. For 
frequencies higher than 100 kHz, a value 


rent flows In the antenna short-circuit. 

Considering the decoupling resistanc 
between (1) and (3), the orthodecoupling i 
complementary to the decoupling by resistanc 
previously observed. 

Decoupling by independent mult id is 
charges [2] - The principle considered her 
does not require any outside intervention i 
the discharge mechanism ; it only consists i 
appropriating the emission zone so as to pro 


(3) (cf. fig. 3d) is made up 
Lve material and its config 
o produce discharges on a c 
a plane perpendicular to 


axis of the discharger. Pla 

field, we notice that the whole circumference 
of the edge shows a great number of dis¬ 
charges emitting a quasi-continuous total 
current, which implies, on the one hand, that 
each discharge emits a small quantity of 
charges and, on the other hand, that these 
discharges are decorrelated in time. 

This discharger presents therefore three 
complementary methods of decoupling : 

- a decoupling by resistance, 

- an orthodecoupling, 

- a decoupling by multidischarges decor- 
related in time. 

It Is this type of discharger which has 
been used in the tests performed in the lab¬ 
oratory and on aircraft and analysed in the 
two following paragraphs. 

The discharger is made up of a resistive 
cylinder whose total resistance between ex¬ 
tremities is of the order of 100 Mil ; the 
external part of the cylinder is organized, 
as shown in the sketch (a) of figure 4, so as 
to form an emission zone perpendicular to the 
discharger axis. The picture (b) in figure 4 
shows the envelope of corona discharges emit¬ 
ted while the discharger is in operation. 

DISCHARGER EFFICIENCY 

A discharger efficiency is characterized 


in relation to the radiofrequency 
a level with the discharge. The 
jsed is schematized in figure 3c : 
tor (3) remains isolated from the 
) owing to the optimal resistance 
; ; the conductor (3) is geometri- 
ited to have its emission point 
n such a way that the dipole f 
a direction perpendicular to the 
symmetry of (1). We notice that, 
;ly, charges of both signs dis- 
(1) by the inductive effect of the 
ire equal in absolute value ; their 


discharg 

of space charge are influenced by the 
trie field configuration of the experi 
set-up. 

LABORATORY CHARACTERIZATION - T 
perimental configuration which is mos 









quently used is that of point-plan arrange¬ 
ment (cf. "corona" discharge studies [5] 
[6]). We have used an electric field con¬ 
figuration around the discharger which aims 
at approximating closely the flight con¬ 
ditions (fig. 5) ; the discharger is in¬ 
stalled on a support representing a partial 
structure of the trailing edge ; this support 
is connected with the electric ground through 
a device for measuring the mean discharge 
current i d . The electric field is applied 
by means of an external electrode brought up 
to a potential V with respect to ground ; the 
polarization electrode is strapped to co¬ 
incide with an equlpotential surface, con¬ 
sidering the trailing edge geometry [1]. The 
electric field used as reference is measured 
at a point L of the support surface in front 
of the discharger location. 

This configuration allows to compare the 
discharger parameters defined previously with 
a "corona" point used as a reference or per¬ 
haps with other types of dischargers. This 
reference point whose length is similar to 
that of the discharger, possesses at its ex¬ 
tremity a known radius of curvature r and 
gives the characteristic i d =■ f(EL) cor¬ 
responding to curve (1) in figure 5. Its 
threshold is E s - 10 kVm“l and the 
evolution of i d vs E s is quasi quad¬ 
ratic ; this result is in agreement with the 
relation given by Loeb [5] 


, = K(E - E ) 


[2 - |E] 


( 6 ) 


where K is a coefficient related to the ge¬ 
ometry of the test and where d is the dis¬ 
tance between the point extremity and the 
polarization electrode. In our case 
8r/d « 1 therefore s 


= K(E - E ) 


(7) 


The physical justification of this vari¬ 
ation law is explained by the fact that the 
quantity of charge Q liberated in each 
Trichel pulse is proportional to the electric 
field existing at the point and that the fre¬ 
quency of these pulses 1/At is proportional 
to the drift velocity of the charges m_e 
( u- : mobility of negative ions) which is 
also proportional to the field, thus : 

{ d = At ' >E2 

We notice that in the configuration of 
tests which has been used, the field E very 
rapidly becomes small and constant outside 
the extremity of the discharger. 

In the configuration of tests which have 
been used, the threshold E s i for "corona" 
at the point used as a reference depends only 
on local pressure. 

Under the same experimental conditions 
we have shown in figure 5, curve 2, the vari¬ 
ations of i d vs E l for the discharger 
defined above. 

A threshold E s2 , measured on the 
sensor L, of 35 kV.m - ^, is higher than 
that of the reference point (10 kV.nf 1 ). 


This discrepancy is due to intensification of 
the field which is more important for the 
first case than for the ring edge of the dis¬ 
charger. 

After a first zone of quadratic vari¬ 
ation, i d tends towards a straight line 
starting at the functioning piint P. For the 
same electric field, the discharger current 
is lower than that of the reference point. 

Since the discharger has multidischarges 
not correlated in time, the space charge 
appears in the form of a continuous flood of 
negative ions driven by a drift velocity pE. 
The space charge field E p reduces the pro¬ 
duction of primary electrons at the dis¬ 
charger tip. From equation (7) : 

i d = K(E - E p - E g2 ) 2 

Ep being proportional to the current 
i d we have : 

i d = K(E - ui d - E g2 ) 2 


As soon as i d is sufficiently large we 
tend towards the linear relation : 

= a E (8) 

In the presence of a flow of air with a 
velocity u, the relation (8) can be written 
in the form : 


id = — E + euE 


(9) 


where £ is a constant [5]. 

AIRCRAFT OPERATIONS - a) Ground testing 

The aircraft used in a Meteor NF11. The 
instrumentation intented for the evaluation 
of potential dischargers has already been 
described in document [4]. 

The first series of measurements con¬ 
sisted in polarizing the aircraft, on the 
ground, by means of a high voltage generator 
and measuring the emission currents of a set 
of dischargers. These dischargers, here 21 of 
them, have been Installed at the wing ex¬ 
tremity (2x7), at the elevator extremity (2 
x 3) and at the top of the rudder (1) ; the 
dischargers are spaced 30 cm apart. 

The electric field used as a reference 
is measured at a point A of the structure on 
the aircraft fuselage. 

In figure 6 we show the current vari¬ 
ation for the discharger at the top of the 
rudder i d H vs E A . Beyond a threshold 
E s of the order of 3.3 kV.m -1 , we ob¬ 

serve the expected variation with E A : a 
quasi-quadratic variation up to the point P, 
followed by a linear variation. We must point 
out a very good reproductibility of measure¬ 
ments in all these tests. Knowing that the 
threshold Eg measured in the laboratory and 
on the aircraft are not directly comparable 
we have indicated in figure 7 the variations 
of i d vs E/Eg for both tests ; we can no¬ 
tice between the two types of measurements a 
very important discrepancy. If we compare the 
slopes of the linear characteristics, we ob¬ 
tain respectively 1/> * 31.76 pA/kV.m”' 

for the laboratory tests and 1/ < • 6.90 









fiA/kV.m - *’ for the aircraft tests, i.e. a 
ratio 4.6. This difference is due to the be¬ 
haviour of the space charge created by the 
emitted ''corona” current ; indeed, both test 
series being performed at atmosphere pressure 
and the field conditions around the emission 
zone of the discharger being identical, only 
the distribution of electric fields outside 
the discharger are different ; consequently 
the drift velocity p E of negative ions 
changes and therefore, the space charge den¬ 
sities. We reach very rapidly a constant 
field in the case of the laboratory configur¬ 
ation (not far from the values measured on 
the discharger support) whereas it tends 
towards zero in the case of the polarized 
aircraft. It is therefore the term a of the 
relation (8) which practically determines the 
value of the emission current as soon as it 
becomes Important. This result is qualitat¬ 
ively corroborated if we examine the charac¬ 
teristics of the dischargers at the wing ex¬ 
tremity (dV) and at the elevator extremity 
(dE) which have an operating threshold ident¬ 
ical to dH for the configuration of the air¬ 
craft chosen (3.3 kV.m - ^ measured at A). 

The corresponding characteristics i<j - 
f(E/E s ) as shown in figure 7 are slightly 
different from those of the discharger dH. At 
last, a final result allows to confirm the 
essential part played by the space charge in 
the discharge efficiency of a discharger ; if 
we compare the operating conditions of the 
seven dischargers placed on the aircraft 
wings, we notice that the discharger in the 
middle systematically emits the smallest 
quantity of current and gives, whatever the 
field E A may be, a discrepancy of the order 
of 60 X in relation to the two dischargers 
nearby. This median discharger is therefore 
essentially Influenced by the space charge 
liberated by the other six. 

b) Tests in flight - Electric field ref¬ 
erence - Measurement of the emission currents 
associated with the dischargers does not pro¬ 
vide any difficulty because of specific in¬ 
strumentation attached to each discharger. 
However, the determination of an electric 
field reference is proving to be much more 
complicated. Indeed, an aircraft in flight 
may be affected by the effects of its own 
electric charge Q due to triboelectric im¬ 
pacts or by the charge due to an external 
field E for example, the local atmospheric 
electric field. 

We know that if the external electric 
field may be considered as uniform in a vol¬ 
ume greater than that of the aircraft, the 
evaluation of its own potential V “ Q/C, (C 
is the aircraft capacity) and the three 
components of the field E may be obtained 
from 4 measurements of the electric field 
Ef on the structure [4]. In these con¬ 


ditions we have : 



where M is a matrix with g ^ coefficients 
corresponding to the location of the sensors 
A, B, C and D on the structure. 

In the particular case of an external 
field equal to zero, the expression (10) is 
reduced to : 

E. = B.V (ID 

The first test we are going to analyse 
deals with a flight (No 1) in which the air¬ 
craft does not possess potential dis¬ 
chargers ; this flight was performed at high 
altitude in clouds made of ice particles and 
without any external field. 

In figure 8 we have drawn two curves 
corresponding to the variations of the impact 
current SI measured on a capture electrode on 
the aircraft structure and those of the field 
E A vs time. 

We notice a strong correlation between 
the time evolution of both parameters. The 
values E A » 0 are observed simultaneously 
with the values SI ■ 0 tfhus clearly indi¬ 
cating the absence of any noticeable external 
field. For the whole sequence we have shown 
in figure 9 the correlation between the elec¬ 
tric field measured at the point A and the 
value of the potential V of the aircraft cal¬ 
culated by means of the relation (10) ; we 
clearly observe (cf. curve flight 1) a lin¬ 
earity between both terms, in agreement with 
the relation (11) ; however the slope 1/B A - 
2.783 m of the straight line is rather dif¬ 
ferent from the theoretical value 1/8 ATH - 
4.149 m. A 

The test having been performed without 
an external field we have to take into ac¬ 
count the space charge effects. The fact that 
the field E A is proportional to the poten¬ 
tial V shows that the aircraft charge is re¬ 
lated to charging and emission (discharging) 
currents of the aircraft. 

Indeed, these space charges can have two 
origins : 

- the former corresponds to positive 
particles liberated after collision with the 
front structures of the aircraft ; 

- the latter is related to the negative 
ions emitted by discharges appearing on the 
structure of the aircraft. 

Since the aircraft is negatively charged 
its distribution of charge can be schematized 
as shown in figure 10. The Impact surfaces 
are acted upon by positive external charge 
produced in the collisions. It then follows 
that a part of the negative charge of the 
aircraft is drawn into these areas ; the 
trailing edge zones, affected by "corona" 
discharges, are subjected to the effects of 
the negative external space charge which 
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reduce the discharge efficiency of 


ure 9 the curves V = f(E A ) obtained during 
two other flights without a discharger (Nos 2 
and 3). We can notice slopes 1/2 A slightly 
different from the previous one. For high po¬ 
tentials a change in the slope of the curve 
V = f(E A ) may be observed ; this is the 
case of flight No 2 in which potentials ex¬ 
ceeding 300 kV were obtained. As shown in the 
curve of figure 9 for flight No 2 a change in 
the slope (point I) is situated around E A = 
65 kV.nf* ; this modification may be due 
to the occurrence of new "corona" discharges 
in different areas of the structure therefore 
to a new redistribution of charges on the 
aircraft. 

Functioning of the dischargers - Three 
flights are analysed that correspond to three 
different installations of the dischargers. 
Their main characteristics are summarized in 
the table below ; 

Flight Number location aircraft local 
No of dis- of dis- velocity pressure 
chargers chargers | 




HD : top of the rudder. 

The curves V = f(E A ), obtained for 
these three flights, exhibit the same charac¬ 
teristics as those obtained in situations 
without dischargers. We can observe again 
(figure 9) linear evolutions of the potential 
vs E a for flights Nos 4 and 5, in which 
high values of the potential have been re¬ 
corded. Changes in the slopes (point J) take 
place at a field value of 42 kV.m" 1 . Let 
us notice, however, that the slopes 1/: A of 
these curves are lower than those obtained 


emitted by the whole set ot dis- 
in figure 11 we have drawn three 
ing the variations of If vs E A 
Nos 4, 5 and 6. 

i notice that the evolutions of 
Lso linear in E A but with a 
slope (points I) in all three 
change of slope detected during 
flight No 5 corresponds to the value 1 
42 kV.m~l mentioned above. 

The general repartition of the cur 
curves Ij = f(E A ) in figure 11 is sim 
to that for the corresponding voltage cu 
V = f(E A ) in figure 9. Except that the 
charge thresholds are not approximatively 
lated to the values of the local pressure 
the three flights considered. Whereas the 
haviour of the discharger installations in 
not notably different at low currents, dis¬ 
charge characteristics tend to vary a good 
deal at strong currents. For example, com¬ 
paring flights Nos 4 and 5, we notice a cur¬ 
rent value of If = 350 pA corresponding to 
field values E A of 50 and 107 kV/m (ratio 
2). Since for the three tests, the velocities 
of the aircraft are close to one another this 
dispersion is essentially due to the space 
charge created by the dischargers whose num¬ 
ber varies (23, 12 and 9 dischargers) accord¬ 
ing to the zones of the structure. 

Figure 12 allows to illustrate the 
specific performance of a discharger namely 
the one dH, at the top of the rudder. The 
starting thresholds are different from those 
to be expected by the changes in local press¬ 
ure particularly for the flights Nos 4 and 6. 
It must be mentioned that the threshold cor¬ 
responding to the flight No 5 correlates with 
the value E A of 42 kV.m - ^ already ob¬ 
served, which would mean that the changes in 
the slope noticed in the characteristics V » 
f (E a ) or I.j. » f(E A ) could be related to 
the thresholds of the discharges at the top 
of the rudder ; besides this point coincides 
also with a change in the slope of the 
characteristic i^jj for the flight No 4. 

Another important fact is to be men¬ 
tioned : it deals with hysterisis of the 
characteristics ij “ f (E) observed during 
several tests ; we notice that when the cur¬ 
rent decreases it exhibits a threshold field 
value that is lower than the one for its in¬ 
itial onset. 

The curve in figure 13 was obtained in 
another sequence of measurements in flight 
No 4. It clearly illustrates this hysterisis 















for decreasing currents the operating point 
follows at first the same straigth character¬ 
istic, but eventually deviates from it and 
leads to lower threshold field value of 17.5 
kV.m -1 . 

The interpretation of this phenomenon 
must also be explained by the effect of space 
charge that lowers the discharge onset of the 
dischargers located at the top of the rudder. 
The space charge delay the starting of the 
discharger but modify its discharge ef¬ 
ficiency until a new electrostatic balance of 
the aircraft is obtained. 

We can sum up the characteristics of 
this top of the rudder discharger by means of 
different curves of figure 14 in which we 
have described the emitted current ijjj in 
relation to the reduced field E A /Eg, Eg 
being the threshold field measured at the 
point A for each test. 

- When the number of dischargers is 
large (23 for the flight No 4) and for fields 
larger than Eg, we observe a character¬ 
istic current i - f(E A /Eg) close to that 
obtained on the ground tests. This result 
shows an electrostatic configuration around 
the aircraft which is nearly the same in both 
cases. The greater observed efficiency is due 
to the aircraft velocity. 

- For a small number of dischargers (9 
for flight No 6) the slope of the linear dis¬ 
charge characteristic is close to that ob¬ 
tained in the laboratory, which indicates a 
relatively homogeneous field in areas remote 
from the discharger tips. 

- For intermediate cases (flight No 5 
with 12 dischargers) the discharge character¬ 
istic of the discharger will be located be¬ 
tween the two extreme operating configur¬ 
ations which have just been defined. 

RESIDUAL ELECTROMAGNETIC RADIATION 

PERFORMANCE OBTAINED IN THE LABORATORY - 
We have shown that the current emitted by a 
"corona" discharge is related to the move¬ 
ments of positive and negative charge that 
give rise to very brief pulses (of the order 
of 5 ns) and whose repetition frequency in¬ 
creases according to the intensity of the ap¬ 
plied electric field. The radiation noise as¬ 
sociated with emission at the reference point 
defined previously is analysed in relation to 
frequency by using the setting up described 
in figure 5. Here the discharger and its sup¬ 
port can be considered as the reception 
antenna itself. The spectral distribution of 
noise is shown in picture (a) of figure 15 : 
we see a frequency distribution corresponding 
to the repetition frequency of current 
pulses. Experimentally and for average cur¬ 
rents exceeding 10 pA, noise can be detected 
up to 100 MHz. 

During the same laboratory tests we have 
shown that the decrease in radiation noise 
obtained through a discharger with coupling 
resistance is of the order of -30 dB and 


rises to about -50 dB if an orthodecoupling 
is added. In the case of a discharger with 
multidischarges a white noise is practically 
obtained and we can state that the reduction 
(by comparison with a "corona" point oper¬ 
ating with a similar mean current) goes to 
-80 dB for frequencies ranging from some 10 
Hz to several MHz. 

AIRCRAFT TESTS - Radio Electric Noise at 
350 kHz - The receiver used is a classical 
ADF system operating as a detector with dy¬ 
namic compression. The reception frequency is 
350 kHz and the bandwidth is 3 kHz. 

The aim of the test is to determine the 
interference level affecting the receiver in 
the presence of a "corona" discharge that is 
well defined and detected (reference point). 
The aircraft is isolated from ground by insu¬ 
lating blocks. It is connected to the earth 
by a resistance of 30 Mil in order to create 
coupling conditions similar to those in 
flight, between the discharger and the corre¬ 
sponding antenna. 

The reference point is at the wing tip 
and a suitable polarization device allows to 
vary the intensity of the average current 
emitted by this point. The radiation noise 
measured with respect to the mean current is 
represented in figure 15. The coupling func¬ 
tion between the location of the point and 
the antenna is small and the noise level is 
rather low. No clearly established variation 
can be observed. This demonstrates the impul- 
sional nature of the discharge ; noise maxima 
are recorded when an harmonic of the rep¬ 
etition frequency of current pulses is situ¬ 
ated in the bandwidth of the receivers. This 
point is ascertained by a discrepancy ob¬ 
served, for the same emission current, depen¬ 
ding on whether the latter is obtained in an 
increasing phase or a decreasing one ; the 
exact harmonic analysis of pulses shows in¬ 
deed a frequential discrepancy between both 
experimental methods. 

Discharger Effects on Aircraft - In or¬ 
der to measure the performance of the dis¬ 
chargers from the point of view of radiation 
we can consider again the three flights (Nos 
4, 5 and 6) performed with installation of 
23, 12 and 9 dischargers, and we compare the 
results obtained with reference to flight (No 
1) performed without a discharger. 

The radiation noise observed corresponds 
to the levels recorded in the ADF system at 
350 kHz and whose detection threshold is of 
the order of some 10“^ V in the case of 
the aircraft in flight without electrifi¬ 
cation phenomenon. 

Flight No 1 (fig. 17) - For the flight 
used as a reference, an important increase in 
noise is observed at a field of the order of 
25 kV.m -1 . This coincides with the 

threshold for the wing or elevator ex¬ 
tremities ; the noise reaches a maximum of 
300 pV for E a = 40 kV.m -1 and then re¬ 

sumes its equilibrium at a value 100 pV. The 
highest level detected in the low field re- 
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glon may be explained by the fact the rep¬ 
etition frequency of the various discharges 
appearing on the structure are located In the 
area of the ADF receiver, whereas for higher 
fields only harmonics of frequency rep¬ 
etitions are taken into account. 

Flight No 6 - For E A ^ 25 kV.m -1 , 

the discrepancy observed In relation to 
flight No 1 is not significant ; the noise 
below 25 kV.m -1 Is undist inguishable from 
flight No 1. But, on the other hand, a very 
rapid evolution of Interference at 350 kHz 
occurs at 35 kV.m - , explaining the ap¬ 
pearance of one or several "corona" dis¬ 
charges strongly coupled with the antenna ; 
the areas taken Into account may be the 
trailing edges or the top of the rudder. A 
rapidly decreasing interference depending on 
the field, clearly Indicates the presence of 
a "corona" source that is very localized and 
strongly coupled with the ADF antenna. 
Because of the low level of electrostatic ac¬ 
tivity observed during this flight, we do not 
have any data for values of E A higher than 
41 kV.m -1 . 

Flight No 5 - The Installation of an ex¬ 
tra discharger on each elevator and at the 
top of the rudder allows to suppress the 
overall interference appearing around the 
value E A ■ 40 kV.m -1 . We can observe a 

noticeable increase in noise for higher 
values of the field around 60 kV.m -1 . At 
this level the 12 installed dischargers have 
almost identical operating conditions. This 
indicates that the unprotected discharges ap¬ 
pear in the strong electric field areas that 
require an increase in the number of dis¬ 
chargers. 

Flight No 4 - Three extra dischargers 
have been installed on the wings, two on the 
elevator and one at the top of the rudder. No 
increase in noise can be noticed up to values 
E a reaching 75 kV.m -1 . This demon¬ 

strates nearly complete radiation protec¬ 
tion ; the recorded noise level is very close 
to the detection threshold of the receiver. 

The Important conclusions of these tests 
are expressed on the following statements : 

a) Noticeable radiation interference ap¬ 
pears in the detection system only it an un¬ 
protected "corona" discharge is created on 
the aircraft structure. This result is ascer¬ 
tained by tests performed on the ground where 
the establishment of a single "corona" dis¬ 
charge is sufficient to disturb the ADF re¬ 
ceiver. 

b) In strong electric field areas, the 
number of dischargers must be sufficient to 
avoid the formation of discharges directly on 
the structure. We are of the opinion that we 
must not exceed a maximal current of 50 )oA 
for each discharger. 

c) The performance of a discharger is 
sufficient if the decrease in radio noise for 
a current of 50 pA, reaches -50 dBm when com¬ 
pared to an unprotected discharge of the same 
mean current. 


CONCLUSION 

The following conclusions can be drawn 
from this analysis : 

a) the electrostatic configuration of an 
aircraft in flight in the presence of tribo- 
electric impact is very complicated ; in ad¬ 
dition to the aircraft potential due to nega¬ 
tive charges accumulated on the structure, 
Important local electric fields due to par¬ 
ticles charged after collision and negative 
ions liberated by the "corona" discharges 
have to be evaluated ; 

b) this electrostatic configuration mod¬ 
ify the discharge characteristics of the air¬ 
craft ; we show that they are linked to the 
presence of space charge which raise their 
operating threshold and make the evolution of 
the emitted current linear vs the potential 
or vs the electric fields on the structure ; 
hysteresis phenomenon in the functioning of 
the dischargers are frequently observed ; 

c) electromagnetic radiation inter¬ 
ference may be practically suppressed on an 
aircraft if the number of the dischargers is 
sufficient ; from the maximal value of the 
emission current for an aircraft of a given 
size and speed [1] the number of dischargers 
may be evaluated by limiting the maximal 
value of current per discharger to be of the 
order of 50 pA ; 

d) the dischargers must exhibit a radi¬ 
ation coupling of at least -50 dBm in re¬ 
lation to an unprotected "corona" discharge 
of the same mean current. 
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